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Abstract  of  Thesis  Presented  to  the  Graduate  Council 
of  the  University  of  Florida  in  Partial  Fulfillment  of 
the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 

TEMPORAL  AND  SPATIAL  TEMPERATURE  VARIATIONS 
IN  UNHEATED  AND  SOLARLY  HEATED  SWIMMING  POOLS 

by 

John  Theodorou  Triandafyllis 
March,  1977 

Chairman:  Erich  A.  Farber 

Major  Department:  Mechanical  Engineering 

The  purpose  of  this  investigation  was  to  determine 
analytically  and  experimentally  the  spatial  water  tempera- 
ture profile  and  the  time  history  of  the  bulk  water 
temperature  in  an  unheated  swimming  pool  and  in  a swimming 
pool  which  is  heated  by  solar  flat  plate  collectors. 

The  water  temperatures  were  measured  in  two  identi- 
cal circular  swimming  pools.  One  pool  was  kept  uncovered 
and  was  unheated  by  any  external  means.  The  other  pool 
v?as  covered  with  various  transparent  plastic  covers,  and, 
in  addition,  it  was  heated  by  flat  plate  collectors. 
Measurements  in  both  pools  were  taken  for  a period  of 
twelve  months.  The  experimental  data  were  compared  and 
were  found  to  be  in  excellent  agreement  with  the  water 
temperatures  as  were  calculated  by  various  mathematical 
models . 
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The  spatial  temperature  variations  were  obtained 
analytically  by  solving  the  continuity,  momentum,  and 
energy  equations  of  an  infinitesimal  volume  of  water,  and 
then  integrating  over  the  whole  volume  of  the  water. 

The  boundary  conditions  were  of  a time  changing  solar 
radiation  source , and  time  changing  amounts  of  energy , 
which  were  transferred  from  the  water  to  its  environment. 

A fast  converging  series  solution  was  obtained  using  Finite 
Cosine  Fourier  Transforms. 

Analytical  considerations  indicated  that  the  motion 
of  the  water  of  the  swimming  pool  approximated  a torroid. 
The  torroid  during  the  daytime  moved  in  opposite  direction 
to  the  torroid  during  the  nighttime.  The  motion  of  the 
water  was  induced  by  free  convection  currents  which  were 
generated  from  the  heating  or  cooling  of  the  water  next 
to  the  swimming  pool  walls.  The  analytical  and  experi- 
mental results  of  this  investigation  confirmed  the  pres- 
ence of  a flow  regime  in  which  (a)  there  is  zero  vorticity, 
(b)  there  is  a vertical  temperature  gradient,  and  (c)  the 
water  temperature  is  almost  constant  in  the  horizontal 
direction.  The  water  moved  with  a maximum  velocity  next 
to  the  wall;  this  velocity  was  calculated  analytically 
and  was  verified  with  other  free  convection  mathematical 
models  from  the  literature.  In  the  derivation  of  the 
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velocity  equation,  the  effects  of  viscous  dissipation  as 
compared  to  the  buoyancy  forces  were  considered  negligible. 

The  time  variations  of  the  bulk  water  temperature  were 
obtained  from  an  overall  accounting  of  all  energy  fluxes 
which  cross  the  boundaries  of  the  control  volume  of  the 
swimming  pool.  The  time  histories  of  the  bulk  temperatures 
were  calculated  and  were  verified  experimentally  in  three 
time  frames:  annual,  monthly,  and  daily.  The  analysis 

and  the  experiments  indicated  that,  in  Gainesville,  Florida, 
the  water  of  a swimming  pool  which  was  heated  by  flat  plate 
collectors  was  kept  over  80°F  for  a period  of  twelve 
months.  The  required  area  of  the  flat  plate  collectors 
for  the  winter  was  equal  to  1.5  times  the  surface  area  of 

i 

the  pool. 

Another  method  of  raising  the  water  temperature  of 
the  swimming  pool  was  by  covering  the  top  surface,  thus 
eliminating  evaporation  losses.  This  method  extended 
the  swimming  season  by  almost  three  months.  Various 
plastic  covers  and  different  numbers  of  layers  were 
tested . 
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CHAPTER  I 


INTRODUCTION 

In  recent  years,  man  has  become  painfully  aware  that 
the  supply  of  the  earth's  fossil  fuels  is  not  inexhausti- 
ble. The  nuclear  fission  power  plant,  man's  greatest 
hope,  is  still  facing  serious  technical  design  problems. 

It  is  no  wonder , then , that  people  have  started  consider- 
ing the  utilization  of  solar  energy  as  an  alternative 
source  of  energy.  The  solar  radiation  falling  on  earth 
is  free,  abundant,  and  non-polluting. 

In  recent  years  scientists  have  conducted  a great 
amount  of  research  in  solar  energy  and  have  built  the 
foundations  for  the  economical  utilization  of  solar  energy 
to  man.  They  have  demonstrated  the  utilization  of  solar 
energy  in  heating  and  air  conditioning  buildings,  dis- 
tilling and  recycling  water,  operating  solar  engines, 
furnishing  power  for  spacecraft  operations,  and  in 
many  other  applications  [1,  2,  3]. 

One  application  of  solar  energy  is  in  the  heating 
of  swimming  pools . Even  in  the  southern  parts  of  the 
United  States,  there  is  a need  for  swimming  pool  heating. 
Although  the  sun  is  the  natural  heating  agent,  there 
are  high  heat  losses  from  the  water  that  prevent  swimming 
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for  many  months  during  the  year.  If  one  were  able  to 
reduce  the  heat  losses  and/or  to  supplement  the  heating 
effect  by  the  sun,  then  the  swimming  season  could  be 
extended.  Heating  the  water  in  a swimming  pool  by  electric 
resistance  heaters,  oil,  and  natural  gas  is  for  the 
average  man  an  expensive  proposition.  New  energy  sources, 
such  as  solar,  must  be  employed  to  reduce  the  cost  of 
heating  swimming  pools . 


CHAPTER  II 


PREVIOUS  INVESTIGATIONS 

There  are  two  methods  of  raising  the  water  tempera- 
ture of  a swimming  pool  without  the  direct  or  indirect  use 
of  fossil  fuels:  (1)  By  reducing  the  heat  losses  from  the 

water  of  the  swimming  pool  to  its  environment,  and  (2) 

By  supplementing  the  heating  effect  of  the  sun  by  using 
external  energy  sources.  A literature  survey  has  shown 
that  scientists  have  researched  both  methods. 

The  heat  losses  from  the  watdr  are  due  primarily  to 
evaporation,  convection,  and  night  sky  radiation  of  the 
water  from  the  top  surface  of  the  swimming  pool.  The 
easiest  and  most  economical  way  to  reduce  the  losses  due 
to  evaporation  and  convection  is  to  cover  the  top  surface 
of  the  water  with  a material  which  will  allow  the  sun  to 
pass  through  it  and  is  impervious  to  water.  Brooks  [4] 
in  1955,  suggested  using  a transparent  plastic  sheet. 
Robson  [5]  preferred  using  a black  plastic  sheet  to 
maximize  the  heating  effects  of  the  sun  (black  paint  us- 
ually has  a high  absorptivity  value  for  solar  radiation) . 
He  also  suggested  lining  the  bottom  side  and  the  walls 
of  the  swimming  pool  with  a black  plastic.  The  black 
lining  maximizes  the  amount  of  solar  energy  which  is 
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absorbed  by  the  bottom  and  the  walls.  Robson  further 
compared  the  temperature  rise  of  the  water  in  a swimming 
pool  covered  with  a black  sheet  and  one  which  was  uncovered 
but  had  its  bottom  side  and  the  walls  lined  with  a black 
sheet.  He  found  a much  greater  temperature  rise  in  the 
covered  pool.  This  suggests  that  the  heat  lost  due  to 
evaporation  is  far  greater  than  the  heat  gained  by  the  in- 
creased absorptivity  of  the  walls  and  bottom  side  of  the 
swimming  pool.  Lof  [6]  used  a cover  made  of  mylar  plastic. 
Root  [7]  suggested  using,  instead 'of  a plastic  cover,  an 
oil  film  floating  on  the  water.  Czarnecki  [8]  inflated 
transparent  plastic  covers  with  air  and  floated  them  on 
the  water.  The  trapped  air  reduces  the  overall  heat 
transfer  coefficient  from  the  top  surface  of  the  water 
to  its  environment,  thus  reducing  the  convection  losses. 

In  addition,  the  evaporation  losses  have  been  minimized 
due  to  the  floating  of  the  plastic  covers  on  the  water 
surface . 

To  reduce  the  night  sky  radiation  losses  from  the 
water,  Robson  [5]  suggested  using  a plastic  cover  on  the 
water  which  has  a low  emissivity  value,  such  as  a 
silver,  or  aluminized  plastic  surface. 

The  method  of  supplementing  the  heating  effect  of 
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the  sun  by  using  external  energy  sources  is  to  circulate 
the  water  of  the  swimming  pool  through  a solar  flat  plate 
collector.  A flat  plate  collector  consists  of  a flat 
metallic  surface  which  is  coated  with  a paint  highly 
absorbent  of  sunlight.  The  metallic  surface  is  enclosed 
in  a frame  which  reduces  the  heat  losses  from  the  metallic 
surface  to  its  environment,  while  at  the  same  time  allowing 
solar  radiation  to  enter  the  collector  through  a "window." 
Part  of  the  solar  radiation  is  absorbed  by  the  metallic 
surface  and  its  temperature  is  raised.  If  the  water  of 
the  swimming  pool  comes  in  direct  contact  with  the 
metallic  surface,  then  heat  is  transferred  from  the  hot 
metallic  surface  to  the  water.  In  1956  Farber  [9]  designed 
a flat  plate  collector  which  consisted  of  a black  painted 
galvanized  sheet  wrapped  in  plastic.  As  an  added 
benefit  he  suggested  placing  the  flat  plate  collectors  to 
form  the  fence  around  the  swimming  pool.  Andrassy  [10] 
constructed  a flat  plate  collector  which  consisted  of 
corrugated  blackened  aluminum  sheets  with  plastic  tubing. 
These  collectors  extended  the  swimming  season  from  50 
days  to  152  days  (with  the  water  at  a minimum  temperature 
of  70 °F) . Another  type  of  a flat  plate  collector  was 
designed  by  Yellott  [11]  which  consisted  of  water  flowing 
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down  a blackened  masonite  backing;  a Tedlar  film  covered 
the  water.  He  measured  a temperature  rise  of  7°F  across 
the  collector  at  an  average  insolation  rate  of  250  BTU/(hr. 
sq.  ft.),  1 gallon/minute  flowrate  from  a four  foot  square 
collector.  Czarnecki [12]  suggested  using  the  corrugated 
iron  roof  of  a building  as  the  absorber  plate  of  a flat 
plate  collector.  The  corrugations  serve  as  conduits  for 
the  water  to  be  heated.  He  also  suggested  using  both  meth- 
ods to  raise  the  water  temperature  simultaneously — plastic 
covers  inflated  with  air,  and  a corrugated  iron  roof.  The 
swimming  season  could  be  extended  by  six  weeks  at  each 
end  of  the  season  in  Melbourne,  Australia.  Thomason 
[13]  designed  a plastic-sand  flat  plate  collector.  Water 
from  the  swimming  pool  flows  by  gravity  through  a layer 
of  sand.  The  layer  of  sand  is  sandwiched  in  a transparent 
plastic  bag  which  rests  on  the  roof  of  the  building. 
Whillier  [14]  concluded  that  in  order  to  raise  the  water 
temperature  10 °F,  one  needs  flat  plate  collectors  of  an 
area  equal  to  the  surface  of  the  swimming  pool. 

It  is  possible  to  gain  significant  knowledge  about 
tiie  heating  of  swimming  pools  with  solar  energy  by  studying 
^-he  flow  patterns  of  the  water  currents  which  occur  during 
the  heating  of  the  water  by  the  sun,  or  the  cooling  of  the 
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water  by  many  heat  transfer  processes,  such  as  natural 
convection.  Investigators  have  studied  these  water  cur- 
rents which  can  be  attributed  to  natural  convection,  and 
will  be  referred  henceforth  as  natural  convection  currents. 
Foster  [15]  performed  experiments  in  layers  of  fluid  in 
which  the  bottom  surface  temperature  was  increased  at 
a constant  rate.  The  side  walls  were  insulated.  He 
observed  the  formation  of  natural  convection  cells  which 
varied  in  horizontal  spacing  in  proportion  to  the  1/5 
power  of  the  rate  of  the  bottom  surface  temperature  increase 
Foster  [16]  also  developed  a mathematical  model  which 
describes  the  formation  of  natural  convection  currents 
in  the  upper  layers  of  the  ocean.  The  solar  radiation 
was  treated  as  a time-dependent  heat  source.  The 
theoretical  calculations  predict  the  formation  of  natural 
convection  currents  only  during  the  night,  provided  that 
the  temperature  of  the  air  is  lower  than  the  temperature 
of  the  water.  No  experimental  verification  of  his 
theoretical  results  was  presented. 

Eckert  and  Carlson  [17]  performed  experiments 
in  determining  the  temperature  variations  in  an  air 
layer  enclosed  between  two  vertical  walls  with  different 
uniform  temperatures.  The  horizontal  walls  were  adiabatic 
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to  heat  and  mass  flow.  They  observed  two  distinct  flow 
regimes.  In  one  regime  heat  is  transferred  by  conduction 
in  the  central  part  of  the  air  layer;  convection  contrib- 
utes only  in  the  corner  regions.  In  the  second  flow 
regime  boundary  layers  exist  along  the  surfaces  of  the 
enclosure  whereas  in  the  central  part  there  is  a constant 
vertical  temperature  gradient  and  the  temperature  is 
constant  in  the  horizontal  direction.  There  is  also  zero 
vorticity  in  the  central  part  of  the  air  layer.  Eckert 
and  Carlson  [17] , found  that  the  critical  dimensionless 
parameter  is  the  Rayleigh  number  (Ra) , and  that  for  an 
aspect  ratio  of  1,  a value  of  the  Rayleigh  number  of  500 
is  the  upper  limit  of  the  first  flow  regime.  Elder 
[18,  19]  calculated  and  measured  the  critical  Rayleigh 
number  to  be  1000.  Gill  [20]  calculated  the  temperature 
field  in  the  second  flow  regime.  Brooks  and  Ostrach  [21] 
performed  an  experimental  investigation  in  a horizontal 
cylinder  where  they  supplied  the  heating  at  various 
circumferential  locations. 


CHAPTER  III 


PURPOSE  AND  SCOPE  OF  PRESENT  INVESTIGATION 

The  literature  survey  indicates  a considerable  lack 
of  information  of  the  temperature  history  of  the  water 
in  a swimming  pool . The  temperature  history  of  the  water 
and  the  equations  which  predict  the  temperatures  can  be- 
come significant  design  aids  to  an  engineer  for  a variety 
of  applications.  These  applications  include  the  heating 
of  swimming  pools  on  a year  round  basis,  and  the  control 
of  the  temperatures  in  a pond  where  fish  are  raised,  or 
in  a pond  used  as  a heat  source  for  heat  pumps. 

The  purpose  of  the  present  investigation  is  to  pro- 
vide the  temperature  history  of  the  water  in  a swimming 
pool  and  to  gain  understanding  of  the  physical  phenomena 
occurring  during  the  heating  or  cooling  of  the  water. 
Further,  to  evaluate  the  heating  of  swimming  pools  by 
solar  energy,  the  temperature  history  of  the  water  in 
a solarly  heated  swimming  pool  will  be  compared  to  the 
temperature  history  in  an  unheated  pool  of  the  same 
design. 

It  is  intended  to  theoretically  calculdate  and 
experimentally  verify: 

(1)  The  hourly  variations  of  the  bulk  water 
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temperature . 

(2)  The  monthly  variations  of  the  bulk  water  tempera- 
ture . 

(3)  The  annual  variations  of  the  bulk  water  tempera- 


ture . 

(4)  The  spatial  variations  of  the  water  temperature 


at  any  instant. 


CHAPTER  IV 


THEORY 


Development  of  Equation 
for  Spatial  Temperature  Variations 
One  of  the  objectives  of  this  investigation  is  to 
calculate  the  temperature  variations  in  a swimming  pool 
as  a function  of  time  and  depth  of  water.  To  accomplish 
this,  it  is  necessary  to  develop  the  continuity,  momentum, 
and  energy  equations  for  an  infinitesimal  volume  of  water, 
and  then  to  integrate  them  over  the  whole  volume  of  the 
water . 

Considering  an  infinitesimal  volume  of  water  in  a 
cylindrical  swimming  pool  (Fig.  4-1) , the  continuity 
equation  may  be  written  as 


3v  v . 3v , . 3v 

_itL  + 

3r  r r 9<f>'  3z 


(4-1) 


where  v^,  v^ , and  v^  are  the  velocity  components  of  the 
water  velocity  vector  in  the  cylindrical  coordinates  r, 
4> ' , ^md  z respectively.  In  the  continuity  equation 
(4-1)  it  was  assumed  that  water  is  incompressible.  The 
continuity  equation  can  be  further  simplified  by 
considering  the  axial  symmetry  of  the  swimming  pool; 
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Figure  4-1 


An  infinitesimal  volume  in  a cylindrical 
swimming  pool . 
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therefore,  all  derivatives  with  respect  to  <J>'  are  zero. 

Then  the  continuity  equation  (4-1)  becomes 


3v  v 3v 

— - + — + — - 
3r  r 3z 


0. 


(4-2) 


The  momentum  equation  for  an  infinitesimal  volume  of 
water  is  the  Navier-Stokes  equation  in  the  axial  direction 
z as  was  developed  in  ref.  [22] . It  can  be  rewritten 
as 


3v  3v  3v 

z , z , z 

+ v + v 

3t  r 3r  z 3z 
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r 

3pl 

r32v 

z 

1 

3v  3 2v 

z z 

p 

F - 
z 

ifj+v 

3r  2 

+ — 
r 

3r  + 3 z 2 

J 

(4-3) 


In  the  momentum  equation  (4-3)  it  was  assumed  that  water  is 
incompressible,  and  the  axial  symmetry  was  incorporated, 
the  only  body  force  acting  on  the  water  is  the  gravita- 
tional force;  therefore,  we  have 


Fz  = -pg  (4-4) 

2 

where  g is  the  acceleration  of  gravity  = 32.17  ft/sec  . 
Since  the  gravitational  force  is  the  only  external 
force  and  it  is  acting  in  the  z direction,  it  is  evident 


that. 
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9p  _ d£ 
8z  dz  ’ 


. (4-5) 


Another  consideration  that  can  be  made  is  that  in  hydro- 
static equilibrium 

apz 

■3?  = <4~6) 


where  is  the  water  density  when  the  water  velocity  is 
zero.  By  substituting  equations  (4-4) , (4-5) , (4-6)  in 
the  term 


1 

P 


dzj 


of  equation  (4-3)  we  have 


(4-7) 


But  from  the  definition  of  the  volumetric  expansion  of 
water  (B)  we  may  write 


B 


1 

i 

p-p  ^ 

oo 

p 

MP 

p 

T-T 

°°J 

P 


(4-8) 


Substituting  equation  (4-8)‘  into  (4-7)  we  have 
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- |J|=  ^ (p  - p ) = ^[-PBCT-T  ) ] = gB  (T-T  ) (4-9) 

P[ZdZjp  oo  p oo  J oo 


Then  the  momentum  equation  (4-3)  can  be  written  as 


3v 


3v 


3t 


+ v 


+ v 


r 3r  z 3z 


= gB (T-T  ) + v 


r32v  , 3v  3 2v 

z . 1 z z 


dr*  + r 3r 


3 z ' 


(4-10) 

An  additional  assumption  can  be  made  regarding  the 

relative  value  of  the  gradient  of  the  velocity  component 

in  the  axial  direction  (v  ) , in  the  r,  and  z directions. 

z 

From  physical  observation,  changes  of  v^  are  much  greater 
in  the  r direction  than  in  the  z direction,  or 


3v  3v 

z z 

>>  , 

3r  3z 


and 


3 2v 
z 

3r  2 


>> 


Therefore,  the  momentum  equation  can  be  further  simpli- 
fied 
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The  development  of  the  energy  equation  is  similar 
to  the  one  derived  on  page  255,  ref.  [23]  except  for  the 
source  term  H which  is  induced  here  by  the  incident 
solar  radiation  on  the  water.  The  energy  equation  is 


pc 


3T 

3t 


3T  , 3T 
+ v - — + v — 
r 3r  z 3z 


k 


3 2T  1ST  3 2t" 

3r2  r 3r  + 3z2j 


+ y$+H. 


(4-12) 

It  is  assumed  that  the  frictional  effects  on  the  energy 
equation  are  negligible,  and  y$=0. 


Description  of  Source  Term 
The  source  term  H is  the  amount  of  solar  radiation 
absorbed  by  the  water  at  different  depths.  The  source 
term  H is  not  only  a function  of  depth  but  also  a 
function  of  the  wavelength  of  solar  radiation.  Stated 
differently,  the  water  absorbs  different  amounts  of 
different  wavelengths  and  at  different  depths.  From  the 
available  incident  solar  radiation  water  absorbs  a fraction 
Y at  the  surface  (z=0) . The  rest  of  the  solar  radiation 
(1-y)  is  absorbed  exponentially  throughout  the  depth 
of  the  water  volume.  The  source  term  H inside  the  water 
volume  is 


H 


(l-Y)  n 


I,  . e 
honz 


-riz 


(4-13) 
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It  will  be  useful  to  examine  the  variation  of  the 
extinction  coefficient,  n»  with  wavelength  of  solar 
radiation  (Table  4-1) . The  data  in  Table  4-1  have  been 
plotted  in  Fig.  4-2.  It  is  evident  that  water  is  very 
transparent  for  radiation  wavelengths  between  0.4y  to 
0.6y,  and  in  the  infrared  the  transparency  is  approximately 
zero.  Therefore,  the  fraction  y represents  all  the  solar 
radiation  other  than  the  visible  which  is  absorbed  by  the 
top  layer.  From  a table  of  black-body  radiation  functions 
[24],  the  visible  portion  (0.4y  to  0.6y)  is  calculated 
to  be  24  percent  of  the  solar  spectrum.  Therefore, 
y = 76  percent. 

Substituting  equation  (4-13)  into  (4-12)  and 
rearranging,  we  have 


3T 

at 


[ar2'  r 3r 


32T  1 3T 

9 * 


(4-14) 


Description  of  Motion  of  Water 


Inside  the  Swimming  Pool 


During  the  daytime  and  nighttime  the  water  inside 


the  pool  moves  due  to  density  differences  which  are  caused 


EXTINCTION  COEFFICIENT  PER  METER  OF  PURE  WATER  AS 
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Wavelength  (y) 


Figure  4-2.  Extinction  coefficient  for  water  as  a func- 
tion of  wavelength;  plotted  from  [25]. 
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The  most  important  heat  transfer  mode  in  these  regions 
is  convection.  The  water  which  is  located  away  from 
the  pool  walls,  the  center  of  the  pool,  and  the  top  and 
bottom  surfaces  of  the  pool  moves  very  slowly,  if  any, 
because  there  are  no  density  gradients  within  these 
regions.  Thermal  energy  is  transferred  only  by  heat 
conduction.  In  these  regions  it  is  assumed  that  the  water 
moves  with  zero  velocity. 

Based  on  this  description  of  the  physical  model  of 
the  water  motion  inside  the  pool,  typical  velocity  profiles 
during  daytime  and  nighttime  are  shown  in  Figure  4-3. 

These  velocity  profiles  suggest  that  the  water  motion 
is  a torroid.  This  is  shown  in  Figure  4-4. 


During  Daytime  During  Night 


Figure  4-4. 


Circulatory  patterns  of  water. 
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Figure  4-3. 


Velocity  profiles  during  day  and  night. 
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It  should  be  noted  that  the  motion  of  the  water 
next  to  the  walls  is  confined  to  a very  small  radial 
distance.  From  experimental  measurements  that  this  author 
conducted  in  a pool  which  was  heated  by  the  sun  the 
maximum  radial  distance  next  to  the  wall  where  there  was 
a noticeable  motion  of  water  was  less  than  1 inch. 

It  is  proposed  then  to  calculate  the  spatial  tempera- 
tures of  the  water  as  a function  of  time  inside  the 
torroid.  This  region  contains  the  major  volume  of  the 
water.  An  assumption  is  next  made  in  the  relative  magni- 
tude of  the  water  temperature  gradient  in  the  radial  and 
axial  direction.  It  is  assumed  that 

9T  _ 9T 
3z  3r 

This  is  a very  good  assumption  during  the  day  when  the  top 
surface  is  always  hotter  than  the  rest  due  to  absorption 
of  solar  radiation.  During  the  night  this  assumption  is 
also  good  because  the  top  surface  is  cooled  a lot  more 
than  the  sides  of  the  pool.  The  top  surface  is  cooled 
primarily  by  evaporation  cooling,  the  sides  are  cooled 
by  convection  and  conduction;  evaporation  cools  the  water 
far  greater  than  convection  or  conduction. 

Since  the  assumption  of  zero  velocity  inside  the 
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torroid  was  already  made,  the  velocity  components  v and 

v in  Equation  (4-14)  are  also  zero.  Incorporating  all 
z 

these  assumptions  into  (4-14)  we  have 


3T  , n /n  \ -nz 

— = a - — y + — (1-y)  I,  . e 
3t  3z  pc  horiz. 


(4-15) 


The  solar  radiation  at  the  top  surface,  I,  . , is  a func- 

^ horiz 

tion  of  time,  and  it  can  be  represented  by  a Fourier 
series  expansion.  Therefore 


I,  • = Cn  + C cos(C  t)  + C sin  (C  _t) 

horiz  1 2 18  3 18 


+ C4  cos  (2Clgt)  + C5  sin  (2Clgt) 


+ Cg  cos  (3Clgt)  + C^sin  (3Clgt) 


+ Cg  cos  (4Clgt)  + Cg  sin  (4Clgt) 


+ C Q cos  (5Clgt)  + C sin  (5Clgt) 


+ C 2 cos  (6Clgt)  + C13  sin  (6Clgt) 


+ C14  cos  (7Clgt)  + Cl5  sin  (7Clgt) 


+ C g cos  (8Clgt)  + C17  sin  (8Clgt) 


(4-16) 
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Substituting 

3 2T  _ 1 3T 
3z  2 a 3t 


(4-16)  into 
-n ( 1~y) e n 


k 


(4-15)  and  rearranging  we  have 
(°1  + •••  + C17Sin(8C18t)j 


(4-17) 

Equation  (4-17)  is  the  Fourier  Heat  Conduction 
Equation  in  the  spatial  variable  (z) ; it  also  contains  a 
source  term.  Equation  (4-17)  can  be  solved  by  using 
Finite  Cosine  Fourier  Transforms  on  the  spatial  variable 
z . Let  the  Finite  Cosine  Fourier  Transform  of  the  water 
temperature  T(z,t)  be  equal  to 


F [T(z,t)  ] E T(w,t) 
\/ 


l 

f T(z,t)  cos 
0 


"wTrzl 

* J 


dz . 


Then  the  general  solution  according  to  ref 


T (z  , t)  = j T ( 0 , 0 ) + j-  £ T (w,  t ) cos 

w=l 


[40]  is 


The  following  boundary  conditions  apply  to  the  water  inside 
the  torroid: 

(1)  Initially  there  is  a spatial  temperature 
distribution  at 

2 

t = 0 -*  T ( z , 0 ) = E z + E z + E . 

•L  Cm 
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(2)  At  the  top  surface  of  the  pool  the  net  energy 
which  is  conducted  to  the  water  is  a function  of  time  and 
can  be  represented  by  a Fourier  Series  expansion  at 

k 3T 

z = ° * TT  = D1  + °2  COS  'V1  + °3  Sin  (CX8t) 

+ D4  cos  (2C18t)  + D5  sin  (2Clgt) 

+ Dg  cos  (3C18t)  + D?sin  (3C  gt) 

+ Dg  cos  (4Clgt)  + D^sin  (4C  t) 

+ D1()  cos  (5Clgt)  4-  Du  sin  (SC^t) 

+ D12  cos  (6Clgt)  + D sin  (6C  t) 

+ D14  cos  (7Clgt)  + D15  sin  (7Clgt). 

(3)  There  is  no  energy  conducted  from  the  water 

to  the  ground . 

ST 

at  z = Z k — = 0 
3z 


After  solving  equation  (4-17)  by  Finite  Cosine 
Fourier  Transforms  in  Appendix  A,  and  applying  the 
three  boundary  conditions,  the  following  equation 
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calculates  the  water  temperature  T(z,t)  inside  the  torroid: 


E1  2 E2  2°°- 
T(z,t)  - — t + — l + + j E T(w,t)  cos 

w=l 


WTTZ 


(4-18) 

where  T(w,t)  is  defined  by  equations  (A-10)  and  (A-13) 
in  Appendix  A. 


Development  of  Equation  for  Bulk 

Water  Temperature  Variation 

The  bulk  temperature  of  the  water  in  the  swimming 

pool  is  raised  or  lowered  according  to  whether  there  is 

a net  gain  or  loss  of  energy  from  or  to  the  water  at  that 

instant.  If  all  energy  transferred  to  or  from  the  water 

can  be  accounted  at  any  interval  of  time  (t  -t  ) , the 

2 1 

bulk  temperature  of  the  water  can  be  predicted  from 


^qgain 


^OSS^ 


mcAT 


(4-19) 


• • 

where  C>GAIN » ^L0SS  = amount  of  energy  gained,  or  lost, 

respectively 

m = mass  of  water  in  swimming  pool 

c = specific  heat  of  water 

AT  = bulk  temperature  change  of  water 


during  time  interval 
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To  account  and  describe  all  the  energy  transferred 
to  or  from  the  water,  it  is  convenient  to  draw  a control 
volume  around  the  total  volume  of  the  swimming  pool  (Figure 
4-5) . The  energy  transferred  to  the  water  as  gain  is  from 
the  solar  collectors  and  from  the  insolation.  The  energy 
transferred  from  the  water  as  loss  during  the  night  is 
the  sky  radiation.  Other  modes  of  energy  that  can  be 
transferred  as  loss  or  gain  are  convection  and  conduction. 
Finally,  evaporation  losses  can  occur  during  night  or  day. 
Once  these  modes  of  energy  have  been  quantitatively 
described,  then  the  time  history  of  the  bulk  temperature 
of  the  water  can  be  calculated. 

The  energy  which  is  transferred  to  the  water  from 
the  solar  collectors  can  be  calculated  from  the  equa- 
tion 

*2 

^collector  " [ Collector  V*’  d*  (4-20> 

*1 


where 


V*> 


’ horiz 


cos (L-s)  cos  6 cos  h + sin  (L-s)  sinfi'l 
cos  L cos  6 cos  h + sin  L sin  6 j 


and  <f>  , <f>  are  initial  and  final  values  of  the  sun's 
incidence  angle  during  the  time  interval  (t^-t  ) . 
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Sky  Radia- 

Connection  Energy  Loss  tion  Gain  from  Gain  from 


Figure  4-5.  Energy  fluxes  crossing  the  control  volume 

of  a swimming  pool. 

i . 
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The  water  in  the  swimming  pool  absorbs  a certain 
amount  of  energy  from  the  incident  solar  radiation.  This 
energy  can  also  be  treated  as  a gain.  The  amount  of  the 
solar  energy  which  causes  a change  in  the  water  tempera- 
ture depends  on  the  absorptivity  of  the  water  and  cover 
(if  any) . The  absorptivity  of  the  water  depends  mainly  on 
the  angle  of  incidence  of  the  solar  rays  on  the  water 
surface.  This  angle  of  incidence  is  a function  of  time. 

The  amount  of  solar  energy  which  causes  a change  in  the 
water  temperature  can  be  expressed  by 

*2 

°sun  “ ' a(*>  Asurf  Jhoriz  <*>  d«  <4-21) 

♦l 

The  energy  which  is  transferred  due  to  convection 
can  be  regarded  as  a gain  or  a loss.  In  a swimming  pool 
both  free  and  forced  convection  occur.  A temperature 
difference  between  the  water  and  the  air  produces  a density 
difference  in  the  water.  This  density  difference  causes 
the  water  to  move.  During  this  motion  the  water 
exchanges  energy  with  the  surrounding  air.  As  a result, 
the  water  temperature  changes.  This  energy  exchange  is 
caused  by  free  convection.  In  forced  convection  the 
driving  force  is  the  wind  which  cools  or  heats  the  top 
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and  side  surfaces  of  the  swimming  pool  (in  case  the 
swimming  pool  is  above  ground) . This  again  causes  the 
water  temperature  to  fall  or  rise.  The  energy  exchange 
due  to  free  and  forced  convection  for  the  vertical  walls 
of  the  swimming  pool  above  ground  can  be  presented  as 


Q 

conv. 

wall 


f2  U A (T 

, wall  water 

1 


T . ) dt 
air 


(4-22) 


Since  the  energy  is  transferred  perpendicular  to  the 
walls  of  the  pool,  the  combined  coefficient  of  heat 
transfer  U can  be  calculated  from 


wall 


(4-23) 


The  magnitude  of  the  convective  heat  transfer  coeffi- 
cient inside  the  swimming  pool  (h_^)  depends  mainly  on 
the  magnitude  of  the  velocity  of  the  water  in  the  case 
when  the  swimming  pool  circulation  pump  is  not  in  use. 

The  water  velocity  is  caused  by  a density  difference  in 
the  water;  this  density  difference  depends  directly  on 
the  temperature  difference  between  the  water  and  the  wall 
of  the  swimming  pool.  Therefore,  the  magnitude  of  h. 

l 

depends  on  this  temperature  difference  between  the  water 
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and  the  wall  of  the  swimming  pool.  The  equation  express- 
ing this  relationship  is  derived  as  follows.  Assume 
the  temperature  of  the  vertical  wall  of  the  swimming 
pool  to  be  T^;  the  temperature  of  the  bulk  of  the  water 
Tco'  a^so  T > T^.  The  average  water  temperature  near 
the  vertical  wall  is  (T  + TJ/2.  The  temperature 
difference  between  the  water  near  the  wall  and  the  bulk 
of  the  water  is 


T + T T,  - T 

1 00  _ 1 00  AT 

2 00  2 2 * 

Since  the  water  near  the  wall  is  warmer  than  the  bulk 
of  the  water,  the  water  near  the  wall  is  lighter  and 
rises  due  to  the  buoyancy  force  F, 


F 


(4-24) 


From  equation  (4-8)  we  may  write  that 


F 


pB 


f AT  _g_ 

l 2 J 9C  ' 


(4-25) 


When  this  buoyancy  force  acts  through  a height  l,  then 
work  is  done.  The  work  is- 
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z=Z 

W = f F 
z=0 


dz  = ~ P$  l. 

q 2 


Assuming  that  all  the  work  is  converted  into  kinetic 
energy,  then 


This  assumption  is  discussed  in  Chapter  IX,  Discussion 
and  Evaluation  of  Results.  Solving  for  V we  have 


V = /g3AT£  (4-26) 

The  well-known  Nusselt  Equation  relates  the  water 
velocity  with  the  inside  convective  heat  transfer  coeffi- 
cient as  follows:  Nusselt 's  Equation 


h.l 

l 

k 


(4-27) 


Solving  for  h,  and  rearranging  Equation  (4-27)  we 

have 


h. 

i 


Ck 

\p2  2 2 

l V p 

b/2 

r > 

c y 
p 

l 

2 

y J 

-<  J 

(4-28) 


Substituting  Equation  (4-26)  into  Equation  (4-28)  we 


have 
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h. 

1 


Ck 

l 


3 2 

gBAT l p 

b/2 

c p'' 
P 

2 

k 

W J 

d 


or 


h. 

1 


Ck 

l 


(Gr) 


b/2 


(Pr)  d 


(4-29) 


Brown  and  Marco  [26]  report  that  the  coefficients  b/2  and 
d are  numerically  the  same  from  free  convection  tests. 
Then  equation  (4-29)  may  be  rewritten  as 

f 

(4-30) 


, Ck  . f k 

hi  = T (Gr'Pr)  = CI 


( 3 2 1 

ggATf  p c 

E 


yk 


From  tests  Saunders  [27]  estimates  the  exponent  f to  be 

3 9 

f = 0.25  in  the  range  of  (Gr  x Pr)  =10  -*  10  . The 

results  of  many  investigators  were  used  to  calculate 
the  empirical  constant  C of  equation  (4-30)  . From  ref. 
[26]  the  value  of  C may  vary  as  following: 

C = 0.55  for  vertical  walls  or  plates 
C = 0.71  for  horizontal  plates,  warm  side  up 
C = 0.35  for  horizontal  plates,  warm  side  down. 

It  has  been  found  experimentally  that  a length  of  a 
surface  which  is  greater  than  2 feet  has  a negligible 
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effect  on  h^.  Therefore,  the  length  1 in  equation  (4-30) 
should  not  exceed  2 feet  [26]. 

From  equation  (4-30)  the  value  of  h^  can  be  calcu- 
lated and  substituted  into  equation  (4-23).  However,  the 

convective  heat  transfer  h needs  to  be  calculated  in 

o 

order  for  the  value  of  U in  equation  (4-23)  to  be  calcu- 
lated. Empirical  data  are  reported  in  ref.  [28]  for  the 
convective  heat  transfer  coefficient  for  stucco  and  glass 
versus  wind  speed.  These  data  have  been  replotted  in 
Figure  4-6  after  removing  the  radiant  portion  of  the 

convective  heat  transfer  coefficient.  The  radiant  portion 

2 

was  about  0.7  BTU/(hr  ft  F)  for  the  conditions  of  the 
tests  as  reported  by  ref.  [29].  Once  the  value  of  hQ  is 
found,  the  combined  coefficient  of  heat  transfer  (U)  in 
equation  (4-23)  can  be  calculated.  It  should  be  noted  that 
the  value  of  U is  time  dependent,  because  the  values  of 
h^  and  hQ  are  time  dependent. 

Equation  (4-22)  calculates  the  energy  exchange  due  to 
free  and  forced  convection  for  the  vertical  sides  of  a 
swimming  pool  above  the  ground.  There  is  a similar 
equation  for  the  energy  exchange  due  to  free  and  forced 
convection  for  the  top  surface  of  a swimming  pool.  This 
is  given  by 


■V 


Convective  Heat  Transfer  Coefficient,  h , (BTU/ (hr-sq . f t . F) ) 


14 


Wind  Speed  (Miles/Hour) 

Figure  4-6.  Convective  heat  transfer  coefficient  vs. 
wind  speed. 
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Q. 


conv 

top 


/ U A (T  - T . ) dt  (4-31) 

1 surf  water  air 


where 


(4-32) 


The  value  of  U,  can  be  calculated  once  the  values  of  h. 

1 i 

and  h are  found.  It  should  be  noted  that  h.  now  is 
o 1 

the  value  of  the  convective  coefficient  near  the  top 
surface  of  the  swimming  pool.  The  value  of  h^  is  calculated 
from  Equation  4-30  after  selecting  the  proper  value  for 
the  constant  C.  The  value  of  hQ  can  be  found  again  from 
Figure  4-6. 

Equations  4-22  and  4-31  calculate  the  energy  exchange 
due  to  convection  from  the  vertical  walls  and  the  top 
surface  of  the  swimming  pool.  Another  energy  exchange 
is  by  conduction  from  the  bottom  side  of  the  swimming  pool 
to  the  ground.  This  energy  exchange  can  be  calculated 
from  Fourier's  equation  for  conduction 


Q 


cond 


/ 

tn 


k A 


dT 
surf  dz 


dt 


(4-33) 


bottom 


The  temperature  gradient  in  the  z-direction  refers  to 
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the  rate  of  change  of  temperature  T of  the  water  tempera- 
ture at  the  bottom  surface  of  the  pool. 

Another  energy  exchange  is  due  to  the  water  emitting 
energy  in  the  form  of  radiation  towards  the  sky.  The 
sky  temperature  is  usually  lower  than  the  temperature  of 
t the  water,  especially  during  the  night.  There  is  not  a 
general  agreement  among  scientists  in  determining  the 
value  of  the  sky  temperature.  Fischenden  and  Saunders 
[30]  reported  the  value  of  -50°F  for  the  temperature  of 
the  sky  during  cold  clear  nights.  Czarnecki  [8]  estimated 
the  mean  sky  temperature  to  equal  the  mean  air  temperature 
less  20°F.  Root  [31]  reported  the  sky  temperature  to 
vary  from  -100°F  for  clear  cold  nights  to  the  value  of 
the  air  temperature  for  cloudy  humid  nights.  Root's 
estimates  were  based  on  experimentation  by  the  Florida 
Citrus  Industry  in  radiation  losses  from  orange  groves  in 
1952  [31].  The  effective  sky  temperature  is  higher 

during  cloudy  humid  nights.  The  water  in  the  swimming 
pool  loses  less  energy  due  to  sky  radiation  losses  during 
cloudy  humid  nights  because  the  water  vapors  in  the 
clouds  radiate  back  to  the  ground.  Due  to  the  geo- 
graphical proximity  of  Root's  research,  his  results  will 
be  adopted  in  calculating  the  sky  temperature  in  this 
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investigation.  The  equation  which  gives  the  sky  radiation 
losses  is  derived  from  the  following  considerations. 

The  amount  of  radiation  emitted  by  the  top  surface  of  a 
swimming  pool  is 


Q . = f 

emi  t 

t. 


A ea  T dt 

surf  top 


(4-34) 


The  amount  of  radiation  absorbed  by  the  top  surface  from 
the  sky  is 


Q 


absor 


4 

A acr  T 
surf  sky 


dt 


(4-35) 


Then  the  net  amount  of  radiation  crossing  the  control  volume 
(Figure  4-5)  during  the  time  interval  ( t 2 - t^)  is 
obtained  by  combining  Equations  (4-34)  and  (4-35)  as 
follows  : 


^sky  ^emit  ^absor^ 

rad 


t t 

2 4 2 4 

/ A sax  dt  - / A a a T4  , dt 

surf  top  surf  sky 

1 fcl 


= - A a / (eT4  - a T4  ) dt 
surf  ■ top  sky 


(4-36) 
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Assuming  that  the  water  can  be  regarded  as  a gray  body, 
then  a = e.  Equation  (4-36)  becomes 

t 


sky 

rad 


A a 
surf 


I 

t. 


e (T 


top 


dt 


(4-37) 


The  energy  exchange  due  to  mass  transfer  out  of  the 
control  volume  (Figure  4-5)  pertains  to  the  evaporation 
of  water  from  the  top  surface  of  the  swimming  pool.  Water 
evaporates  whenever  the  value  of  the  relative  humidity  in 
the  surrounding  air  is  less  than  100  percent.  The  water, 
in  order  to  evaporate , removes  heat  from  the  surrounding 
water.  This  removal  of  heat  from  the  top  surface  of  the 
water  always  occurs  whenever  evaporation  is  permitted. 

Many  researchers  have  investigated  the  rate  of  evaporation 
of  water.  The  appropriate  mass  transfer  equation  is 


G = C (p  - p ) 
w a 

or 


G = c'Ap  (4-38) 

Where  many  of  these  researchers  disagree  is  on  the 

value  of  the  experimental  constant  c'.  The  reason  of 

the  disagreement  may  be  due  to  different  experimental  set- 
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ups,  different  temperature  ranges  of  the  water  and  differ- 
ent air  velocity  ranges.  The  variation  of  c'  is  given  as 
follows : 


Reference 

[30]  : 

c ' 

= 22.9 

+ 

17.4 

u ’ 

Reference 

[30]  : 

c ' 

= 22.0 

+ 

16.8 

u 

Reference 

[31]  : 

c ' 

= 17.8 

+ 

15.2 

u 

to  be  used  in 

Equation  (4-38) 

Reference 

[32]  : 

c ' 

= 31.0 

+ 

13.5 

u 

with  metric  units 

Reference 

[33]  : 

c 1 

= 18.3 

+ 

11.0 

u 

0 

J 

. 8 

Reference 

[29]  : 

c ' 

= 1.08 

(up.  ) 

to  be  used  in 

Equation  (4-38) 

Reference 

[8]  : 

c ' 

= 0.20 

h 

n 

with  English 

j 

units 

The  cooling  effect  due  to  vaporization  can  be  calculated 
by  multiplying  the  evaporation  rate  (G)  times  the 
enthalpy  of  vaporization  of  water.  The  equation  is 


Q 


evap 


- S 

t., 


G H A _ dt 
fg  surf 


(4-39) 


Since  all  energy  terms  in  Figure  4-5  have  been 
discussed  and  equations  have  been  presented  to  evaluate 
them,  it  is  now  appropriate  to  write  the  energy  balance 
for  the  water  in  a swimming  pool  above  ground.  Equation 
(4-19)  gives  the  general  form  of  the  energy  balance 
in  a time  interval  (t2  - t^)  . 
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(QGAIN  ~ qloss) 


mcAT 

(t2"tl) 


(4-19) 


Substituting  the  above  derived  equations  into  (4-19)  we 
have 

( +Q  , + Q + Q + Q + Q , - Q 

collector  sun  conv  conv  cond  sky 

wall  top  rad 


- Q_  ) = 


mcAT 


evap  (t2_tl) 


(4-40) 


or 
*2 

+/  <M4>)  A 


4>. 


, , . I.  (4>)  d(j)  + / a (<p)  A _ i.  . (4))  d<|)+ 

collector  t r x T surf  horiz  Y ^ 

*1 


f u A (T 

wall  water 

1 


T . )dt  +/  U A _ (T  ^ -T  . ) dt  + 
air  , 1 surf  water  air 

tl 


2 

/ k A 
t„ 


dT 
surf  dz 


bottom 


t„ 

^ 4 

dt  - A c if  £ (T 

surf  top 

1 


T , ) dt 

sky 


2 

■/  G H 
t. 


.A  dt 

fg  surf 


mcAT 

<vv 


(4-41) 


It  should  be  noted  that  during  the  night 
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O =0 

^collector 


Q = 0. 

sun 


Equation  (4-41)  can  be  readily  solved  for  the  only  unknown 


variable  T. 


water 


CHAPTER  V 


THEORETICAL  RESULTS 

For  theoretical  results  to  be  meaningful,  they  have 
to  be  compared  to  data  which  were  directly  measured  by 
experiments,  or  data  which  were  directly  derived  from 
experimental  measurements.  The  climatic  data  which  are 
required  in  calculating  the  spatial  or  the  bulk  tempera- 
ture variations  will  be  provided  from  measured  climatic 
data  at  the  test  site.  In  this  Chapter,  the  spatial 
and  bulk  temperature  variations  are  calculated  for  a 
given  set  of  climatic  data. 

Spatial  Water  Temperature  Variations 
The  climatic  data  and  initial  boundary  conditions, 
which  are  to  be  provided  as  input  in  Equation  (4-18)  , 
are  (a)  the  temporal  variations  of  the  solar  radiation 
incident  on  a horizontal  surface,  (b)  the  temporal  varia- 
tions of  the  net  energy  transferred  to  the  water  at  the 
top  surface  of  the  pool,  and  (c)  the  initial  spatial 
water  temperature  distribution.  These  data  were  provided 
from  measurements  taken  at  the  University  of  Florida, 
Solar  House  on  12/1/72;  measurements  were  taken  each 
hour  from  0600  hours  to  2100  hours  of  the  same  day. 
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The  temporal  variations  of  the  solar  radiation  and  .of 
the  net  energy  transferred  to  the  water  are  shown  in 
Figures  5-1  and  5-2,  respectively.  The  initial  spatial 
water  temperature  distribution  at  0600  hours  is  shown  in 
Figure  5-3. 

A Fortran  program  was  written  for  equation  (4-18) . 

This  Fortran  program  appears  in  Appendix  B.  The  converg- 
ence of  the  solution  was  quite  quick.  Only  the  first 
200  terms  in  the  series  of  equation  (4-18)  were  considered. 
The  error  which  is  introduced  by  considering  the  first 
200  terms  is  0.0033%.  Water  temperatures  are  calculated 
as  a function  of  depth  for  each  hour  during  a 24  hour 
interval.  These  temperatures  are  plotted  in  Figures 
5-4  and  5-5. 

It  should  be  noted  that  the  water  temperatures  in 
the  upper  6 inches  of  a 40-inch  high  swimming  pool  were 
not  plotted  because  they  were  not  calculated.  Equation 
(4-18)  calculates  the  spatial  water  temperatures  inside 
the  torroid.  The  cooling  and  heating  periods  of  the 
water  can  be  clearly  seen  in  Figures  5-4  and  5-5.  On 
12/1/72,  the  cooling  ended  around  1000  hours,  and  it 
started  again  after  1700  hours  local  time.  This  is  as 
expected  because  of  the  short  duration  of  winter  days. 
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Figure  5-2.  Net  energy  transferred  to  water  at  the  surface 
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Water  Temperature  (F) 


Figure  5-3. 


Initial  temperature  distribution  at 
0600  hours  on  12/1/72. 
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Figure  5-4.  Spatial  Water  Temperature  Profile  from  0600  hours-  to 

1800  hours  on  12/1/72. 
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Figure  5-5.  Spatial  water  temperature  profiles  from 

1800  hours  on  12/1/72  to  0500  hours  on 
12/2/72. 
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Equation  (4-26)  calculates  the  maximum  velocity  of 
the  water  which  is  located  next  to  a heated  or  cooled  pool 
wall.  To  gain  a quantitative  measure  of  this  velocity, 
the  climatic  data  and  water  temperatures  on  12/1/72  were 
used  in  Equation  (4-26) . The  results  were  plotted  in 
Figure  5-6. 


Bulk  Water  Temperature  Variations 
Equation  (4-41)  calculates  the  temporal  variations 
of  the  bulk  water  temperature  in  an  unheated  swimming 
pool  and  in  a pool  heated  by  solar  collectors . These 
temporal  variations  can  be  examined  in  three  time  frames, 
annual,  monthly,  and  daily.  In  an  annual  time  frame 
the  bulk  water  temperature  for  each  month  can  be  predicted 
if  the  monthly  average  information  in  climatic  conditions 
is  provided  as  input  to  Equation  (4-41) . In  a monthly 
time  frame  the  bulk  water  temperature  can  be  predicted 
during  several  days  for  each  month.  In  a daily  time 
frame  the  bulk  water  temperature  can  be  predicted  every 
hour . 

Annual  Variations 

To  calculate  the  variation  of  the  bulk  water 


temperature  on  an  annual  time  frame , two  actions  need  to 
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Figure  5-6.  Calculated  maximum  water  velocity  next  to  the  pool  walls  vs.  time 
for  uncovered  pool  on  12/1/72. 
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be  taken.  One  is  to  transform  Equation  (4-41)  into  a more 
useful  and  convenient  form.  The  second  is  to  gather 
adequate  climatic  data  to  provide  them  as  input  to  Equa- 
tion (4-41) . 

Equation  (4-14)  is  the  general  heat  balance  equation 
which  could  be  used  for  either  the  unheated  or  the  solarly 
heated  swimming  pool  during  the  day  or  night.  This 
equation  can  be  simplified  by  considering  the  conduction 
losses  to  the  ground  to  be  negligible.  This  is  a reason- 
able assumption  as  it  will  be  shown  later  in  Chapter 
VIII  under  Experimental  Results.  A second  simplification 
is  achieved  by  transforming  the  quantities  under  the 
integral  signs  to  time-averaged  quantities.  A third 
simplification  is  achieved  by  assuming  that  the  average 
monthly  ambient  air  temperature  is  lower  than  the  average 
monthly  bulk  water  temperature.  This  is  also  true  as 
it  will  be  shown  later  in  Chapter  VIII  under  Experimental 
Data.  Since  the  time  frame  is  the  annual,  it  can  be 
stated  that  the  average  equilibrium  temperature  of  the 
bulk  of  the  water  is  constant.  This  implies  that 
Equation  (4-19)  can  be  modified  as 


'GAIN 


- Q 


LOSS 


= 0 


(5-1) 
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Utilizing  the  three  simplifications  as  discussed  above, 
along  with  Equation  (5-1)  , Equation  (4-41)  can  be 
transformed  to  the  following: 


tl  Acollector  It  + a Asurf  "^horiz  + U Awall  ^Twater  Tair^ 


+ UlAsurf(Twater-Tair)1  ' t2a  Asurfe  ('I'4top-T4sky) 


- t_  G H.  A _ = 0 
3 fg  surf 


(5-2) 


where 


t 


1 


t 


2 


t 


3 


= average  number  of  daylight  hours  for  a certain 
month 

= average  number  of  night  time  hours  for  a 
certain  month 
= 24  hours 


All  bar  quantities  in  (5-2)  are  average  monthly  values 

For  Equation  (5-2)  to  be  solved  for  the  average  monthly 

bulk  water  temperature  (Twater) , an  initial  assumption 

needs  to  be  made  on  the  value  of  T . Once  T 

top  water 

has  been  calculated,  a simple  heat  balance  equation  can 

calculate  r7  , and  the  calculation  of  T is 

top  water 

repeated  as  necessary.  For  Equation  (5-2)  to  be 

solved , input  data  are  also  needed.  The  average  monthly 
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air  temperature  is  obtained  from  Table  5-1.  The  average 
daily  air  temperatures  were  calculated  from  the  record 
of  a Belfort  hygrometer  which  was  located  at  the  test  site. 
The  average  monthly  and  daily  insolation  data  are  obtained 
from  Table  5-2.  These  values  of  insolation  are  the  sum 
of  the  direct  and  diffuse  insolation  falling  on  a hori- 
zontal surface . Values  for  the  apparent  absorptivity 
are  obtained  from  data  under  Experimental  Results  in 
Chapter  VIII.  The  combined  heat  transfer  coefficient 
for  the  side  walls  (U)  can  be  calculated  from  (4-23) 

. T 1 

t] 

J plastic 


*xT 

kJ 


X o 


steel 


where 


^x'l 


steel 


32x12 


26 


0.00010 


fx'l 


^ Aplastic 


32x12 

0.05 


= 0.05208 


Then  Equation  (4-23)  becomes 


U = 


~ + — + 0.0522 
h . h 
1 o 
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TABLE  5-1 

AVERAGE  AIR  TEMPERATURE  (F) 
MEASURED  AT  THE  TEST  SITE 


Date 

April 

’72 

May 

•72 

June 

'12 

July 

'12 

Aug. 

'12 

Sept. 
1 72 

Oct. 

'12 

Nov. 

'72 

Dec. 

'12 

' Jan . 

' 73 

Feb. 

’73 

March 

’73 

1 

58.0 

. 71.5 

7 3.5 

82.5 

78.5 

' 77.0 

73.5 

72.5 

50.0 

67.5 

64  .5 

58.0 

2 

60.0 

71.0 

73.5 

82.5 

78.5 

75.5 

75.0 

73.0 

52.0 

68.0 

66.5 

65.5 

3 

62.0 

72.5 

73.0 

83.0 

79.0 

— 

75.5 

74.5 

59  .0 

69.5 

53.5 

67.0 

4 

68.0 

72.5 

76.5 

80.5 

82.5 

— 

77.0 

74.5 

63.5 

72.0 

48.5 

6 7.5 

5 

71.0 

71 .0 

75.5 

82.0 

83.5 

83.0* 

77.5 

74.5 

69.0 

71.0 

50.5 

70  .0 

6 

69.0 

73.5 

77.5 

76.5 

81.0 

84.0* 

74.5 

73.0 

73.0 

68.0 

56.0 

71.5 

7 

69.5 

73.0 

79.0 

77.0 

81.0 

81.0* 

70.0 

71.0 

69.0 

58.5 

60.0 

72.5 

8 

65.5 

75.0 

77.0 

78.0 

81.0 

81.0* 

69.0 

65.5 

70.5 

58.0 

62.5  . 

67.5 

9 

58.0 

75.0 

79.0 

73.0 

80.5 

83. C* 

70.0 

61.0 

— 

39.0 

56.0 

68.5 

10 

63.5 

73.0 

80.0 

76.0 

80.5 

82.0* 

73.0 

66.5 

70.5 

41.0 

34.5 

69,0 

11 

66. C 

73.5 

77.0 

75.0 

80.5 

7S.0* 

73.5 

73.5 

69.0 

— 

40.5 

73.0 

12 

69.0 

73.0 

75.5 

30.5 

78.5 

82.0* 

73.5 

72.5 

70.5 

40.5 

44.5 

73.0 

13 

70.5 

72.5 

79.0 

82.0 

79.0 

82.0* 

72.0 

75.0 

73.0 

41.0 

55.0 

72.0 

14 

78.0 

76  ,0 

79.5 

79.0 

79.0 

82.0 

— 

70.5 

71.5 

42.5 

65.0 

75  ,0 

15 

77.0 

74.5 

78.5 

78.5 

81.5 

83.5 

— 

57.0 

73.5 

48.5 

60.5 

75.5 

16 

75  ,.5 

76.5 

79.0 

78.0 

81.5 

83.5 

— 

56.0 

48.5 

54.5 

46.0 

74.5 

17 

75.5 

75.0 

80.0 

78.5 

81.5 

83.5 

75.0 

49.0 

37.5 

59.0 

40.5 

59  .0 

18 

75.0 

75.5 

74.5 

80.5 

83.5 

82.5 

75.0 

57.0 

57.0* 

66.0 

43.5 

55.0 

19 

75.5 

72.0 

77.5 

76  .5 

82.0 

80.5 

76.0 

64.5 

63.0* 

64.5 

49.0 

61.0 

20 

74  .5 

70.5 

76.0 

80.5 

81.0 

79.0 

75.0 

61.0 

65.0* 

55.0 

51.5 

61.0 

21 

73.0 

71.5 

83.5 

81.5 

78.5 

78.0 

61.5 

51.5 

61.0* 

61.5 

51.0 

66.5 

22 

73.5 

73.5 

31.0 

81.5 

78.5 

73.5 

62.5 

51.5 

65.0* 

66  .5 

51.0 

56.5 

23 

70.5 

76.0 

79.0 

83.0 

75.0 

75.5 

71.0 

49.0 

57.0* 

56.0 

53.0 

59.5 

24 

72.5 

75.0 

76.0 

82.5 

72.5 

78.5 

73.0 

51.5 

— 

58.5 

55.0 

59.0 

‘Estimated  from  maximum  and  minimum  temperature  readings  taken  at  the 
Gainesville  airport 


Table  5-1  continued 
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Date 

April 

•72 

May 

'll 

June 

'll 

July 

'll 

Aug  . 
' 72 

Sept . 
'll 

Oct . 
'll 

Nov. 

'll 

Dec . 
'll 

Jan . 
'73 

Feb. 

'73 

March 

'73 

25 

64.5 

72.0 

75.5 

75.5 

72.0 

79.0 

74.5 

61.0 

50.0* 

■53.5 

60.0 

62.5 

26 

63.0 

71.0 

79.0 

76  .5 

72.0  . 

80.5 

72.5 

54.5 

50.0* 

63.0 

64.5 

6 3.5 

27 

63.0 

63.0 

78.0 

83.0 

77.0 

81.5 

69.0 

— 

50.0* 

62.0 

62.5 

60.5 

28 

66  .5 

70.0 

81.0 

84.5 

79  .0 

80.5 

75.0 

61.0 

58.0* 

60.0 

57.0 

63.0 

29 

65.5 

12.5 

82.0 

83.0 

80.5 

80.5 

73.5 

64.0 

62.0* 

42.0 

64.5 

30 

68.5 

75.0 

82.0 

81.0 

77.5 

77 .0 

65.0 

59.0 

68.0* 

41.5 

68.0 

31 

75.0 

75.5 

74.5 

70.0 

67.0 

50. a 

73.0 

MEAN 

68.85 

72.95 

77.93 

82.85 

79.06 

80.27 

72.23 

63.62 

61.83 

56.62 

53.66 

66.21 

‘Estimated  from  maximum  and  minimum  temperature  readings  taken 
Gainesville  airport. 


at  the 


TABLE  5-2 
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HEMISPHERIC  SOLAR  RADIATION 
ON  A HORIZONTAL  SURFACE  (.LANGLEYS) 


Station:  University  of  Florida 

Gainesville , Florida 


Year  - 1972 


Year  - 1975 


Month 


Month 


4 

5 

6 

7 

8 

9 

10  . 

11 

12 

1 

2 

3 

Date 

1 

374 

430 

743 

565 

419 

410 

293 

321 

286 

238 

115 

39  7 

2 

660 

461 

687 

627 

476 

573 

329 

314 

224 

272 

45 

413 

3 

660 

448 

704 

667 

541 

513 

288 

262 

39 

227 

405 

458 

4 

668 

594 

677 

613 

525 

527 

346 

236 

210 

215 

395 

443 

5 

598 

627 

722 

543 

590 

470 

344 

34  2 

159 

293 

425 

367 

6 

537 

417 

744 

346 

336 

313 

451 

17 

169 

— 

400 

308 

7 

555 

256 

744 

579 

4 30 

438 

456 

73 

104 

— 

349 

290 

8 

114 

255 

589 

631 

344 

475 

469 

39  8 

186 

39 

36  3 

156 

9 

341 

411 

563 

429 

372 

452 

461 

38S 

299 

66 

42 

147 

10 

607 

49  8 

480 

488 

427 

405 

411 

222 

304 

45 

46 

282 

11 

626 

486 

458 

57  3 

518 

4 80 

327 

287 

288 

105 

438 

318 

12 

671 

268 

464 

555 

533 

400 

417 

244 

284 

97 

358 

447 

13 

668 

2 78 

6 35 

457 

610 

504 

446 

236 

277 

333 

334 

330 

14 

652 

727 

615 

339 

592 

488 

415 

149 

208 

358 

270 

488 

15 

637 

238 

65C 

515 

518 

522 

408 

335 

162 

332 

400 

410 

16 

484 

606 

535 

359 

548 

507 

334 

299 

231 

372 

438 

453 

17 

722 

624 

694 

460 

517 

489 

367 

173 

296 

310 

393 

5 38 

18 

632 

648 

92 

607 

591 

457 

397 

297 

242 

316 

135 

577 

19 

502 

735 

— 

328 

462 

354 

335 

19 

231 

293 

178 

567 

20 

660 

597 

171 

644 

544 

431 

249 

274 

220 

356 

374 

215 

21 

685 

7 45 

568 

580 

480 

44  6 

337 

313 

46 

202 

308 

556 

22 

544 

717 

443 

593 

506 

415 

329 

180 

106 

350 

480 

572 

23 

153 

660 

— 

5 85 

314 

428 

43  2 

176 

104 

168 

333 

586 

24 

402 
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The  combined  heat  transfer  coefficient  for  the  top 
surface  of  the  swimming  pool  from  Equation  (4-32)  is 
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Two  types  of  plastic  covers  were  primarily  used.  Both 
are  made  of  polyethylene.  The  first  one  will  be  referred 
to  as  "Sears  Plastic"  because  it  is  readily  available 
at  a Sears  Department  Store,  the  second  one  as  "Bubble 
Plastic"  because  it  has  3/4  inch  air  bubbles  (Figure  5-7) . 
The  Bubble  plastic  has  twice  the  thickness  of  the  Sears 
plastic.  Four  arrangements  of  these  plastic  covers  were 
used.  The  value  of  for  each  cover  arrangement  is  as 
follows : 

(a)  One  layer  of  Sears  plastic 
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rr-  + r-  + 0.05208 
h . h 
1 o 


(b)  Two  layers  of  Sears  plastic  with  a 3/4"  air 

space  between  the  two  layers.  According  to  ref. 

[29 J the  resistance  to  heat  flow  for  a 3/4"  air  space 

2 

is  R=0 . 76  (hr-ft  -F)/BTU.  Therefore, 
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Figure  5-7. 


"Bubble"  plastic. 
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jp  + jp  + (2)  (0.05208)  + 0.76 
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r~  + ^-  + 0.864 
h . h 


i o 


(c)  One  layer  of  Bubble  plastic;  U1 , will  be 
approximately  the  same  as  for  the  case  of  two  layers 
of  Sears  plastic,  and  therefore. 


It  is  true  that  the  smaller  volumes  of  the  entrapped 
air  in  the  Bubble  plastic  should  provide  greater  resis- 
tance to  heat  transferred  through  the  plastic.  However, 
this  benefit  is  outweighed  by  the  absence  of  any  entrapped 
air  in  the  region  around  the  air  bubbles  (Figure  5-7) . 
Therefore,  the  net  effect  on  the  resistance  is  approxi- 
mated to  that  of  two  layers  of  Sears  plastic. 

(d)  Two  layers  of  Bubble  plastic.  The 
conservative  approach  is  to  calculate  when  the 
two  layers  of  Bubble  plastic  provide  minimum  resis- 
tance to  heat  transfer.  This  happens  when  the  air  bubbles 
of  the  top  layer  touch  the  spacing  between  the  air 
bubbles  of  the  bottom  layer.  Maximum  benefits  would 
be  realized  when  the  air  bubbles  of  the  top  layer  touch 
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the  air  bubbles  of  the  bottom  layer.  Under  the  minimum 
benefits  condition  will  be  approximately  the  same  as 
for  four  layers  of  Sears  plastic  with  a 3/4"  air  space. 
Therefore , 
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Computer  Program 

To  facilitate  the  computation  of  the  average  monthly 
temperature,  a computer  program  was  written  in  Fortran 
language.  The  listing  of  the  program  appears  in  Appendix 
C.  The  following  flow  chart  was  followed  for  this 
program. 

The  average  monthly  water  temperature  was  computed 
for  the  following  six  different  cases: 

1.  A swimming  pool  which  has  no  plastic  cover  and 
is  not  heated  by  solar  collectors. 

2.  A swimming  pool  which  is  covered  with  a single 
layer  of  Sears  plastic  and  is  not  heated  by 
solar  collectors. 

3.  A swimming  pool  which  is  covered  with  two  layers 
of  Sears  plastic  and  is  not  heated  by  solar 


collectors . 
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FLOW  CHART 


Start 


| Display  calculated  quantities 


End 


4*  A swimming  pool  which  is  covered  with  a single 
layer  of  Bubble  plastic  and  is  not  heated  by 
solar  collectors. 

5.  A swimming  pool  which  is  covered  with  two  layers 
of  Bubble  plastic  and  is  not  heated  by  solar 
collectors . 

6.  A swimming  pool  which  is  covered  with  a single 
layer  of  Sears  plastic  and  is  heated  by  ten 
solar  collectors  (the  area  of  the  ten  collectors 
is  equal  to  270  sq . f t . ) . 

Initially,  the  monthly  water  temperature  was  assigned 
an  assumed  value.  This  was  necessary  to  help  calculate 
the  water  properties  which  are  to  be  used  in  the 
calculation  of  Ik  , tL^  and  U.  After  the  monthly  water 
temperature  is  finally  calculated,  then  it  is  checked 
with  the  assumed  water  temperature  value.  If  there  is 
not  an  agreement  within  0.1°F,  then  the  mean  value  of 
the  calculated  and  of  the  assumed  water  temperatures  is 
used  in  the  calculation  of  the  water  properties.  It 
was  found  that  there  was  a very  quick  convergence  (less 
than  four  iterations) . 

In  the  computer  program  (Appendix  C) , there  are 
nine  property  tables  for  water  and  air  which  cover  a 
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temperature  range  of  40°F  to  130°F  in  increments  of  10°F. 
Within  this  10°F  increment,  the  water  and  air  properties 
were  expressed  as  second  degree  polynomials  of  the  fluid 
temperature . 

The  calculation  of  h^  is  accomplished  by  employing 
Equation  (4-30)  with  the  constant  C = 0.55  for  vertical 
walls  and  f = 0.25  in  Equation  (4-30) . Once  h.  values 
have  been  calculated  for  all  six  cases , the  combined 
heat  transfer  coefficients  are  calculated.  Then  the 
incoming  energy  flux  from  the  sun  and  the  outgoing  flux 
to  the  sky  can  be  calculated.  Finally,  application  of 
Equation  (5-2)  calculates  the  monthly  water  temperature. 

Results 

The  annual  variation  for  the  monthly  water  tempera- 
ture is  shown  in  Figure  5-8.  All  curves  exhibit  a maximum 
water  temperature  around  the  month  of  July  and  a minimum 
water  temperature  around  January.  The  highest  water 
temperature  is  achieved  by  the  swimming  pool  with  single 
Sears  cover  and  being  heated  by  the  ten  collectors.  The 
minimum  water  temperature  is  achieved  by  the  uncovered 
pool.  The  water  temperature  of  the  pool  with  the  single 
Sears  cover  and  of  the  pool  with  the  double  Bubble  cover 
vary  almost  identically.  The  water  temperature  of  the 
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Figure  5-8. 


Calculated  monthly  bulk  temperatures. 
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pool  with  the  double  Sears  cover  and  of  the  pool  with  the 
single  Bubble  cover  vary  identically.  Generally,  the 
water  of  the  pool  with  the  double  Sears  cover  is  at  a 
few  degrees  higher  temperature  level  than  of  the  pool  with 
the  single  Sears  cover. 

To  understand  the  reasons  of  a bulk  temperature 
increase  or  decrease,  one  has  to  consider  all  the  energy 
transferred  across  the  control  volume  of  the  swimming  pool. 
This  is  accomplished  in  Figures  5-9,  10,  11  where  the 
solar  energy  absorbed  by  the  water,  the  convection,  and 
sky  radiation  losses  from  the  water  are  all  shown 
separately . 

In  Figure  5-9,  the  calculated  solar  energy  which  is 
absorbed  by  the  water  is  shown  as  a function  of  time  for 
various  swimming  pools . The  greatest  amount  of  solar 
energy  is  absorbed  by  the  pool  which  is  heated  by  the  ten 
flat  plate  collectors  and  which  is  also  covered  with  a 
single  Sears  cover.  The  reason  is  that  it  collects 
more  solar  energy  more  efficiently  and  also  occupies  greater 
area  than  a 15  foot  diameter  uncovered  pool.  Next,  the 
uncovered  pool  absorbs  the  greatest  amount  of  solar 
radiation  compared  to  any  of  the  covered.  The  covered 
pools  absorb  according  to  their  apparent:  transmissivity 


Vi 

(3 

<U 

X 


5 

c 


•H 


s- 


(Aea/nxs) 


OT  x js'-jem  Aq  paqaosqv  A6aaua  usios 


Figure  5-9 . Solar  energy  absorbed  by  water  as  a function  of  time 
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Figure  5-10  . Convection  losses  vs . time . 
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Figure  5-11. 


Sky  radiation  losses  vs  time 
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values.  As  expected,  the  single  Sears  covered  pool 
absorbs  more  of  the  solar  energy  than  the  other  covered, 
while  the  double  Bubble  covered  absorbs  the  least.  Both 
the  double  Sears  covered  pool  and  the  single  Bubble  covered 
absorb  the  same  amount  of  solar  energy. 

Figure  5-10  shows  the  total  convection  losses  from 
the  pool  to  its  surroundings.  The  convection  losses  are 
the  sum  of  side  surface  plus  top  surface  convection  losses. 
The  greatest  amount  of  convection  losses  is  registered 
by  the  solarly  heated  pool,  and  the  least  by  the  uncovered 
pool.  The  unheated  covered  pools  show  the  same  amount 
of  convection  losses.  There  is  a maximum  value  of 
convection  losses  in  June  of  1972  and  a minimum  in  December 
of  1972. 

Figure  5-11  shows  a similar  trend  of  sky  radiation 
losses  to  the  convection  losses.  The  only  difference 
between  the  two  trends  is  that  sky  radiation  curves  exhibit 
a plateau  from  April  of  1972  to  September  of  1972. 

After  examining  Figures  5-9,  10,  and  11,  it  is 
clear  that  the  solarly  heated  pool  should  have  the  highest 
bulk  water  temperature  because  the  amount  of  energy 
transferred  to  the  pool  is  the  highest  compared  to  any 
of  the  other  pools.  This  transferred  energy  reaches  a 
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maximum  value  in  the  summer  months  because  the  available 
solar  radiation  has  its  maximum  value.  Therefore,  the 
bulk  water  temperature  reaches  its  maximum  value  in  the 
summer  months.  The  unheated  covered  pools  absorb  less  of 
the  solar  radiation  because  of  their  covers,  and  therefore, 
their  bulk  water  temperatures  is  lower. 

From  Figures  5-9,  10,  and  11,  the  bulk  water 
temperature  for  the  uncovered  pool  should  have  been  at 
a higher  level  than  is  shown  in  Figure  5-8.  This  expec- 
tation is  due  to  the  higher  amounts  of  solar  radiation 
the  uncovered  pool  absorbs  and  the  smaller  amounts  of 
convection  and  sky  radiation  losses  as  compared  to  the 
unheated  covered  pools.  The  expectation  is  not  warranted 
because  the  uncovered  pool  has  an  additional  mode  of 
energy  transfer  that  none  of  the  other  pools  has;  that 
is,  energy  transferred  from  the  pool  by  evaporation  cool- 
ing. The  evaporation  cooling  and  the  absorbed  solar 
radiation  are  the  most  significant  modes  for  the  uncovered 
pool  as  shown  in  Figure  5-12.  In  contrast,  the  most 
significant  energy  modes  for  a covered  pool,  for  example, 
the  Sears  covered,  are  the  absorbed  solar  radiation  and 
convection  losses,  as  is  shown  in  Figure  5-13.  The  least 
significant  energy  mode  for  all  the  pools  is  the  sky 
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Figure  5-12. 


Energy  transferred  vs.  time  - uncovered  pool. 


Energy  Transferred  x 10  (BTU/Day) 
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TIME  (Month/Year) 

Figure  5-13.  Energy  transferred  vs.  time  - single  Sears 
covered  pool. 
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Temperature  Gains 

The  temperature  level  of  the  bulk  of  the  water  is 
one  measure  of  the  effectiveness  of  the  different  plastic 
covers  to  reduce  the  heat  losses.  Another  measure  of 
effectiveness  is  the  temperature  gain  of  the  water 
during  a specified  time  interval.  The  temperature  gain 
is  an  indication  of  how  well  the  swimming  pool  with  a 
particular  cover  can  absorb,  store,  and  radiate  back  to 
the  atmosphere  the  solar  radiation  it  receives  during  a 
time  interval . The  temperature  gain  can  be  found  by 
calculating  the  net  energy  transferred  to  the  swimming 
pool  and  dividing  it  by  the  product  of  the  total  mass  of 
the  water  times  the  specific  heat  of  the  water. 


AT 

GAIN 


Q 

net 

cm 


(5- 


The  net  energy  transferred  is  calculated  by  using  all 
the  theoretical  and  semi-empirical  equations  outlined 
in  Chapter  IV.  The  air  and  water  temperatures  can  be 
obtained  from  Table  5-1  and  Figure  5-8,  respectively. 

Since  the  objective  is  the  comparison  of  calculated 
values  of  temperature  gain  for  the  different  pools,  one 
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would  choose  a day  in  a month  in  which  the  temperature 
gain  value  will  be  most  pronounced.  This  is  accomplished 
by  choosing  the  day  with  the  highest  solar  irradiation. 

In  addition,  it  can  also  be  accomplished  by  choosing  as 
a time  interval  in  that  day  the  daytime  interval,  instead 
of  a 24-hour  period  of  time . The  experimentally  determined 
climatic  conditions  and  bulk  water  temperatures  measured 
in  the  morning  of  these  days  are  used  as  input  data  in 
the  calculation  of  the  temperature  gain  values . The 
calculated  values  of  the  temperature  gain  vs.  time  are 
presented  in  Figure  5-14. 

Monthly  Variations 

Now  that  the  variation  of  the  water  temperature  on 
an  annual  basis  was  presented  in  the  previous  section, 
it  would  be  of  interest  to  present  the  water  temperature 
variation  on  a monthly  basis. 

A time  period  of  seven  days  in  the  month  of  November 
1972  was  chosen.  The  objective  is  to  calculate  the  bulk 
water  temperature  of  the  various  pools  on  a continuous 
basis.  The  calculated  temperature  at  the  end  of  the 
daytime  will  be  used  as  input  to  calculating  the  water 
temperature  at  the  end  of  the  nighttime.  Other  data 
used  as  input  is  the  climatic  data  which  was  actually 
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Figure  5-14.  Calculated  temperature  gain  values  on  a high 
insolation  day  for  each  month. 
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measured  during  those  seven  days. 

Figure  5-15  presents  the  calculated  time  variation 
of  the  bulk  water  temperature  of  the  pools  as  the  climatic 
conditions  vary  from  day  to  day  during  the  selected 
seven  days  in  November  1972.  The  water  temperature  indi- 
cated on  a certain  date  is  the  calculated  morning  tempera- 
ture; the  next  water  temperature  indicated  is  the  calcu- 
lated evening  temperature  on  the  same  day  at  the  end  of 
the  day.  The  uncovered  pool  and  the  single  Sears  covered 
pool  are  assigned  the  water  temperature  values  actually 
measured  on  11/13/72.  The  other  covered  pools  are  arbi- 
trarily assigned  the  water  temperature  of  the  single 
Sears  covered  pool.  All  covered  pools  have  the  same 
water  temperature  as  a starting  point  in  the  calculations . 
As  it  can  be  seen  in  Figure  5-15,  it  takes  until  the 
evening  of  11/15/72  for  the  trend  of  the  water  temperature 
to  emerge.  That  is,  from  the  covered  pools  the  single 
Sears  covered  pool  maintains  the  lowest  water  temperature, 
while  the  double  Sears  covered  pool  maintains  the  second 
highest  water  temperature.  The  uncovered  pool  has  the 
minimum  temperature,  and  the  single  Sears  covered  pool 
which  is  heated  by  the  ten  solar  collectors  has  the 
maximum  temperature  as  shown  in  Figure  5-16  . While  all 
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Figure  5-15. 


Calculated  bulk  water  temperature  vs.  time. 
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Figure  5-16 . 


Calculated  water  temperature  vs.  time. 


81 


other  pools  show  a general  decrease  in  the  bulk  water 
temperature , the  heated  pool  by  the  solar  collectors  shows 
an  increase  in  its  temperature.  This  fact  establishes 
beyond  doubt  the  great  effectiveness  of  heating  the  swim- 
ming pool  with  the  solar  collectors  and  covering  its  sur- 
face with  a plastic  cover. 

An  interesting  fact  emerges  from  examining  Figure 
5-15;  the  uncovered  pool  registers  a net  water  temperature 
increase  during  the  24  hour  period  of  11/19/72  while 
all  other  pools  show  a decrease.  In  order  to  demonstrate 
this  fact  better,  the  temperature  gain  values  for  this 
time  period  are  plotted  in  Figure  5-17.  • It  should  be 
noted  that  a negative  temperature  gain  value  implies  a net 
water  temperature  decrease.  Indeed,  the  temperature 
gain  value  is  positive  for  the  uncovered  pool,  while  it 
is  negative  for  the  covered  ones  during  the  24  hour  period 
of  11/19/72.  This  can  be  explained  by  noting  that  the 
air  temperature  for  that  time  period  was  warmer  than  the 
water  temperature  of  the  uncovered  pool,  and  it  was  colder 
than  the  water  temperature  of  the  covered  pools.  But 
generally  the  uncovered  pool  has  the  greatest  temperature 
losses  and  the  single  Sears  covered  pool  which  is  heated 
by  the  solar  collectors  has  the  greatest  temperature 
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Figure  5-17. 


Calculated  temperature  gains  vs.  time. 
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gains . 

Daily  Variations 

The  temporal  variations  of  the  water  temperature 
of  the  pool  can  be  examined  in  a third  time  frame,  the 
daily  time  frame.  Given  the  climatic,  conditions  for  a 
certain  day  the  bulk  water  temperature  of  the  pools  under 
several  covers  is  predicted.  The  objective  is  to  predict 
how  the  bulk  water  temperature  varies  from  one  hour  to 
another  for  a certain  day. 

To  meet  that  objective  certain  days  in  the  year 
of  1972  were  selected  for  the  different  pool  covers  as 


follows : 

DAY 

COVER 

12/1/72 

No  cover 

12/12/72 

Single 

Sears  cover 

12/11/72 

Double 

Sears  cover 

12/18/72 

Single 

Bubble  cover 

12/19/72 

Double 

Bubble  cover 

12/13/72 

Single  Sears  cover  and  ten 
solar  heaters 

These  days  were  chosen  because  water  temperature  and 
climatic  data  were  taken  on  these  particular  days,  so  a 
comparison  between  predicted  and  experimental  values  can 
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be  accomplished.  Equation  (5-2)  was  utilized  again  in 
the  prediction  of  the  bulk  water  temperature.  The 
calculations  were  again  accomplished  by  using  a computer 
program.  The  calculated  water  temperature  at  the  end 
of  an  hour  was  used  as  an  initial  temperature  in  the  next 
hour.  The  first  initial  water  temperature  assigned  is 
the  experimentally  measured  one  for  a certain  day.  The 
calculated  time  variations  of  the  bulk  water  temperature 
for  the  various  covers  are  presented  in  Figures  5-18,  19, 
20  , 21,  22,  and  23.  Figure  5-18  and  Figure  5-23  present 

4 

quite  different  profiles  from  the  rest  of  the  figures. 

The  time  variation  of  the  water  temperature  in  the  un- 
covered pool  (Figure  5-18)  exhibits  a much  greater 
temperature  decrease  from  0600  to  1000  hours,  a slower 
temperature  increase  from  1000  to  1600  hours,  and  a 
greater  temperature  decrease  from  1600  to  2100  hours 
than  Figures  5-19,  20,  21,  and  22.  The  reason  is  that 
there  are  heat  transfer  losses  due  to  evaporation  cooling 
in  the  uncovered  pool  throughout  the  day,  and  there  are 
none  in  the  other  pools.  The  plastic  cover  impedes  the 
transfer  of  heat  out  of  the  water;  that  explains  the 
flatter  temperature  profiles  during  the  periods  of  0600 
to  1000  hours  and  1600  to  2100  hours  for  the  pools  with 
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(d)  oangeaedtuaj. 


Figure  5-18.  Hourly  variation  of  calculated  bulk  water  temperature  of  uncovered 
pool  on  12/1/72. 


86 


(d)  aangeuoduia,! 


Figure  5-19.  Hourly  variation  of  calculated  bulk  water  temperature  of  single  Sears 

covered  pool  on  12/12/72. 


87 


(d)  oungraoduioj, 


Figure  5-20.  Hourly  variation  of  calculated  bulk  v/ater  temperature  of  double  Sears 
covered  pool  on  12/11/72. 
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Figure  5-21.  Hourly  variation  of  calculated  bulk  water  temperature  of  single  bubble 

covered  pool  on  12/18/72. 
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igure  5-22.  Hourly  variation  of  calculated  bulk  water  temperature  of  double  bubble 

covered  pool  on  12/19/72. 


Water  Temperature  (F) 
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Figure  5-23.  Hourly  variation  of  calculated  bulk  water 

temperature  of  single  Sears  covered  and 
solarly  heated  pool  with  flat  plate 
collectors  on  12/13/72. 
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the  plastic  covers . The  pool  which  is  covered  with  a 
single  Sears  cover  and  is  heated  by  ten  solar  collectors 
exhibits  a much  greater  increase  in  its  time  variation 
of  the  water  temperature  (Figure  5-23) . This  is  the  result 
of  the  greater  amounts  of  absorbed  solar  energy  due  to 
the  pool's  greater  effective  surface  area.  Generally 
all  temperature  profiles  in  Figures  5-18  through  5-22 
show  a temperature  decrease  from  0600  to  1000  hours, 
an  increase  from  1000  to  1600  hours,  and  a decrease  from 
1600  to  2100  hours.  These  time  periods  demonstrate  the 
cooling  and  heating  periods  for  the  pools  during  a day. 


CHAPTER  VI 


EXPERIMENTAL  EQUIPMENT 

To  compare  the  temperature  of  the  water  of  a solarly 
heated  swimming  pool  to  the  temperature  of  an  unheated 
pool,  it  was  decided  to  place  two  identical  pools  next 
to  each  other.  Two  swimming  pools  fifteen  feet  in  diameter 
and  forty  inches  high  were  placed  above  the  ground  on  the 
west  side  of  the  Solar  House  of  the  University  of  Florida. 
The  capacity  of  each  pool  was  4400  U.S.  gallons  of  water. 
The  swimming  pools  were  equipped  with  circulating  pumps 
which  filtered  the  water  (Figure  6-1) . 

The  experimental  measurements  which  are  necessary 
to  support  any  theoretical  solution  obtained  by  Equation 
(4-41)  or  (4-18)  are  water  temperature,  dry  bulb  air 
temperature,  relative  humidity,  water  evaporation  rates, 
wind  speed,  and  insolation.  It  was  also  deemed  necessary 
to  measure  the  transmissivity  to  sun  light  of  the  various 
plastic  covers  used  for  the  swimming  pools. 

Experimental  Equipment  for 
Water  Temperature  Measurements 
There  were  two  reasons  for  measuring  the  water 
temperature  of  the  swimming  pools.  One  was  to  provide 
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Figure  6-1. 


Two  swimming  pools  at  University  of 
Florida  Solar  House . 
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a record  of  the  annual  temperature  variations  between  the 
solar ly  heated  and  unheated  swimming  pool ; for  this  reason 
the  bulk  water  temperature  was  monitored.  The  other  was 
to  determine  the  spatial  temperature  distribution  at 
different  time  intervals. 

The  bulk  water  temperature  was  measured  with  a 
mercury  glass  thermometer  with  an  accuracy  of  +0.5°F. 

The  description  of  the  procedure  to  measure  the  temperature 
appears  in  the  following  chapter.  Chapter  VII,  Experimental 
Procedure . 

The  spatial  and  hourly  variations  of  the  water 
temperature  were  measured  with  #24  gauge  copper-cons tantan 
thermocouples  which  were  recorded  by  three  Honeywell  multi- 
point recorders.  The  recorders  were  located  inside  the 
instrument  room  of  the  University  of  Florida  Solar 
House  (Figure  6-2) . The  thermocouples  were  placed  in  a 
crosslike  pattern  along  the  north-south  and  east-west 
axes  of  both  pools.  The  center  of  the  cross  pattern 
lay  in  the  center  of  the  swimming  pool.  There  were 
twenty-four  locations  in  the  cross  pattern  and  four  thermo- 
couples in  each  location  measuring  the  temperature  at 
four  different  depths.  Therefore,  in  each  pool,  there 
were  a total  number  of  eighty-four  thermocouples  which 
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were  measuring  the  water  temperature. 

In  the  construction  of  the  thermocouples,  care  was 
taken  to  let  the  thermocouples  have  the  least  possible 
electrical  contact  with  the  water  in  order  to  avoid 
generation  of  EMF  currents  due  to  a possible  electrolysis 
of  the  water  because  of  chlorine  and  minerals  present 
in  the  water.  The  thermocouples  were  tied  to  a wooden 
dowel  and  inserted  in  a water  resistant  plastic  tube; 
the  ends  of  the  tube  were  sealed  with  a water  resistant 
sealant  (Figure  6-3) . The  locations  of  the  thermocouples 
in  the  north  and  south  pools  are  shown  in  Figures  6-4,  and 
5. 

The  calibration  of  the  thermocouples  was  accomplished 
by  calibrating  the  Honeywell  recorders  to  which  they 
were  connected.  The  voltage  output  of  the  thermocouples 
which  was  monitored  on  the  recorders  was  plotted  versus 
a known  water  temperature  range.  Several  thermocouples 
were  first  inserted  in  a container  with  ice  and  water 
(32°F)  and  then  the  mixture  was  heated  to  the  boiling 
point  of  the  water  (212°F) . The  plots  of  thermocouple 
voltage  output  in  a known  water  temperature  range  became 
the  calibration  curves  of  the  Honeywell  recorders, 
and  also,  of  the  thermocouples. 
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Figure  6-2.  Honeywell  multi -point  recorders  used  to 

measure  the  temperatures . 


Figure  6-3 


Location  of  Thermocouples 
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'igure  6-4.  Location  of  thermocouples  in  north  swimming  pool. 


South  Side  North  Side 

Figure  6-5.  Location  of  thermocouples  in  south  swimming  pool. 
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Experimental  Equipment  for  Air  Temperature 
and  Relative  Humidity  Measurements 
The  dry  bulb  air  temperature  and  relative  humidity 
were  both  measured  by  the  same  instrument,  a Belfort 
mechanical  hygrometer  (Figure  6-6) . As  can  be  seen  from 
a typical  weekly  chart  (Figure  6-7) , the  relative  humidity 
and  air  temperature  can  be  read  on  an  hourly  basis. 

The  accuracy  for  the  relative  humidity  is  ±1%,  and  for 
the  air  temperature  is  +1.0F.  The  instrument  was  housed 
inside  a standard  United  States  Weather  Bureau  shelter 
which  was  located  near  the  swimming  pools  (Figure  6-8)  . 

The  hygrometer  was  calibrated  with  a sling  psychrometer . 

Experimental  Equipment  for 
Water  Evaporation  Rates  Measurements 
Early  it  was  attempted  to  measure  water  evaporation 
rates  by  measuring  the  depth  of  the  water  in  the  swimming 
pools  once  a day.  The  data  did  not  prove  conclusive 
due  to  certain  leaks  in  the  filter  loop  of  the  swimming 
pools.  An  alternate  method  was  adopted  by  using 
Equation  (4-38)  . The  constant  C'  in  Equation  (4-38)  was 
determined  experimentally,  as  will  be  shown  in  Chapter 


VIII,  Data  and  Results. 
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Figure  6-6.  Belfort  Mechanical  Hygrometer. 


Figure  6-7.  Typical  weekly  chart. 
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Experimental  Equipment  for 
Wind  Speed  Data  Measurements 
Wind  speed  data  has  been  taken  by  the  Agricultural 
Experiment  Station  of  the  University  of  Florida  at  a 
location  which  is  two  miles  west  of  the  Solar  House. 

Because  of  the  close  proximity  of  the  two  locations,  their 
wind  speed  data  were  utilized  in  this  investigation. 

The  error  in  the  measurements  is  estimated  to  be  ±5%. 

Experimental  Equipment  for 
Solar  Insolation  Measurements 
Insolation  measurements  were  made  with  a Matrix- 
Sol-a-meter,  Mark  XVIII  (Figure  6-9) . The  calibration 
of  this  meter  was  performed  by  its  manufacturer  which 
was  accomplished  by  comparing  it  to  a standard  pyranometer 
within  +1%.  It  has  an  array  of  silicon  cells,  and  it  can 
provide  an  instantaneous  measurement,  as  well  as  the 
total  daily  insolation. 

Experimental  Equipment  for 
Transmissivity  Measurements 

The  transmissivity  measurements  of  the  various  plastic 
covers  were  conducted  at  the  Solar  Laboratory  of  the 
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Figure  6-9.  Sol-A-Meter  Mark  XVIII. 
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University  of  Florida  utilizing  a "black  box."  This  will 
be  described  in  more  detail  in  Chapter  VIII,  Data  and 


Results . 


CHAPTER  VII 

EXPERIMENTAL  PROCEDURE 

The  measurement  of  the  bulk  water  temperature  was 
made  twice  a day  from  April  1972  to  April  1973  in  both 
swimming  pools  to  provide  a record  of  the  annual  water 
temperature  variations.  One  swimming  pool  was  always 
uncovered  and  was  not  heated  by  any  solar  collectors; 
this  pool  will  be  referred  to  as  the  unheated  pool. 

The  other  pool  was  always  covered  with  a plastic  cover, 
and  sometimes  was  heated  with  solar  collectors;  this 
pool  will  be  referred  to  as  the  solarly  heated  pool. 

Each  day  after  sunset,  the  filtering  pumps  were 
turned  on,  and  were  in  operation  ail  night.  Just  before 
the  sun  rays  struck  the  top  surface  of  the  swimming 
pools,  the  filtering  pumps  were  turned  off.  The  purpose 
of  the  filtering  pumps  was,  other  than  the  filtering  of 
the  water,  the  thorough  mixing  of  the  water  so  that 
accurate  measurements  of  the  bulk  water  temperature  could 
be  made.  The  bulk  water  temperature  was  measured 
twice  a day  in  both  swimming  pools,  early  in  the  morning 
and  late  in  the  afternoon.  Occasionally,  the  water 
temperature  was  measured  in  several  locations  in  the  same 
pool  to  verify  that  the  measured  temperatures  were  truly 
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the  bulk  water  temperatures;  no  differences  ever  arose. 
The  bulk  water  temperatures  were  measured  with  a mercury 
glass  thermometer.  One  of  the  reasons  for  not  measuring 
bulk  water  temperatures  with  the  thermocouples  was  the 
ease  and  short  calibration  time  which  was  required  to 
calibrate  a single  mercury  thermometer  as  compared  to  the 
calibration  of  several  thermocouples. 

During  the  same  time  interval  of  April  1972  to  April 
1973,  continuous  measurements  were  made  of  the  relative 
humidity  and  dry  bulb  temperature  of  the  ambient  air. 

If  a need  arose  to  know  the  instantaneous  values  for  the 
relative  humidity  and  air  temperature,  one  had  to  simply 
read  them  off  the  Belfort  hygrometer  chart.  If  the  need 
arose  to  know  their  average  values  during  a certain  time 
interval,  one  had  to  obtain  the  chart  of  that  time  inter- 
val, and  with  the  aid  of  a planimeter,  to  calculate 
them. 

The  same  procedure  was  followed  with  the  insolation 
measurements.  The  Sol-a-meter  could  provide  the 
instantaneous  insolation  value,  or  the  value  of  the 
total  insolation  for  one  day. 

There  were  no  instantaneous  data  for  the  wind  speed. 
As  it  was  described  in  Chapter  V,  the  daily  average 
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values  for  the  period  of  April  1972  to  April  1973  were 
obtained  from  the  Agricultural  Experiment  Station  of  the 
University  of  Florida. 

Spatial  water  temperatures  were  measured  with  the 
copper-constantan  thermocouples  for  six  days  in  December 
of  1972.  Each  day  the  solarly  heated  pool  had  a different 
plastic  cover,  and,  or,  was  heated  by  the  solar  collec- 
tors . The  night  before  the  spatial  temperatures  were 
taken,  the  filtering  pumps  were  not  turned  on.  The  reason 
was  to  map  the  initial  spatial  distribution  of  the 
temperatures.  The  first  measurements  started  at  0600 
hours.  Since  sunrise  time  was  not  until  about  0800  hours 
on  December  1,  1972,  the  spatial  measurements  provided 
a temperature  record  of  the  cooling  of  the  water  by  the 
ambient  air.  The  spatial  measurements  were  taken  in 
hourly  increments,  thus  providing  a detailed  record  of  the 
heating  or  cooling  of  the  water.  The  last  measurement  of 
the  spatial  temperatures  was  made  at  2100  hours.  Imme- 
diately after  the  spatial  temperatures  were  recorded, 
the  average  values  for  the  relative  humidity  and  dry 
bulb  temperature  of  the  air  were  read  off  the  Belfort 
hygrometer.  The  Sol-a-meter  provided  the  cumulative 
insolation  during  each  hour. 


CHAPTER  VIII 


EXPERIMENTAL  RESULTS  AND  DATA 

For  the  successful  application  of  the  theoretical 
equations  which  were  developed  in  Chapter  IV,  certain 
experimental  information  was  required.  If  the  experi- 
mental information  was  not  readily  available  in  the  litera 
ture,  or  it  was  in  a form  which  was  not  directly  applica- 
ble to  this  investigation,  then  experiments  were  planned 
and  executed.  The  results  of  these  experiments  are 
presented  in  this  Chapter  under  Experimental  Results. 
Finally,  the  experimental  data  of  the  water  temperatures 
of  the  uncovered  and  covered  swimming  pools,  along  with 
all  the  measured  climatic  data  are  presented  under 
Experimental  Data. 

Experimental  Results 

To  understand  the  physical  phenomena  taking  place  in 
the  solarly  heated  and  unheated  swimming  pools,  the 
following  parameters  were  experimentally  investigated: 

(a)  Evaporation  rates  of  water 

(b)  Conduction  losses  from  the  swimming  pool  to 
the  ground 

(c)  Transmissivity  measurements  of  various  plastic 
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covers 

(d)  Absorptivity  measurements  of  water 

(e)  Wind  speed  data 

(f)  Surface  conductance  data 

(a)  Evaporation  Rates  of  Water 

The  need  for  the  determination  of  the  evaporation 

rate  of  water  is  twofold.  As  was  shown  in  Chapter  IV, 

there  is  disagreement  among  researchers  as  to  the  value 

of  the  constant  C*  in  Equation  (4-38) . It  was  decided 

that,  through  a simple  experimental  procedure,  the  value 

of  C'  could  be  calculated.  The  second  need  arose  from 

the  leaking  of  water  from  the  covered  swimming  pool  which 

made  it  impossible  to  calculate  any  evaporation  rate  of 

water  by  measuring  the  water  level. 

One  oil  drum  was  cut  radially  in  two;  the  two 

cylindrical  ends  had  a height  of  14  inches  and  a diameter 

of  22.5  inches.  They  were  filled  with  water  and  placed 

on  the  ground  on  the  west  side  of  the  two  swimming  pools. 

One  container  was  covered  with  a plastic  cover.  The 

temperatures  of  water,  dry  bulb  air  temperature,  relative 

humidity,  and  water  levels  in  both  cylindrical  containers 

were  taken  twice  daily.  The  driving  force  in  the 

evaporation  process  Ap  = p -p  in  Equation  (4-38)  can  be 

w a 
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calculated  by  first  calculating  p from  the  measured  water 

w 

temperature  in  the  container.  Next  p can  be  calculated 

a 

by  the  following  equation 


p = <f>p 
a w,s 


(8-1) 


The  evaporation  rate  G in  Equation  (4-38)  can  be  calculated 
by 


G = pL' 


(8-2) 


The  results  of  the  measurements  are  presented  in 
Table  8-1  and  plotted  in  Figure  8-1  as  evaporation  rate 
G vs.  driving  force  Ap . The  slope  of  the  best  curve 
through  the  data  points  was  calculated  by  the  method  of 
least  squares  and  found  to  be  equal  to  0.241.  Therefore, 
the  constant  C'  in  Equation  (4-38)  is  equal  to  0.241. 

The  scatter  in  the  data  in  Figure  8-1  could  be  attributed 
to  possible  experimental  errors  from  measuring  the  water 
and  air  temperature,  the  relative  humidity,  and  the  length 
L' . No  evaporation  took  place  in  the  covered  container. 

The  measured  evaporation  data  was  compared  against 
evaporation  data  measured  by  the  Agronomy  Department 
of  the  University  of  Florida  at  a location  two  miles 
west  of  the  test  site.  Their  data  is  presented  in  Table 


EVAPORATION  DATA  TAKEN  AT  THE  SOLAR  HOUSE 
OF  UNIVERSITY  OF  FLORIDA 
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Figure  8-1.  Experimental  data  of  water  evaporation  rate  vs 

press ure  di f f e ren ce  . 
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8-2  and  plotted  also  in  Figure  8-1.  There  is  good  agree- 
ment between  the  two  independently  taken  evaporation  data. 

The  constant  C' =0.241  which  was  calculated  from  the 
data  in  Figure  8-1  was  then  compared  against  the  constant 
C'  as  calculated  from  various  sources  in  the  literature. 

The  results  are  presented  in  Table  8-3. 

As  can  be  seen  from  Table  8-3,  most  of  the  calculated 
values  are  higher  than  C' =0.241,  even  when  the  wind  speed 
is  1 mph.  Since  the  measured  average  speed  at  the  time 
of  the  data  taken  at  the  test  site  was  about  seven  mph, 
it  is  evident  that  there  is  quite  a disagreement.  This  can 
be  attributed  to  differences  in  the  experimental  arrange- 
ment, water  temperature  range,  and  wind  speed  range. 

(b)  Conduction  Losses  to  the  Ground 

The  following  experimental  procedure  was  followed 
to  measure  the  relative  importance  of  conduction  losses  to 
the  other  modes  of  heat  transfer  and  mass  transfer. 

A copper-constantan  thermocouple  was  located  six 
inches  deep  in  the  ground  under  thermocouple  #49  (see 
Figure  6-4  for  location  of  thermocouple  #49) . A second 
copper-constantan  thermocouple  was  placed  five  inches 
deep  in  the  ground  fifteen  feet  away  from  both  swimming 
pools.  All  thermocouples  were  connected  to  a multi-point 


EVAPORATION  DATA  MEASURED  BY  AGRONOMY  DEPARTMENT 
OF  UNIVERSITY  OF  FLORIDA 
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TABLE  8-3 

VALUE  OF  CONSTANT  C' 
IN  EQUATION  (4-38) 


Literature 

Source 

Wind  Speed  = 
1 0 mph 

Wind  Speed  = 
1 mph 

‘[32] 

1.025 

0.324 

[32] 

1.062 

0.312 

[33] 

0.905 

0 . 26 

[34] 

0.966 

0.392 

[35] 

0.713 

0.245 

[31] 

1.23 

0.185 

[8] 

0.80 

~0 . 20 
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Honeywell  Recorder.  Data  were  taken  continuously  for  25 
hours  on  July  2 and  3,  19  73.  The  data  were  recorded 
in  Table  8-4. 

The  heat  losses  due  to  conduction  were  calculated 
by 


q 

cond 

loss 


bottom 


7T7T 


7^ 


TJT 

I 


water 


♦ *J.  + 


wall 


ground 


r 'i 

T -T 
#49  ground 

6"  under 

#49 


(197.3939) 


f 

T 

#49 


ground 
6"  under 
#49 


(8-3) 


The  thermal  conductivity  of  the  ground  at  the  Solar 
House  of  the  University  of  Florida  was  measured  by  Dickert 
[36]  to  be  k=0.82  BTU/(hr.  ft.  F) . The  results  of  the 
calculations  and  the  data  are  plotted  in  Figure  8-2  and 
8-3. 

These  figures  show  that  the  conduction  losses  during 
24  hours  were  2249.9  BTU  from  the  water  to  the  ground. 

The  conduction  losses  proved  to  be  two  orders  of  magnitude 
smaller  than  the  evaporation  or  convection  losses , as 


shown  in  Figure  8-4 . 
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TABLE  8-4 
CONDUCTION  DATA 


Date 

Time 

T 4- 

water 

(F) 

(T/C#  49 ) 

Tground 
6 " unde  r 
#49  (F) 

ground 
15 1 away 
from 
#49  (F) 

(Tw-Tgr) 

(F) 

^conduction 
loss,  BTU 

Air 

T.emp 

(F) 

7/2 

1300 

. 90.25 

89.35 

84.9 

0.9 

177.6 

91 

1400 

91.1 

89.8 

86.2 

1.3 

256.6' 

92 

1500 

91.6 

91.2 

88.0 

0.4 

78.9 

92 

1600 

92.75 

91.9 

89.05 

0.85 

167.8 

92 

1700 

93.5 

92.45 

90.1 

1.05 

207.3 

91 

1800 

94.1 

93.25 

90.85 

0.85 

167.8 

91 

1900 

95  .0 

93.8 

91.1 

1.2 

236.9 

89 

2000 

94.85 

93.95 

90.7 

0.9 

177.6 

87 

2100 

93.95 

93.5 

89.8 

0.45 

88.8 

83 

2200 

93.5 

92.9 

88.9 

0.6 

118.4 

81.5 

2300 

92.9 

92.3 

88.15 

0.6 

118.4 

78 

2400 

92.3 

92.0 

87.7 

0.3 

59.2 

77 

7/3 

0100 

91.9 

91.7 

86.8 

0.2 

39.5 

76 

0200 

91.3 

91.2 

86.5 

0.1 

19.7 

74 

0300 

91.1 

91.0 

86.1 

0.1 

19.7 

73 

0400 

90.7 

90.4 

85.7 

0.3 

59.2 

72.5 

0500 

90.25 

90.25 

85.5 

0.0 

0.0 

72 

0600 

89.5 

89.5 

85.2 

0.0 

0.0 

72 

0700 

89.35 

89.35 

84.45 

0.0 

0.0 

71.5 

0800 

88.75 

83.9 

84.3 

-0.15 

-29.6 

71 

0900 

88.75 

88.75 

84.15 

0.0 

0.0 

76 

1000 

89.00 

88.60 

84.07 

0.4 

78.9 

80.5 

1100 

89.2 

88.60 

83.55 

0.6 

118.4 

84 

1200 

89.8 

89.35 

84.3 

0.45 

88.8 

87 

1300 
Avg . 

91.1 

91.475 

89.8 

91.0 

85.5 

86.92 

1.3 

256.6 

90 

81.417 

qtotal 

2249.9  BTU/Day  Cond. 

loss  to 

ground 
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(j)  oangujroduiox 


Figure  8-2.  Conduction  data 
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Heat  Losses  x 10  (BTU/Day) 
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7/2/73  7/3/73 

TIME  (Hours) 

Figure  8-4.  Heat  losses  vs.  time. 
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(c)  Transmissivity  Measurements  of  Various 

Plastic  Covers 

Various  plastics  were  used  as  a cover  for  the 
swimming  pool  to  eliminate  heat  losses  due  to  evaporation 
of  water.  Since  no  data  could  be  found  as  to  the  plastics' 
transmissivities  to  solar  radiation,  experiments  were 
undertaken  to  measure  them. 

Two  plastics  were  used.  Both  are  made  of  polyethy- 
lene. The  first  one  will  be  referred  to  as  "Sears 
Plastic"  because  it  is  readily  available  from  a Sears 
Department  Store,  the  second  one  as  "Bubble  Plastic"  because 
it  has  3/4  inch  air  bubbles  (Figure  5-7) . The  Sears 
plastic  is  2 mils  thick,  while  the  Bubble  is  4 mils 
thick . 

The  transmissivity  tests  were  conducted  in  the 
Solar  Energy  Laboratory  of  the  University  of  Florida. 

The  instrument  used  was  a "Black  Box"  (Figure  8-5) . The 
plastic  to  be  measured  is  mounted  on  the  "window"  of 
the  "Black  Box."  With  the  aid  of  a pyranometer  the 
transmitted  energy  through  the  plastic  is  measured 
(Figure  8-6) . Also  the  total  available  solar  radiation 
is  measured  with  the  pyranometer.  The  ratio  of  the 
transmitted  to  the  total  incident  solar  energy  gives  the 
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Figure  8-6 


In  operation  "Black  Box 
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transmissivity  of  the  plastic  for  a particular  incidence 
angle  of  the  sun.  (The  incidence  angle  of  the  sun  is 
the  angle  between  the  sun's  rays  and  the  normal  to  the 
test  surface.) 

The  incidence  angle  of  the  sun  is  calculated  by 

<p  = cos  1 [cos  L cos  h cos  6 + sin  L sin  6]  (8-4) 

Figure  8-7  shows  the  variation  of  the  incidence  angle 
of  the  sun  with  the  time  of  the  day  and  the  season  of 
the  year  for  a 30  degrees  north  latitude. 

It  is  evident  that  the  amount  of  energy  transmitted 
through  the  plastic  depends  very  strongly  on  the  incidence 
angle  of  the  sun,  and  since  the  incidence  angle  varies 
considerably  with  time  (Figure  8-7) , there  is  a need  for 
an  "apparent  transmissivity"  of  the  plastic  to  be  defined. 
The  apparent  transmissivity  is  the  value  of  the  trans- 
mitted energy  through  the  plastic  at  a certain  time  of  the 
day  and  certain  seasons  of  the  year  over  the  incident 
solar  energy. 

The  Sears  Plastic  and  the  Bubble  Plastic  were 
measured  in  two  arrangements;  first,  a single  layer  of 
plastic,  and  second,  two  layers  of  plastic  with  a 2 inch 
air  space  between  them.  The  data  obtained  is  plotted  in 


Incidence  Angle  of  Sun  (Degrees) 
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Incidence  angle  of  sun  vs.  solar  tine  for 
Gainesville,  Florida 


Figure  8-7. 
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Figure  8-8  and  Figure  8-9.  It  is  apparent  that  for  both 
plastics  the  apparent  transmissivity  decreases  with  in- 
creasing incidence  angle,  and  decreases  with  increasing 
the  number  of  layers  of  plastic.  The  Sears  Plastic 
seems  to  have  a better  ability  to  transmit  energy.  In 
both  figures,  the  apparent  transmissivity  of  plate  glass 
is  also  drawn.  It  seems  that  both  plastics  are  inferior 
to  glass  concerning  transmission  of  solar  energy. 

Figures  8-8  and  8-9  show  how  the  apparent  trans- 
missivity of  the  two  plastics  varies  with  the  incidence 
angle  on  a daily  basis.  But  since  the  incidence  angle 
of  the  sun  varies  with  the  season  of  the  year,  it  follows 
that  the  apparent  transmissivities  will  also  vary  with  the 
season  of  the  year. 

Specifically,  there  is  the  need  to  know  the  average 
value  of  the  apparent  transmissivity  of  a plastic  (x) 
during  a certain  day  of  a month  of  the  year.  During  this 
day  of  that  month,  there  is  an  average  value  of  solar 
energy  which  is  transmitted  through  the  plastic;  this 
transmitted  energy  is 


n 

Q,  = E x . H 

transm.  . , i su, 
•i=l  l 


(8-5) 


Also  the  total  available  solar  energy  incident  on  the 


Apparent  Transmissivity 
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Incidence  Angle  of  Sun  (Degrees) 

Figure  8-8.  Apparent  transmissivity  of  Sears  plastic  vs. 
incidence  angle  of  sun. 


Apparent  Transmissivity 
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Incidence  Angle  of  Sun  (Degrees) 

Figure  8-9.  Apparent  transmissivity  of  bubble  plastic  vs. 
angle  of  incidence  of  sun. 
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plastic  during  this  day  of  the  month  is 


n 

H = E H 
su  . , su. 
1=1  1 


(8-6) 


The  average  value  of  the  apparent  transmissivity  during 
a day  of  a certain  month  is  then 


T 


transm. 


H 

su 


n 


E x . H 
1 su 

n 


E 

i=l 


H 

su . 

l 


(8-7) 


The  incident  solar  energy  H 
that  is 


Sll 


has  two  components , 


H 

su 


I , . * + I , . 
dif  dir 


But 


(8-8) 


I,.  = I,  * cos<j)  • cf  . (8-9) 

dir  dn 


There  are  published  data  for  I,._  and  In  for  a certain 

di  t dn 

day  of  each  month.  Then  by  using  Equations  (8-8)  and 


(8-9)  one  can  calculate  the  time  variation  of  H . These 

su 

calculations  are  shown  in  Appendix  D [37]. 


Once  the  available  solar  energy  (H  ) has  been 

su 

calculated  for  a certain  day  of  the  month,  the  transmitted 
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energy  needs  to  be  calculated  for  the  same  day  of  the 
month.  This  can  be  accomplished  by  calculating  the 
product  of  x H for  each  hour  of  the  day  and  then 

•L  b u , 

1 

summing  up  all  the  products  for  all  the  hours  of  the 
day  in  Equation  (8-5) 5 


Q 


transm. 


n 

l 

i=l 


T.  H 

1 su.  • 
1 


(8-5) 


The  hourly  variations  of  the  transmissivity  are 
taken  from  Figures  (8—8)  and  (8-9) . Thus  the  value  of 
the  average  apparent  transmissivity  is  calculated  for  a 
certain  day  of  a month  by  using  Equation  (8-7) . The  same 
procedure  can  be  repeated  for  the  same  day  for  the  rest 
of  the  months  of  the  year.  These  calculations  are  shown 
in  Appendix  E and  F where  the  apparent  transmissivities 
of  a single  layer  of  a Sears  plastic  (x  ) and  of  a 

O -L 

Bubble  Plastic  (x  ) r and  also  of  a double  layer  of  a 
Sears  Plastic  (xcc_)  and  of  a Bubble  Plastic  (x  ) 
are  calculated  as  a function  of  time . These  results 
are  plotted  in  Figure  8-10  as  apparent  transmissivity  of 
plastic  covers  vs.  time.  All  four  curves  in  Figure  8-10 
exhibit  a maximum  point  near  the  middle  of  June  since 
then  the  incidence  angle  of  the  sun  has  its  minimum  value 
(see  Figure  8-7) . One  layer  of  Sears  Plastic  allows  the 


Apparent  Transmissivity 
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Figure  8-10. 


Apparent  transmissivity  vs.  time. 
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greatest  amount  of  solar  energy  to  go  through,  while  two 
layers  of  Bubble  plastic  the  least.  There  is  a very  small 
difference  in  the  variation  of  the  transmissivity  near 
June  between  one  layer  of  Bubble  Plastic  and  two  layers 
of  Sears  Plastic.  This  is  expected  because  of  the  same 
material  and  the  thickness  of  the  two  plastics. 

(d)  Absorptivity  Measurements  of  Water 

The  absorptivity  of  water  is  needed  in  calculating 
the  amount  of  energy  absorbed  by  the  water.  Obviously, 
the  value  of  the  absorptivity  depends  on  the  incidence 
angle  of  the  sun.  Therefore,  there  is  a similar  need  to 
define  an  apparent  absorptivity  and  calculate  its 
variation  with  time.  Such  a variation  is  obtained  from 
experimental  data  of  absorptivity  of  water  versus  incidence 
angle  of  sun  and  using  Equations  (8-6)  , (8-8)  , and  (8-9) 

to  obtain  the  available  solar  energy. 

The  measured  variation  of  the  apparent  absorptivity 
of  water  as  a function  of  the  incidence  angle  of  the  sun 
is  shown  in  Figure  8-11.  The  apparent  absorptivity 
decreases  with  increasing  incidence  angle  of  the  sun, 
as  expected.  The  absorptivity  data  were  obtained  by 
measuring  the  reflectivity  of  the  water  at  different  times 
for  several  days  with  a single  cell  matrix  sol-a-meter 


Apparent  Absorptivity  of  Water 
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Incidence  Angle  of  Sun  (Degrees) 


Figure  8-11. 


Measured  apparent  absorptivity  of 
water  vs.  incidence  angle  of  sun. 
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(Figure  8 12) . If  the  body  of  the  water  in  the  swimming 
pool  can  be  regarded  as  an  opaque  body,  then 


a = 1 -p'  (8-10) 

To  find  the  variation  of  the  apparent  absorptivity  of 
water  with  time,  calculations  are  shown  in  Appendix  E. 

The  results  are  plotted  in  Figure  8-13.  The  apparent 
absorptivity  reaches  a maximum  point  in  June. 

(e)  Wind  Speed  Data 

No  additional  supportive  experimental  data  were 
taken.  However,  two  other  experimental  sources  were  util- 
ized to  provide  wind  speed  data  and  surface  conductance 
data.  The  Agronomy  Department  of  the  University  of  Florida 
provided  data  on  wind  speed  versus  time;  their  results 
have  been  plotted  in  Figure  8-14  for  the  year  of  1972  and 
part  of  1973.  Data  for  individual  days  are  available 
in  Agronomy  Research  Report  AG  73-2. 

(f)  Surface  Conductance  Data 


The  surface  conductance  of  stucco  and  glass  vs.  wind 
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Figure  8-12. 


Single  Cell  Sol-A-Meter. 


Apparent  Absorptivity 
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Figure  8-13.  Apparent  absorptivity  of  water  vs.  time. 
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speed  have  been  plotted  in  Figure  4-6.  The  two  curves 
represent  the  highest  and  the  lowest  of  the  curves  from 
experimental  data  in  [26]  . These  curves  have  been  replotted 

after  removing  the  radiant  portion  of  the  surface  conduc- 
tance . 


Experimental  Data 

Experimental  data  is  needed  to  support  the  analytical 
results.  Therefore,  bulk  water  temperature  data  were 
measured  twice  every  day  for  a period  of  twelve  months. 

Also  data  were  taken  by  measuring  the  water  temperatures 
at  the  eighty— four  locations  within  the  pool. 

The  bulk  water  temperature  data  are  presented  in  Table 
8-5  for  each  day  of  the  month  and  for  every  month  in  the 
time  period  of  April  1972  to  March  1973.  The  values  of 
the  temperatures  were  obtained  by  forming  the  arithmetic 
average  value  of  the  morning  and  afternoon  readings  for 
each  day.  Unfortunately  there  are  some  days  during  which 
no  temperature  data  were  taken;  this  was  due  to  necessary 
cleaning  of  the  water  and  the  pools,  and  to  periodic 
calibrations  of  the  measuring  instruments.  The  monthly 
mean  bulk  water  temperature  for  each  month  is  plotted 
in  Figure  8-15 . The  lowest  bulk  water  temperature  was 
measured  in  February  1973  and  the  highest  in  July  1972  for 


BULK  WATER  TEMPERATURE  (F) 
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Figure  8-15. 


Measured  bulk  water  temperature  vs.  time. 
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both  the  uncovered  and  the  single  Sears  covered  pool. 

The  reported  measured  bulk  water  temperatures  were 
the  actual  measured  temperatures.  However,  in  case 
of  rain,  the  reported  bulk  temperature  was  raised  to  com- 
pensate for  the  cooling  effect  of  the  pool  water  mixing 
with  the  colder  rain  water. 

For  purposes  of  comparing  the  theoretically  calcu- 
lated with  the  measured  values  for  the  time  period  analyzed 
in  the  Monthly  Variations,  Chapter  IV,  Figure  8-16  was 
prepared;  it  shows  the  variation  of  the  bulk  water  tempera- 
ture for  the  time  period  of  November  13,  1972,  to  November 
20,  1972. 

Daily  water  measurements  v/ere  taken  from  December  1, 
1972,  to  December  19,  1972.  Water  temperature  and  climatic 
data  was  taken  each  hour.  Starting  at  0600  hours  and  end- 
ing at  2100  hours,  data  was  taken  with  a different  cover 
for  each  day  of  measurements.  A sample  of  daily  water 
temperature  measurements  for  the  uncovered  pool  appears 
in  Appendix  G.  The  temperature  at  the  geometric  center 
of  each  covered  pool  is  plotted  in  Figure  8-17  through 
Figure  8—22;  the  value  of  the  water  temperature  is  ob- 
tained by  forming  the  arithmetic  average  value  of  the 
temperature  readings  of  thermocouples  #122  and  #123 


Temperature  (F) 


O Bulk  Water  Temp. -Single  Sears  Covered 
□ Bulk  Water  Temp.-  Uncovered 
O Air  Temperature 


1972 

TIME  (Month/Day) 


Figure  8-16  . Measured  bulk  water  temperature  vs . 

time . 
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Figure  3-17.  Hourly  variation  of  measured  water  temperature  at  the  center  of 

the  uncovered  pool  on  12/1/72. 
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(d)  oangeuaduioj, 


single  bubble  covered  pool  on  12/18/72 
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(d)  ajtnqBJGduioj,  jro:p?M 


double  bubble  covered  pool  on  12/19/72. 


Water  Temperature  (F) 
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TIME  (hours) 

Figure  8-22.  Hourly  variation  of  measured  water 

temperature  at  the  center  of  a single 
Sears  covered  and  solarly  heated  by 
collectors  pool  on  12/13/72 
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(for  location  of  the  thermocouples  see  Figure  6-5) . The 
common  characteristics  in  all  Figures  8-17  to  8-22  are  the 
similar  cooling  and  heating  periods  of  the  water. 

The  measured  spatial  water  temperatures  are  shown 
in  Figures  8-23  and  8-24.  The  plotted  temperatures  were 
obtained  by  averaging  all  thermocouple  readings  at  each 
height.  Measurements  were  taken  on  12/1/72  from  0600  hours 
to  2100  hours  in  the  uncovered  pool.  In  Figure  8-25 
are  shown  the  water  temperature  variations  radially  at  the 
mid-height  of  the  pool;  only  temperatures  from  the  eastern 
portion  of  the  cross-like  pattern  of  the  thermocouples 
have  been  plotted. 
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Figure  8-23.  Measured  spatial  water  temperature  variations  on  12/1/72 
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Figure  3-24.  Measured  spatial  water  temperature  variations  on  12/1/72 
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Radius  (feet) 

Figure  8-25.  Water  temperature  variation  radially 

at  the  mid-height  of  the  pool  on 
12/1/72  - uncovered  pool 


CHAPTER  IX 


DISCUSSION  AND  EVALUATION  OF  RESULTS 


Spatial  Water  Temperature  Variations 


A basic  question  is  which  one  of  the  two  distinct 
flow  regimes,  as  described  by  Eckert  and  Carlson  [17] 
and  Elder  [18,  19],  is  present  in  a solarly  heated  or  in 
an  unheated  swimming  pool.  The  assumptions  which  were 
made  in  the  analysis.  Chapter  V,  were  true  for  the  second 
flow  regime,  where  the  core  temperature  is  constant  in  the 
horizontal  direction  and  there  is ‘a  constant  vertical 
temperature  gradient,  as  described  in  references  [17,  18, 
19,  20]. 

The  calculation  of  the  Rayleigh  number  confirms 
the  presence  of  the  second  flow  regime.  For  a 60°F  water 
temperature,  and  a water  temperature  difference  of  1°F 
we  have 


Pr  = 8.03 


Gr 


eg  AT  l 


2 

V 


0.85x10  4x32.17xlx(3.5) 3 
(1 . 22xl0_5 ) 2 


7.87  x 108. 


Therefore 


8 9 

Ra  = Gr  • Pr  = 7.87x10  x8.03  = 6.32x10 
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g 

The  value  of  Ra  = 6.32  x 10  is  much  greater  than  1000 
which  was  the  upper  limit  of  the  first  flow  regime, 
reference  [18,  19] . Equation  (4—18)  was  developed  to  calcu 
late  the  spatial  water  temperatures  inside  the  core  of  the 
torrus.  The  calculated  core  temperatures  indicate  a 
definite  vertical  temperature  gradient  which  is  not 
necessarily  constant  (Figures  5-4  and  5-5) . The  measured 
temperatures  in  the  vertical  direction  also  confirm  the 
calculated  temperatures  (Figures  8-23  and  8-24).  Also, 
in  Figure  8-25,  it  is  shown  that  there  is  a very  small 
temperature  variation  in  the  horizontal  direction.  In 
concluding  then,  the  analytical  and  experimental  results 
of  this  investigation  confirm  the  presence  of  the  second 
flow  regime  in  a heated  or  cooled  swimming  pool. 

There  is  very  good  agreement  of  calculated  to  measured 
spatial  water  temperatures  every  where  inside  the  torrus 
except  in  the  upper  six  inches,  Figure  9-1.  This  is  not 
surprising  since  the  assumptions  which  were  made  in  the 
derivation  of  equation  (4-18)  are  valid  in  the  core  of  the 
torrus.  The  worst  deviation  was  less  than  -2%  which 
is  well  within  the  experimental  uncertainty  (+5%) . 

The  theoretical  solution  predicts  very  successfully 
the  variations  of  the  spatial  temperatures  inside  the 


Distance  from  top  surface  of  pool  (inches) 
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Measured 

Calculated 


Water  Temperature  (F) 


Figure  9-1.  Comparison  of  calculated  to  measured  water 
temperatures  at  0800  and  1400  hours  on 
12/1/72  - uncovered  pool. 
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torrus . The  cooling  and  heating  periods  of  the  water  are 
shown  clearly  in  Figures  8-23,  8-24.  The  water  tempera- 
ture profile  has  a parabolic  shape  in  early  morning  and 
night.  The  temperature  profile  becomes  distorted  during 
the  day  time  reflecting  the  hotter  temperatures  in  the 
upper  layers  due  to  the  absorption  of  the  solar  radiation. 
The  temperature  profile  flattens  out  right  after  the  sun 
has  disappeared;  it  obtains  its  parabolic  shape  by 
0600  hours.  The  reason  for  the  parabolic  shape  is  that 
colder  water  remains  in  the  bottom  of  the  pool  due  to  its 
density.  Also,  in  the  absence  of  any  incoming  solar  radia- 
tion the  upper  layers  are  cooled  by  evaporation  cooling. 

Maximum  Water  Velocity  Next  to  the  Wall 
Equation  (4-26)  calculated  the  maximum  value  of  the 
velocity  of  the  water  next  to  the  heated  or  cooled  walls 
of  the  pool  (Figure  5-6) . However,  in  the  derivation  of 
(4-26)  a basic  assumption  was  made.  It  was  assumed  that 
the  work  which  is  done  on  the  water  by  the  buoyant  force 
is  all  converted  to  kinetic  energy.  Viscous  dissipation 
was  considered  negligible.  Equation  (4-26)  is  as  follows: 

V = (g0  l AT ) ^ (4-26) 


To  test  the  validity  of  this  assumption,  the  water  velocity 


as  calculated  from  Equation  (4-26)  was  compared  to  the 
velocity  as  calculated  in  ref.  [38].  The  later 
equation  is 


157 


V = 0.766  v (0.952  + -) 


(9-1) 


a 


2 

v 


Equation  (9-1)  calculates  the  f ree-convection  velocity  of  a 
fluid  near  a vertical  heated  wall;  it  was  assumed  that  the 
velocity  has  a parabolic  profile.  The  viscous  dissipation 
effect  was  not  considered  negligible  (Ref.  [38]). 

Equations  (9-1)  and  (4-26)  were  used  to  calculate  the 
water  velocity  with  identical  physical  parameters.  The 
results  are  plotted  in  Figure  9-2.  The  agreement  is 
excellent,  and  therefore,  it  is  concluded  that  the  viscous 
dissipation  can  be  considered  negligible.  The  theoretical 
results  are  accepted  as  more  representative  than  the 
measured  because  of  the  crude  ways  which  were  employed 
to  measure  the  water  velocity. 

It  is  clear  that  f ree-convection  currents  are  present 
throughout  the  24  hours  of  a day.  They  have  zero  velocity 
only  twice;  when  the  heating  of  the  water  is  initiated, 
and  when  the  heating  is  terminated.  It  seems  that  higher 
velocities  are  present  during  the  daytime,  Figure  (9-2). 
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This  is  not  generally  true;  it  all  depends  on  the  relative 
magnitude  of  the  wall  temperature  and  of  the  water  tempera- 
ture during  the  daytime  and  night  time.  The  size  of  the 
f ree-convection  cells  during  the  night  are  more  symmetrical 
than  during  the  day.  This  is  due  to  a more  symmetrical 
cooling  during  the  night  than  heating ' during  the  daytime. 

It  is  also  due  to  the  lack  of  symmetry  in  the  location  of 
the  solarlv  heated  vertical  walls  of  the  pool  during  the 
daytime . 

Temporal  Variations  of  the 
Bulk  Water  Temperature 

Equation  (4-41)  was  examined  in  three  time  frames; 
annual,  monchly,  and  daily.  Application  of  Equation  (4-41) 
in  an  annual  time  frame  revealed  whether  there  was  any 
economic  justification  for  a swimming  pool  owner  to 
reduce  the  thermal  losses  of  the  water,  or  in  addition, 
to  heat  his  pool  with  solar  collectors.  In  a monthly  time 
frame.  Equation  (4—41)  predicted  the  morning  and  afternoon 
water  temperatures  for  several  days  in  a row;  this  proved 
a fine  test  of  Equation  (4-41)  by  comparing  the  theoretical 
results  to  the  measured  results.  Finally,  Equation  (4-41) 
was  examined  in  a daily  time  frame;  temperatures  of  the 
kulk  water  were  calculated  at  the  end  of  each  hour  for  a 


160 


whole  day.  The  daily  time  frame  provided  a good  test  for 
the  validity  of  Equation  (4-41)  since  the  calculated  water 
temperature  must  respond  instantaneously  to  the  continually 
changing  physical  boundary  conditions. 

Annual  Variations 

The  Annual  Variations  of  the  bulk  water  temperature 
have  been  calculated  and  are  shown  in  Figure  5-8.  The 
measured  Annual  Variations  of  the  bulk  water  temperature 
of  an  uncovered  and  of  single  Sears  covered  pool  are  shown 
in  Figure  8—15.  The  measured  temperatures  are  compared 
to  the  calculated  in  Figure  9-3. 

Starting  in  April  of  1972,  the  bulk  water  temperature 
six  cases  rises  to  a maximum  in  July.  This  maximum 
point  does  not  coincide  with  the  maximum  amount  of  solar 
energy  which  is  absorbed  by  the  water  which  occurs 
mid-June  of  1972  in  Figure  5-9.  However,  from  Figure  5-10 
the  convection  losses  are  lower  in  July  than  in  mid-June . 

The  convection  losses  are  lower  in  July  because  the 
ambient  air  temperature  reaches  a maximum  point  in  July,  as 
is  shown  in  Figure  8-15.  From  Figures  5-9  and  5-10,  it 
can  be  seen  that  the  magnitude  of  the  difference  in 
convection  losses  is  much  higher  than  in  the  absorbed 
solar  energy  by  the  water;  therefore,  the  net  effect  is  a 
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Figure  9-3. 


Comparison  of  calculated  with  measured  bulk 
water  temperatures 
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higher  water  temperature  in  July  than  in  mid-June. 

In  August  and  September  of  1972,  the  amount  of  the 
absorbed  solar  radiation  decreases  as  expected,  since 
the  available  solar  energy  flux  decreases.  The  reduced 
amounts  of  absorbed  solar  energy  would  be  expected  to 
cause  a decrease  in  the  bulk  water  temperature . Since  the 
change  of  the  convection  losses  is  not  as  high  in  magnitude 
as  the  change  in  the  absorbed  solar  energy  in  the  period 
of  July,  August,  and  September,  the  net  effect  is  a small 
decrease  in  the  water  temperature.  The  explanation  again 
lies  in  the  small  decrease  in  the  ambient  air  temperature 
during  this  period.  The  air  temperature  does  not  decrease 
as  fast  as  the  available  solar  energy  flux  does  because 
of  the  thermal  capacity  effect  of  the  earth  (the  ambient 
air  is  heated  by  the  earth).  Dickert  [36]  has  measured 
and  calculated  the  soil  temperature  at  various  depths 
at  the  same  site  where  the  measurements  of  the  swimming 
pools  were  taken.  Dickert  compared  the  soil  temperatures 
to  the  measured  ambient  air  temperature  as  a function  of 
time;  this  comparison  appears  in  Appendix  H. 

In  September,  the  amount  of  absorbed  solar  energy 
continues  to  decrease  at  about  the  same  rate  to  a minimum 
point  in  December . The  ambient  air  temperature  decreases 
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at  an  accelerating  rate  in  this  period.  The  earth  can  no 
longer  provide  any  more  heat  to  the  ambient  air.  Therefore 
the  convection  losses  decrease  also  at  an  accelerating  rate 
The  net  result  is  a sudden  decrease  in  the  water  tempera- 
ture of  the  pools. 

At  the  end  of  December,  the  amount  of  absorbed  solar 
energy  starts  increasing.  However,  the  ambient  air 
temperature  does  not  start  increasing  until  February. 

Despite  the  increasing  amounts  of  solar  flux  during 
January,  the  ambient  temperature  continues  to  decrease  due 
to  the  earth's  cool  temperature.  The  earth  is  now  cooling 
the  ambient  air.  The  decreasing  ambient  temperature  and 
the  increasing  amounts  of  solar  flux  result  in  an  increase 
of  convection  losses  in  January.  From  Figures  5-9  and 
5-10  it  can  be  seen  that  the  curves  of  the  absorbed  solar 
energy  and  of  the  convection  losses  follow  the  same  trend; 
therefore,  it  is  expected  that  the  water  temperature  would 
reach  its  minimum  in  mid-December.  From  Figure  8-15, 
the  minimum  was  measured  in  January;  the  difference  can 
probably  be  attributed  to  the  onset  of  increasing  convection 
losses  in  January. 

There  is  an  excellent- agreement  between  calculated 
and  measured  water  temperatures  for  the  uncovered  pool 
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until  August  (Figure  9-3); in  August  the  deviation  is  -2.87%. 
After  August,  the  calculated  values  are  still  lower  than 
the  measured  but  the  deviation  is  smaller,  until  December. 
The  deviation  in  December  is  -3.7%.  From  December  to  April, 
there  is  good  agreement  between  calculated  and  measured 
water  temperatures.  The  lower  calculated  values  can  be 
attributed  to  the  experimental  limitations  of  the  multitude 
of  experimental  data  and  many  assumptions  which  are  re- 
quired for  the  calculations.  The  measured  apparent  water 
absorptivity  could  be  as  much  as  ±10  percent  off,  especially 
at  the  low  altitude  angles  of  the  sun. 

For  the  covered  pool,  the  deviation  between  calculated 
and  measured  values  reaches  a maximum  value  in  July  (+4.12%) 
in  December  (-4.2%),  and  in  February  (-3.1%).  Again,  these 
deviations  are  due  mainly  to  experimental  limitations 
in  measuring  the  apparent  transmissivities  of  the  plastic 
covers . It  was  obvious  to  the  author  from  the  beginning 
that  the  measured  plastic  transmissivity  in  the  "Black  Box" 
is  higher  than  that  in  the  pool.  This  is  due  to  the 
higher  reflection  losses  which  are  caused  by  the  presence 
of  the  water  under  the  plastic  cover;  in  the  "Black  Box" 
there  was  no  practical  way  to  make  transmissivity  measure- 
ments with  water  under  the  plastic  cover.  The  effect 
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of  the  higher  reflection  losses  due  to  the  water  was 
estimated  to  be  15%,  and  therefore,  all  the  apparent 
transmissivity  factors  were  lowered  by  15%.  It  was  obvi- 
ously assumed  that  this  15%  reduction  in  apparent  trans- 
missivities holds  for  all  the  months  in  the  year;  this 
assumption  is  probably  most  responsible  for  the  disagree- 
ments between  calculated  and  experimental  water  tempera- 
tures as  were  stated  in  the  beginning  of  this  para- 
graph . 

The  Annual  Variation  of  the  Temperature  Gains  reveals 
another  view  of  the  heat  transfer  mechanisms  in  a swimming 
pool.  Figure  5-14  shows  the  variation  of  the  temperature 
gains.  A temperature  gain  is  the  total  water  temperature 
increase  of  the  water  temperature;  it  is  the  difference 
between  the  morning  and  early  evening  water  temperature . 
Figure  5-14  shows  a maximum  value  for  the  temperature  gain 
in  May  and  a minimum  in  September . The  maximum  value 
is  due  to  the  larger  amounts  of  solar  energy  absorbed 
by  the  water  and  the  longer  daylight  hours.  During  the 
daylight  hours,  the  convection  losses  are  greater,  for 
example,  in  July  than  May,  because  there  are  more 
daylight  hours  in  July  than  May.  The  opposite  holds  true 
in  September.  Here  the  absorbed  solar  energy  starts 


166 


ing  at  an  accelerating  rats,  and  the  daylight  hours 
become  less  than  the  nighttime. 

The  effect  of  covering  the  water  with  a plastic 
cover  on  reducing  the  heat  losses  of  the  water  can  be 
seen  in  Figures  8-15  and  5-8.  The  presence  of  the  plastic 
cover  reduces  the  amount  of  the  absorbed  solar  energy 
considerably;  the  cover  also  impedes  the  transfer  of  energy 
away  from  the  top  surface  of  the  water  as  convection  losses . 
But  the  greatest  impact  of  the  plastic  cover  is  in  the 
elimination  of  the  evaporation  losses.  The  elimination 
of  the  evaporation  losses  greatly  compensates  for  the 
reduction  of  absorbed  solar  energy,  thus  resulting  in  a 
higher  water  temperature  for  the  covered  pools  . 

Monthly  Variations 

The  morning  and  afternoon  bulk  water  temperatures 
were  calculated  and  measured  for  five  consecutive  days. 

This  time  period  was  chosen  because  the  ambient  air 
temperature  dropped  abruptly  due  to  the  passage  of  a mass 
of  cold  air.  Figure  9-4  compares  the  calculated  tempera- 
tures to  the  measured. 

Both  the  covered  and  the  uncovered  swimming  pools 
had  approximately  the  same  bulk  water  temperature  on 
H/13/72,  which  is  the  beginning  of  the  five  day  period  of 


1972 

TIME  (Month/Day) 

Figure  9-4.  Comparison  of  calculated  with 

measured  bulk  water  temperatures. 
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interest.  It  is  possible  for  both  pools  to  have  approxi- 
mately the  same  bulk  temperature  whenever  the  net  amount 
of  energy  transferred  to  both  pools  is  about  the  same. 

As  can  be  seen  from  Table  5-1  where  the  average  daily  air 
temperatures  are  listed,  the  convection  losses  in  this 
time  period  should  be  minimal  because  the  difference  between 
the  bulk  water  temperature  and  the  air  temperature  is 
minimal  for  both  pools.  This  is  true  up  to  11/15/72  when 
a mass  of  cold  air  swept  the  test  site.  Also  from  Table 
5-1,  it  can  be  seen  that  the  air  temperature  was  about  72°F 
from  10/23/72  to  11/14/72.  It  is  clear,  then,  that  both 
pools  had  minimal  convection  losses  for  a relatively  long 
period  of  time.  The  greater  amounts  of  solar  radiation 
that  the  uncovered  pool  absorbed  were  offset  by  thermal 
losses  due  to  evaporation.  The  result  was  an  approximatelv 
equal  amount  of  net  energy  transferred  to  both  pools. 

The  effect  of  the  plastic  cover  became  very  evident 
after  11/15/72  when  the  mass  of  the  cold  air  swept  through 
the  test  site.  Usually,  after  the  passage  of  a cold  front, 
the  ambient  air  is  cold  and  dry,  and  the  amount  of  solar 
radiation  reaches  a level  which  is  higher  than  the  monthly 
average.  These  observations  are  substantiated  by  the 
measurements  of  ambient  air  and  solar  radiation  for  the 
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period  of  11/13  through  11/18  (Tables  5-1  and  5-2, 
respectively) . As  a result,  the  evaporation  losses  would 
increase  for  the  uncovered  pool  due  to  the  drier  ambient 
air.  The  convection  losses  would  increase  for  both  pools 
due  to  the  lower  ambient  air  temperature.  The  increased 
amounts  of  absorbed  solar  radiation  cannot  offset  the 
increased  amounts  of  convection  and  evaporation  losses. 

The  bulk  water  temperature  would  drop  for  both  pools. 
However,  it  drops  more  in  the  uncovered  pool  because  of  the 
increased  evaporation  losses. 

Specifically,  the  bulk  water  temperature  in  the 
covered  pool  increased  every  day  during  the  daytime  hours 
except  on  the  day  of  11/17.  On  that  day,  the  solar 
radiation  was  below  the  monthly  average  for  the  month  of 
November;  the  thermal  losses  were  greater  than  the  thermal 
gains . During  night  time , the  bulk  water  temperature 
decreased,  as  expected.  In  the  uncovered  pool,  the  bulk 
water  temperature  during  the  daytime  hours  decreased  except 
on  11/13,  11/16,  and  11/18.  On  these  days,  the  thermal 
gains  from  the  absorbed  solar  radiation  were  greater  than 
the  thermal  losses  from  convection  and  evaporation. 

The  calculated  temperature  profiles  follow  quite 
closely  the  measured  (Figure  9-4) . 


In  the  uncovered  pool. 


170 


the  largest  deviation  between  measured  and  calculated  was 
-2.61s;  in  the  covered  pool  the  largest  deviation  was 
+4.92%.  This  excellent  agreement  between  calculated  and 
measured  temperatures  proyes  the  validity  of  Equation 
(4-41)  . 

Dally  Variations 

From  a practical  standpoint,  the  measured  bulk  water 
temperatures  on  an  annual  and  monthly  time  frame  were  of  a 
pool  covered  with  a single  Sears  plastic  cover.  In  a 
long  time  frame,  it  was  not  practically  possible  to  have 
several  pools  to  test  the  several  kinds  of  plastic  covers. 
However,  in  a short  time  frame,  such  as  daily,  it  was 

practical  to  cover  a pool  with  a different  plastic  cover 
for  each  day  of  testing. 

During  December  of  1972,  the  following  plastic  covers 
were  tested:  one  layer  of  Sears  cover;  two  layers  of  Sears 

cover;  one  layer  of  "Bubble”  plastic;  two  layers  of  "Bubble" 
plastic.  Bulk  water  temperatures  were  measured  at  the 
end  of  each  hour.  The  calculated  results  are  compared  to 
the  measured  in  Figures  9-5  through  9-9.  Figure  9-5 
shows  the  temperature  profile  of  a pool  with  no  cover. 

The  greatest  temperature  increase  in  the  water  bulk 
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12/11/72  - Double  Sears  covered  pool. 
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temperature  was  measured  in  the  pool  v/ith  the  single  Sears 
cover  and  also  in  the  uncovered  pool.  However,  the  uncov- 
ered pool  lost  a lot  of  its  thermal  energy  at  the  end  of 
the  day,  and  its  afternoon  temperature  returned  to  its 
morning  level.  Therefore,  it  is  concluded  that  the  most 
beneficial  way  is  to  cover  the  water  of  a pool  with  a 
single  layer  of  Sears  plastic.  The  poor  performance  of 
the  two  layers  of  Sears  plastic  can  be  attributed  primarily 
to  the  great  reduction  of  the  absorbed  solar  radiation 
by  the  water  which  is  due  to  the  low  apparent  absorptivity 
of  the  two  layers  of  the  Sears  plastic.  The  poor  perform- 
ance of  the  "Bubble " plastic  can  be  attributed  to  the 
lack  of  entrapped  air  in  the  regions  of  the  plastic  around 
the  air  bubbles  (Figure  5-7)  . 

Figures  9-5  through  9-9,  other  than  providing  compara- 
tive information  of  the  various  plastic  covers,  present 
also  a picture  of  the  heating  and  cooling  periods  of  the 
water.  During  those  December  days,  there  was  a cooling 
period  of  the  water  from  1600  to  1000  hours,  and  a heating 
period  j.rom  1000  to  1600  hours.  It  is  true  that  the  sun 
rose  around  0800  solar  time.  However,  its  heating  effect 
on  the  water  was  not  measured  until  1000  hours.  This 
time  delay  in  the  morning  sun's  heating  effect  on  the  water 
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is  similar  to  a time  delay  in  the  rise  of  the  ambient  air 
temperature.  However,  this  mass  capacity  effect  is  not 
the  only  reason  for  the  time  delay.  As  was  discussed 
ear^er'  the  apparent  absorptivity  of  the  water  and  the 
apparent  transmissivity  of  the  plastic  covers  all  have 
minimum  values  at  large  sun  incidence  angles  (early  in  the 
morning  and  late  in  the  afternoon) . 

After  1000  hours,  there  is  a steady  increase  in  the 
bulk  water  temperature  until  about  1500  hours.  Again,  it 
is  true  that  the  maximum  amount  of  solar  radiation  occurred 
around  1200  hours  solar  time.  This  time  delay  again 
in  measuring  the  effects  of  the  decreasing  solar  radiation 
on  the  bulk  water  temperature  is  due  to  thermal  capacity 
effects.  The  walls  of  the  swimming  pool  which  had  been 
heated  by  the  sun  do  not  lose  their  thermal  energy  instan- 
taneously with  the  decreasing  amounts  of  solar  radiation. 

It  takes  a certain  time  for  their  cooling.  In  addition, 
the  values  of  the  apparent  absorptivity  and  transmissivi- 
ties are  high  during  this  period  of  the  day;  the  result 
is  a more  efficient  absorption  of  the  solar  radiation, 
even  when  the  solar  radiation  starts  decreasing. 

After  1500  hours,  there  is  generally  a temperature 
drop  in  the  bulk  water  temperature.  This  is  the  result 
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of  decreased  amounts  of  solar  radiation.  This  temperature 
drop  is  especially  evidenced  in  the  uncovered  pool  (Figure 
9-5) . The  bulk  water  temperature  drops  very  slowly  in  the 
covered  pools  (Figures  9-6  through  9-9).  The  explanation 
lies  in  the  fact  that  a plastic  cover  impedes  the  transfer 
of  energy  and  mass  through  it.  The  plastic  cover  reduces 
the  amount  of  solar  radiation  that  can  be  absorbed  by  the 
water.  But  it  also  reduces  the  amount  of  energy  that  can 
be  lost  by  convection  and  evaporation  from  the  top  surface 
of  the  pool.  That  explains  the  flatter  temperature  profiles 
m the  afternoons  for  the  covered  pools  as  compared  to  the 
uncovered  pool. 

Another  interesting  fact  is  that  the  water  temperature 
of  the  uncovered  pool  returns  to  about  the  same  temperature 
level  in  the  evening  as  compared  to  the  morning  level. 

This  implies  that  there  was  not  a net  increase  or  decrease 
in  the  water  temperature  of  the  uncovered  during  the  fifteen 
hour  time  increment  of  that  particular  day.  The  same  does 
not  occur  in  all  of  the  covered,  unheated  pools.  The 
evening  temperature  level  of  the  water  is  always  higher 
than  the  morning's,  which  implies  a net  temperature  increase 
during  the  daylight  hours.'  This  increase  is  the  result 
of  the  plastic  cover  impeding  energy  and  mass  transfers. 
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as  was  explained  in  the  previous  paragraph. 

There  is  an  excellent  agreement  between  calculated 
and  measured  temperatures  as  demonstrated  by  Figures  9-5 
through  9-9.  in  general,  the  calculations  overpredict 
the  bulk  temperatures  by  a maximum  +2.00%.  There  is  a 
similarity  in  the  temperature  profiles  in  Figures  9-5 
and  9-8;  the  calculated  profiles  cross  the  measured  at 
about  1200  hours.  No  conclusions  can  be  truly  drawn  from 
these  crossings  because  of  the  less  than  1%  deviation 

between  the  calculated  and  measured  temperatures  at  1200 
hours . 

Swimming  Pool  Heated  by  Solar  Collectors 
and  Covered  by  Single  Sears  Cover 

Special  attention  should  be  given  to  the  pool  which 
is  heated  by  solar  collectors  and  is  covered  by  a single 
Sears  cover.  As  it  was  shown  earlier,  this  is  the  only 
pool  which  has  its  bulk  water  temperature  above  80 °F 
year  around  in  Gainesville,  Florida.  Water  temperature 
measurements  of  such  a pool  were  made  on  12/13/72;  the 
measured  bulk  temperatures,  along  with  the  calculated 
are  shown  in  Figure  9-10. 

The  temperature  profile  is  much  steeper  than  in  any 
of  the  other  covered  and  uncovered  pools.  The  bulk  water 
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TIME  (Hours) 

Figure  9-10.  Comparison  of  calculated  and  measured 
bulk  water  temperatures  on  12/13/72  - 
Single  Sears  covered  and  solarly  heated 
pool 
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temperature  depends  very  much  on  the  amount  of  solar  radia- 
tion. That  explains  the  very  sudden  leveling  at  1500 
hours.  This  sudden  leveling  reminds  one  of  the  temperature 
drop  at  1500  hours  in  the  uncovered  pool  for  the  same 
reason.  Unlike  the  uncovered  pool,  though,  the  bulk 
temperature  remains  almost  constant  after  1500  hours  simi- 

i 

lar  to  the  other  covered  pools.  The  plastic  cover  helps 
keep  the  water  temperature  almost  constant  for  a couple  of 
hours. 

In  the  calculation  of  the  bulk  water  temperature  for 
the  pool  which  was  covered  with  a single  Sears  cover  and 
was  also  heated  by  ten  flat  plate  collectors,  information 
was  needed  as  to  the  efficiency  of  the  collectors  and 
its  variation  with  incidence  angle.  Since  the  collector 
efficiency  was  not  readily  available  for  the  collectors 
which  were  used,  it  was  decided  that  the  next  best  thing 
to  having  the  efficiency  curve  of  the  actual  collectors 
was  to  use  the  efficiency  curve  of  a similar  collector 
design.  The  collector  which  was  found  similar  in  its 
design  to  the  collectors  used  in  the  measurements  is 
made  by  NASA/Honeywell;  it  is  coated  with  a non-selective 
black  coating  with  e = a = 0.97,  and  has  one  glass  cover. 

Ref.  [39]  experimentally  found  the  efficiency  curve 


for  the  NASA/Honeywell  collector.  The  efficiency  as  a 
function  of  incidence  angle  is  equal  to 
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= 0.85K  - (1.1390  + O.16102) 


(9-2) 


where 


K = 1.0  - 0.078  — 

(p 


l.Oj 


(9-3) 


and 


(9-4) 


Equation  9 ^ was  utilized  to  predict  the  water  temperature 
in  Figure  9-10.  There  is  excellent  agreement  of  the 
calculated  temperatures  to  the  measured  up  to  1700  hours. 
At  1700  and  1800  hours,  the  deviation  is  about  +1.7%. 

It  is  probably  due  to  the  large  uncertainty  of  the  value 
of  the  apparent  transmissivity  of  the  plastic  cover 
at  large  sun  incidence  angles. 

Chillier  [14]  had  concluded  that  to  raise  the  water 
temperature  10°F,  one  needs  flat  plate  collectors  of  an 
area  equal  to  the  surface  of  the  swimming  pool.  The 
present  work  has  found  that  more  area  of  flat  plate 
collectors  is  needed.  As  Figure  9-10  shows,  to  raise 
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the  water  temperature  10°F  in  Gainesville,  Florida, 
one  needs  flat  plate  collectors  of  an  area  equal  to  1.5 
times  the  surface  of  the  pool  for  the  winter  months. 


Comparison  of  Benefits  and  Costs  of 
Heating  a Swimming  Pool 

As  was  discussed  in  the  previous  paragraph,  one  very 
good  way  to  supplement  the  heating  effect  of  the. sun  is 
by  covering  the  water  with  a single  plastic  cover.  its 
cost  is  very  low  (about  $10  to  cover  a 15  foot  diameter 
pool).  The  benefits  of  using  a plastic  cover  can  be 
seen  clearly  in  Figure  9-3.  if  it  is  assumed  that  a water 
temperature  of  80 °F  provides  comfortable  swimming,  then 
the  following  observations  can  be  made.  The  swimming  pool 
with  a plastic  cover  in  Gainesville,  Florida,  permitted 
swimming  in  the  first  week  of  March,  whereas  the  uncovered 
pool  did  not  permit  swimming  until  mid-May.  Also,  the 
covered  pool  extended  the  swimming  season  almost  to  the 
middle  of  October;  with  the  uncovered  pool,  the  swimming 
season  ended  in  the  last  week  of  September.  m conclusion, 
the  swimming  season  in  Gainesville,  Florida,  was  extended 


by  almost  three  months  simply  by  covering  the  water 
with  a single  layer  of  plastic.  Obviously,  this  is  the 
least  expensive  way  to  raise  the  water  temperature  of  the 
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pool . 

Some  methods  of  heating  the  water  are  with  solar 
collectors,  electric  heating,  and  natural  gas  heating. 

From  Figure  5-8,  the  water  temperature  of  an  uncovered  pool 
is  below  80 °F  for  seven  months  (October  through  April) . 

To  raise  the  water  temperature  to  80 °F  during  these  seven 
months,  it  is  necessary  to  supply  8.914x10  BTU  to  the 
water  of  a 15  foot  diameter  pool.  More  energy  is  actually 
needed  to  balance  the  thermal  losses  of  the  water  to  its 
surroundings.  However,  for  comparison  purposes  the  calcu- 
lated  energy  of  8.914x10  BTU  is  used.  For  electric  heat- 
ing with  an  efficiency  of  90  percent  and  the  cost  of  4.22 
cents/kw-hr,  the  total  fuel  cost  is  $12 25/year.  For  natural 
gas  heating  with  an  efficiency  of  60%  and  the  cost  of 
64  cents/10^  BTU,  the  total  fuel  cost  is  $95/year.  It 
should  be  noted  that  these  costs  do  not  include  the  purchase 
of  the  electric  heater  or  the  gas  furnace. 

If  the  swimming  pool  is  heated  during  these  seven 
months  with  solar  collectors  and  these  solar  collectors 
are  also  used  for  summer  air  conditioning,  then  the 
total  materials  cost  for  a 270  sq.  ft.  area  of  collectors 
is  equal  to  270x10x7/12  = $1575.  Due  to  this  high  initial 
cost,  the  return  on  investing  the  $1575  would  have  to  be 
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added  to  the  materials'  cost  of  the  collectors.  Assuming 
that  the  amount  of  $1575  is  invested  at  a compounded 
interest  rate  of  5.75%  annually,  paid  twice  a year,  then 
the  value  of  the  investment  after  15  years  (assumed  life 
of  the  collectors)  is 


1575  fl  + ---°--75  j 
l.  2 J 


$3686 . 


Therefore,  the  annual  cost  of  heating  the  water  with  solar 
collectors  is 


-J5-  = $ 246/year . 

As  is  shown  in  Figure  5-8  an  area  of  270  ft2  of  solar 
collectors  can  keep  the  water  temperature  of  a 15  foot 
diameter  pool  above  80 °F  year  around.  The  cost  of  the  solar 
collectors  was  assumed  to  be  $10/square  foot.  The  $246/ 
year  cost  for  the  flat  plate  collectors  is  reduced  consid- 
erably by  using  the  same  collectors  for  winter  space  heat- 
ing. The  collectors  can  be  used  for  swimming  pool  heating 
when  the  temperature  of  the  water  in  the  collectors  is 
below  the  temperature  of  the  thermal  energy  storage  tank. 

From  these  considerations,  it  is  concluded  that 
heating  the  water  with  solar  collectors  is  not  the  most 
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economical  way.  Electric  heating  is  by  far  the  most 
expensive.  At  today's  prices  heating  the  water  with  natural 
gas  is  quite  economical.  However,  the  market  availability 
of  natural  gas  in  the  future  is  not  certain  at  all.  Already 
some  states  in  the  United  States  have  passed,  or  are 
considering,  laws  prohibiting  natural  gas  heating  for 
swimming  pools  due  to  its  scarcity.  It  is  the  author's 
belief  that  utilization  of  coal  or  oil  for  swimming  pool 
heating  could  also  become  prohibited  in  the  near  future 
by  legislation  or  by  economic  considerations. 


CHAPTER  X 


CONCLUSIONS 

1.  The  analytical  and  experimental  results  of  this 
investigation  confirm  the  presence  in  a pool  of  a flow 
regime  in  which  (a)  there  is  zero  vorticity,  (b)  there  is 

a vertical  temperature  gradient,  and  (c)  the  water  tempera- 
ture is  almost  constant  in  the  horizontal  direction. 

2.  The  spacial  water  temperature  variations  inside 
the  central  part  of  the  above  described  flow  regime  can 
be  calculated  by  solving  equation  (4-18)  by  Finite  Cosine 
Fourier  Transforms. 

3.  The  pattern  of  the  free  convection-induced  circu- 
latory currents  for  a circular  swimming  pool  has  the  shape 
of  a torroid  (Figure  4-4) . 

4.  The  effect  of  viscous  dissipation,  as  compared  to 

the  buoyancy  forces  on  the  water  motion,  can  be  considered 
small. 

5.  The  method  of  calculating  all  energy  and  mass 

transfers  from  the  water  of  the  swimming  pool  has  been 
proven  to  be  successful  in  predicting  the  bulk  temporal 
variation  of  the  water  temperature.  The  bulk  temporal 
variations  can  be  predicted  accurately  in  three  time 
frames:  annual,  monthly,  and  daily. 
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6.  The  most  effective,  and  relatively  economical  de- 
sign to  raise  the  water  temperature  is  to  cover  the  sur- 
face of  the  pool  with  a plastic  cover  and  pump  the  water 
through  flat  plate  collectors . 

7 . The  amount  of  area  of  flat  plate  collectors  which 
is  required  in  Gainesville,  Florida  to  raise  the  water 
temperature  10 °F  is  equal  to  1.5  times  the  surface  area 

of  the  pool  for  the  winter  months. 

8.  The  least  expensive  design  to  supplement  the 
heating  effect  of  the  sun  is  to  cover  the  swimming  pool 
with  a single  plastic  cover  which  floats  on  the  water 
surface.  The  water  temperature  measurements  which  were 
made  in  Gainesville,  Florida,  for  one  year  have  shown  that 
a single  plasric  covered  pool  extended  the  swimming  season 
from  middle  of  September  to  middle  of  October,  and 
permitted  swimming  earlier  from  March  to  May.  This 
extension  of  the  swimming  season  was  based  on  assuming 
that  a water  temperature  of  80 °F  provides  comfortable 
swimming . 

9.  A single  layer  of  polyethylene  plastic  is 
recommended  for  covering  the  surface  of  the  swimming 
pool  because  of  its  good  thermal  performance,  market 
availability,  and  very  low  cost. 
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10.  The  Bubble  plastic  is  not  recommended  for  covering 
the  surface  of  the  swimming  pool  because  of  its  very  poor 
ability  to  transmit  solar  radiation. 


APPENDIX  A 


SOLUTION  OF  EQUATION  (4-17) 

BY  FINITE  COSINE  FOURIER  TRANSFORMS 


Equation  (4-17)  is 
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According  to  ref.  [40]  , if  the  Finite  Cosine  Fourier  Transform 

of  the  water  temperature  T(z,t)  is 
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then  the  general  solution  is 
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Also,  according  to  Ref-  [40] 


we  have 
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Transforming  Equation  (4-17)  by  utilizing  equations  (A-l),  (A-3) 

and  (A-4)  we  have 
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Rearranging  (A-5)  we  have 
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Then  equation  (A-6)  can  be  rearranged  as  follows: 
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Equation  (A  8)  is  an  inhomogeneous  first  order  ordinary  differential 
equation;  the  solution  is 
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where  is  a constant  to  be  evaluated  from  boundary  condition  1. 
Performing  the  integration  in  Equation  (A-9)  we  have  that 
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(A-10) 


To  evaluate  the  constant  G;l  in  equation  (A-10)  , boundary  condition 
1 has  to  be  transformed  by  Finite  Cosine  Fourier  Transforms.  Boundary 
condition  1 stated  that 


at  t = 0 + T (z , 0)  = E z +Ez+E. 

12  3 


By  transforming  boundary  condition  1 we  have 

l , , 

F [T  (z  ,0)  ] = T (w,0 ) = E / z2cos 

1 Q 1 J 


dz  + 
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+ E 
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dz  + E / cos  p^l  dz 
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(A-ll) 


Integrating  (A-ll)  we  have  that 


T(w,0)  = 
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(A-12) 


Equation  (A-12)  is  then  the  transformed  boundary  condition  1.  We  can 
now  apply  Equation  (A-12)  into  equation  (A-10)  to  evaluate  the  constant 


, and  we  have  that 
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To  complete  the  solution  T(0,0)  needs  to  be  evaluated. 


I l 

T(0,0)=  / T(0,0)dz  = /(E^z  +E^z+E^)dz  = E^  — + E^  + E £ (A-14) 


Substituting  equation  (A-14)  into  (A-2)  we  have 


T(z,t) 


E 

_1 
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00 

£ T (w,  t)  cos 

W=1 


f "V 
W7TZ 

U J 


(A-15) 


v/here  T(w,t)  is  defined  by  equations  (A-10)  and  (A-13). 


APPENDIX  B 


FORTRAIN  PROGRAM  OF  EQUATION  (4-18) 


100  Cl=43. 45833 

102  T=0. 

103  P=0.0 

104  TBAR=0 . 0 

105  C2=-25. 64077 
110  03=73.85999 
115  C4=-44. 05532 
120  C5=-35. 03184 
125  06=25.86181 
130  C7=-16 .8228 
135  08=4 .58333 
140  C9=9 .526278 
145  010=1.347427 
150  011=5.132332 
155  C12=-8 . 5 

160  013=3.25 
165  C14=-. 751386 
170  C15=-8. 105816 
180  016=4.416668 
185  017=1.154698 
190  El=-3 .3000126E-01 
195  E2=l. 3160042 
200  E3=59. 105831 
210  Dl=-12. 11464 
220  D2=-23. 07272 
230  D3=61. 635036 
240  D4=-37. 7271176 
250  D5=-32 .319607 
260  D6=25 .2235148 
270  D7=-15. 106789 
280  D8=3. 85744 
290  D9=10. 86322 

300  D10=-. 23689148 

301  Dll=3. 126473 
310  D12=-5. 6588413 
315  D13=2. 5559101 
320  D14=-. 55692188 
325  D15=-6 .003132 
330  Pl=3. 141592654 

340  A=5 . 46292E-03 

341  Y=. 39729 
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345  0K= . 34 
350  OL=.3 . 2.5 
355  B= .76 
360  TI=24  . 

365  Z=0.0 

366  ON=200 . 

370  370  P=0.0 

371  IF ( Z . GT.  OL)  GO  TO  825 
375  C18=2 . *PI/24 . 

380C  FOLLOWING  EQUATION  IS  FOR  TBAR 
381  W=l. 

390  30  A1=A* ( (W*PI/OL) **2) 

391  BBB=A1*TI 

392  IF  ( BBB  .LE.  88.)  GO  TO  31 

393  BB=0 . 0 

394  GO  TO  32 

395  31  BB=1./EXP (BBB) 

396  32  AA=1 ./EXP (Y*OL) 

39  7 IF  (W  .GE.  ON)  TO  TO  20 

398  H=- ( Y * * 2 * ( 1 . -B ) *A/OK) * ( ( (-1 . ) *AA) -1 . ) / (Al/A+Y**2 ) 

405  Rl-~  (H*C1-  (Dl*A/OK)  ) /A1+  ( (H*C2-  (D2*A/OK)  )/(C18**2 
410S+A1**2)  ) * (Al*COS  (C18*TI) +C18*SIN  (C18*TI)  ) + ( (H*C3- (D3*A/OK 
415&/ (C18**2+A1**2) ) * ( A1*SIN (C18*TI) -C18*COS (C18*TI) ) 

420  R2= ( (H*C4- (D4*A/OK) ) / ( (2 . *C18) **2+Al**2) ) * (Al* 

4 25 SCOS (2.*C18*TI)+2.*C18*SIN(2. *C18*TI) ) + 

430& ( (H*C5- (D5*A/OK) ) / ( (2 . *C18) **2+Al**2) ) * (Al* 

435&SIN (2 . *C18  *TI ) -2 . *C1 8 *COS (2 . *C18*TI) ) 

440  R3= ( (H*C6- (D6*A/OK) )/ ( (3 . *C18) **2+Al**2) ) * (Al* 

445SCOS (3 . *C18*TI) +3 . *C18*SIN (3 . *C18*TI) ) + 

450s ( (H*C7- (D7*A/OK) )/( (3.*C18) **2+Al**2) ) *(A1* 

455  SSIN ( 3 . *C18  *TI ) -3.*C18*COS (3. *C18*TI) ) 

460  R4= ( (H*C8- ( D8*A/OK) )/( (4.*C18) **2+Al**2) ) * (Al* 

465SCOS (4 . *C18*TI) +4 . *C18*SIN (4 . *C18*TI) ) + 

470s ( (H*C9- (D9*A/OK) ) / ( ( 4 . *C18) **2+Al**2) ) * (Al* 

475SSIN (4 . *C18*TI) -4 . *C18*COS (4.*C18*TI) ) 

4 80  R5=((H*C10-( D10 * A/ OK ) ) / ( ( 5 . *C18) **2+Al**2) ) * (Al* 

485SCOS (5 . *C18  *Tl) +5 . *C18*SIN (5 . *C18*TI) ) + 

490S ( (H*C11- (Dll*A/OK) ) / ( (5 . *C18 ) **2+Al**2) ) * (Al* 

495 SSIN (5 . *C18*TI) -5 . *C18*COS (5 . *C18*TI) ) 

500  Ro=  ( (H*C12--  (D12  *A/OK)  ) / ( (6  . *C18)  **2+Al**2)  ) * (Al* 

505  SCOS  (6  . *C18*TI)  +6  . *C18*SIN  (6  . *C18*TI)  ) + 

510s ( (H*C1 3- ( Dl 3*A/OK ) )/( (6 . *C18) **2+Al**2) ) * (Al* 

515 SSIN (6 . *C18  *TI ) -6 . *C18*COS ( 6 . *C18*TI) ) 

520  R7=( (H*C14- (Dl4*A/OK) )/ ( (7. *C18) **2+Al**2) ) * (Al* 

5 25 SCOS (7. *C18*TI) +7 .*C18*SIN (7. *C18*TI) )+ 

5 30S  ( (H*C15- ( Dl5*A/OK) )/ ( ( 7. *C18) **2+Al**2) ) * (Al* 
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535&SIN (7 . *C18*TI ) -7 . *C18*COS ( 7 . *C18*TI ) ) 

540  R8=.(H*C16/(  (8  .*C18)  **2+Al**2)  ) * (Al* 

545&COS  (8 . *C18*TI) +8  . *C18*SIN ( 8 . *C18*TI) ) + 

550&(H*C17/ ( ( 8 . *C18) **2+Al**2) ) *(A1* 

555&SIN (8 . *C18*TI) -8 . *C18*COS (8 . *C18*TI) ) 

556  G=- (2 . *0L*E1) / (Al/A) - (2 . *E2) / (Al/A) - (H*C1- (Dl*A/OK) ) /Al- 
55 7&  (H*C2- (D2*A/OK) ) *A1/ ( (C18**2) +A1**2) +C18* (H*C3- 
55 8 & (A* D 3/OK ) ) / { (C18**2) +A1**2 ) - (H*C4- (A*D4/OK) ) *A1/ ( (2 . * 
559&C18) **2+Al**2 ) +2 . *C18* (H*C5- (A*D5/OK) ) /( (2 .*C18) **2  + 
560&A1**2) -Al* (H*C6-(A*D6/OK) )/( (3.*C18) **2+Al**2 ) +3 . *C18* 

56 1& (H*C7- (A*D7/OK) )/ ( (3. *C18) **2+Al**2) -Al* (H*C8- (A*D8/OK) ) 
562&/ ( ( 4 . *C18) **2+Al**2) +4.*C18*(H*C9- (A*D9 /OK) )/((4.*Cl8)** 
563&2+Al**2) -Al* (H*C10- (A*D10/OK) )/( (5.*C18) **2+Al**2) +5 . *C18 
564&* (H*C11- (A*Dll/OK) )/ ( (5 . *C18) **2+Al**2) -Al* (H*C12- (A*D12 
565&/OK)  )/(  (6  . *C18)  **2+Al**2)  +6  . *C18*  (H*C13- (A*Dl'3/OK)  )/(  (6  . 
566  &*C18 ) **2+Al**2) -Al* (H*C14- (A*D14/OK) ) / ( (7. *C18) **2+Al**2) 
567&+7 . *C18* (H*C15- (A*D15/OK) ) /( (7 . *C18) **2+Al**2) -A1*H*C16 / 
568& ( (8 .*C18) **2+Al**2) +8.*C18*H*C17/( ( 8 . *C18) **2+Al**2) 

569  TBAR=G*BB+ 

570&R1+R2+R3+R4+R5+R6+R7+R8 
575  P=TBAR*COS (W*PI*Z/OL) +P 
580  IF ( W .GT.  ON)  GO  TO  20 
585  W=W+2 . 

590  GO  TO  30 
595  20  W=2 . 

600  40  A1=A* ( (W*PI/OL) **2) 

601  BBB=A1*TI 

602  IF ( BBB  .LE.  88.)  GO  TO  61 

603  BB=0 . 0 

604  GO  TO  62 

605  61  BB=1./EXP  (BBB) 

606  62  AA=1 ./EXP (Y*OL) 

610  H=- ( Y**2* (1 .-B) *A/OK) * ( (AA) -1 . ) / (Al/A+Y**2) 

615  Rl= (H*C1- (Dl*A/OK) ) /A1+ ( (H*C2- (D2*A/OK) ) /(C18**2 
620&+Al**2) ) * (Al*COS (C18*TI ) +C18 *SIN (C18*TI) ) +( (H*C3- (D3*A/OK) ) 
625&/(C18**2+Al**2) ) * (A1*SIN (C18*TI) -C18*COS (C18*TI) ) 

630  R2= ( (H*C4-(D4*A/OK) )/( (2.*C18) **2+Al**2) ) * (Al* 

6 35&COS  (2  . *C.18*TI ) +2  . *C18*SIN  (2  . *C18*TI)  ) + 

640&  ( (H*C5-(D5*A/OK) )/( (2.*C18) **2+Al**2) ) * (Al* 

645&SIN (2 . *C18*TI) -2 . *C18*COS (2 .*C18*TI) ) 

650  R3=((H*C6-(D6*A/OK))/( (3.*C18) **2+Al**2) ) * (Al* 

655&COS (3.*C18*TI) +3. *C18*SIN (3 .*C18*TI) ) + 

656 & ( (H*C7-  (D7*A/OK) )/( (3.*C18) **2+Al**2) ) * (Al* 

660 &SIN ( 3 . *C18*TI ) -3 . *C18*COS (3.*C18*TI) ) 

665  R4= ( ( H * C 8 - ( D 8 * A/ OK ) ) / ( ( 4 . *C18) **2+Al**2) ) * (Al* 

670&COS (4 .*C18*TI) +4 . *C18*SIN (4 . *C18*TI) ) + 
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675&  ( (H*C9-  (D9*A/OK)  )/(  (4  .*C18)  **2+Al**2)  ) * (Al* 

680&SIN (4 . *C18*TI) -4 . *C18*COS (4 . *C18*TI) ) 

685  R5=( (H*C10- ( DIO *A/OK) ) / ( ( 5 . *C18 ) **2+Al* *2 ) ) * (Al*' 

690&COS (5.*C18*TI)+5 . *C18*SIN (5 .*C18*TI) )+ 

6 95&  ( (H*C11- (Dll*A/OK) )/( (5 . *C18) **2+Al**2) ) *(A1* 

700&SIN (5 . *C18*TI) -5 . *C18*C0S (5 . *C18*TI) ) 

705  R6= ( (H*C12- (D12*A/0K) ) / ( (6 . *C18) **2+Al**2) ) * (Al* 

710&COS (6 ,*C18*TI)+6 .*C18*SIN (6 .*C18*TI) )+ 

715  & ( (H*C13- (D13*A/0K) ) /( (6 . *C18) **2+Al**2) ) * (Al* 

720&SIN (6 . *C18*TI) -6 . *C18*C0S (6 . *C18*TI) ) 

725  R7=((H*C14-(D14*A/0K) )/( (7.*C18) **2+Al**2) ) * (Al* 

730&COS (7.*C18*TI)+7.*C18*SIN(7. *C18*TI) ) + 

735 & ( (H*C15- (B15*A/0K) ) /( (7. *C18) **2+Al**2) ) * (Al* 

740&SIN (7 . *C18*TI) -7. *C18*C0S (7 . *C18*TI) ) 

745  R8=(H*C16/( (8 .*C18) **2+Al**2) ) * (Al* 

750 SCO S (8.*C18*TI)+8.*C18*SIN(8.*C18*TI) ) + 

755& (H*C17/ ( ( 8 . *C18) **2+Al**2) ) * (Al* 

760&SIN (8 . *C18*TI) -8 . *C18*COS ( 8 . *C18*TI) ) 

761  (2 . *D1*E1) / (Al/A) - (H*C1- (Dl*A/OK) ) /Al- 

76 2 & (H*C2- (D2*A/OK) ) *A1/ ( (C18**2 ) +A1**2 ) *C18* (H*C3- 
763& (A*D3/OK) ) / ( (C18**2 ) +A1**2 ) - (H*C4- (A*D4/OK) ) *A1/ ( (2 . * 
764&C18) **2+Al**2) +2 . *C18* (H*C5- (A*D5/OK) ) /((2.*C18)**2+ 

765  &A1**2 ) -Al* (H*C6- (A*D6/OK) )/( (3 . *C18) **2+Al**2) +3 .*C18* 

766  & (H*C7- (A*D7/OK) )/ ( (3.*C18) * *2+Al**2) -Al* (H*C8- (A*D8/OK) ) 
76  7&/ ( (4 . *C18 ) **2+Al**2 ) +4 . *C18* (E*C9- (A*D9/OK) )/ ( (4 . *C18) ** 
768&2+Al**2 ) -Al* (H*C10- (A-DIO/OK) ) /( (5 . *C18) **2+Al**2) +5 . *C18 
769&* (H  *C11- (A*D11/0K) ) / ( (5 . *C18) **2+Al**2) -Al* (H*C12- (A*D12 
770&/OK) )/( ( 6 . *C18 ) **2+Al**2 ) +6 .*C18* (H*C13- (A*D13/0K) )/( (6 . 
771&*C18) **2+Al**2) -Al* (H*C14- (A*D14/OK) ) / ( ( 7 . *C18) **2+Al**2) 
772S+7 . *C18* (H*C15- (A*D15/OK) ) / ( (7 . *C18) **2+Al**2) -A1*H*C16/ 
773&((8.*C18) **2+Al**2 ) +8 . *C18*H*C17/ ( (8 . *C18) **2+Al**2) 

774  TBAR=G*BB+ 

775&R1+R2+R3+R4+R5+R6+R7+R8 
780  P=TBAR*COS (W*PI*Z/OL) +P 
790  IF  (W  .GE.  ON)  GO  TO  50 
795  W=W+2 . 

800  GO  TO  40 

805  50  P=(2./OL) *P 

810  T=El*OL**2/3.+E2*OL/2 .+E3+P 

815  WRITE (06 ,1000)  TI,Z,T 

820  1000  FORMAT ( 1 TI= 1 F10 . 1 , ' Z= ' F10 . 6 , ' T= ' F10 . 4 ) 

821  IF ( Z .GT.  OL)  GO  TO  825 

822  Z=Z+ . 25 

823  GO  TO  370 
825  825  STOP 
830  END 


APPENDIX  C 


COMPUTATION  OF  BULK  WATER  TEMPERATURE 
ANNUAL  TIME  FRAME 

C_.  PROGRAM  COMPUTES  AVERAGE  MONTHLY  TEMPERATURE 
C 1 REFERS  TO  UNCOVERED  POOL 
C 2 REFERS  TO  SINGLE  SEARS  COVER 

C 3 REFERS  TO  DOUBLE  SEARS  COVER 

C 4 REFERS  TO  SINGLE  BUBBLE  COVER 

C 5 REFERS  TO  DOUBLE  BUBBLE  COVER 

C 6 REFERS  TO  SINGLE  SEARS  COVER  AND  TEN  SOLAR  HEATERS 
DIMENSION  TW(6),  TWA(6),  TWA(6),  TWB(6),  Z(6),  W(6), 
T(6)  , C (6) 

DIMENSION  AW ( 6 ) , FKW(6),  PW(6),  BW(6),  PWS(6), 

P RW ( 6 ) , TAVG ( 6 ) 

DIMENSION  DTW  ( 6 ) , DTA  (6  ) , DE  LT  ( 6 ) , V(6),  HW(6),  UT(6), 
US  (6) 

DIMENSION  QRAD  (6 ) , QEV(6),  QIV(6),  G(6),  H(6),  PA(6) 
DIMENSION  QCT ( 6 ) , QCS(6),  QCDN(6) 

10  READ  (5,109)TW(1)  ,TW(2) 

READ  (5 ,100) APA,TSP,TSSP,TBP,TBBP,HRS 
109  FORMAT  (2F10.5) 

TW  ( 3)  =TW  ( 2 ) 

TW ( 4 ) =TW ( 2 ) 
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TW  (5)  =TW  (2 ) 

QEV  (2)  =0 . 0 
QEV  ( 3)  =0 . 0 
QEV (4) =0 .0 
QEV  (5)  =0.0 

READ  (5 ,101)Q,TA,HO,QEV(1) 

101  FORMAT  (3F10.5 ,E10.2) 

WRITE  (6 ,101) Q ,TA ,HO ,QEV ( 1) 

WRITE  (6 , 105 ) APA,TS? ,TSSP ,TBP ,TBBP , HRS 
QEV (1) =QEV (1) /30 . 

WRITE  (6,102)  (TW(I)  ,T=1;5) 

Q-Q  *3.68 

FKA=0. 014- (32. 0*0 .0014) /6 8.0+ (0 .0014/68.0) *TA 
PRA=0 .72 

AA= (3.16+32.0*(1.4/68.0)-(l. 4/68 . 0 ) *TA) *1 . 0E6 
AT=176 .0 
AS=15  3 . 0 
TWA ( 1 ) =TW ( 1 ) 

TWB (1) - TW (1) 

TWA ( 2 ) =TW ( 2 ) 

TWB ( 2 ) =TW ( 2 ) 

TWA ( 3 ) =TW ( 3 ) 


TWB  (3)  =TW(3) 


TWA ( 4 ) =TW ( 4 ) 

TWB  (4)  =TW(4) 

TWA  ( 5 ) =TW  ( 5 ) 

TWB ( 5 ) =TW ( 5 ) 

DO  26  N=1 , 5 
TAVG  (N)  =TW  (N) 

26  CONTINUE 
11  DO  22  1=1,5 

M= (TWB (I) -30 . 0 ) /10 . 0 

GO  TO  (1,2, 3,4,5, 6,7, 8, 9)  ,M 

1 AW(I) =-20 .5E6  + 0 . 5 7E6  *TWB ( I ) 
FKW (I) =0 .297  + 0 . 0007*TWB (I) 
PW(I) =-0 .1039  + 0.00564*  TWB(I) 
3W  (I)  =-0 . 96E-4+0 .29  E-5  *TWB(I) 
PWS (I) =-0 .1039+0 .00564*  TWC 
PRW(I)  = 19 .8-0 .205*TWB (I) 

GO  TO  22 

2 AW (I) =-44 .0  E6+1.04E6*  TWB(I) 

FKW ( I ) = 0. 292  + 3. E -4  *TW3 (I) 
PW(I)=-0 .2129  + 0.00782  * TWB(I) 
BW ( I ) =-l . 31E-4  + 0.36E-5  *TWB(I) 
PWS (I)=-0. 2129  + 0.00782  *TWC 
PRW (I) = 17.15-0. 152*TWB (I) 
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GO  TO  22 

3 AW (I) =-78 .9  E6+1.62E6  * TWB(I) 

PW(I) =-0 .38  45+0 .0106  3*TWB (I) 

PWS(I)=-0. 3845+0. 01068  * TWC 

FKW  (I)=0.277  + 13.5E-4  *TWB  ( I ) -5  . OE-6  * (TWB  (I)  **2.0) 
PRW(I) =21 . 19-0 . 303*TWB (I) +1 . 4E-3* (TWB (I) **2.0) 

BW(I) =-2 . 3E-4+6 . 75E-6  *TWB (I) -2 .5E-B* (TWB (I) **2.0) 

GO  TO  22 

4 AW (I) =-115 . 4E6  + 2 . 14E6  *TWB ( I ) 

PW(I) =-0 . 6435  + 0.01438  * TWB(I) 

PWS (I) =-0 .6435  + 0.01438  * TWC 

FKW(I)=0.277  + 13 .5E-4*TWB (I) -5 . OE-6* (TWB (I) **2 . 0 ) 
PRW(I) =21 .19-0 . 303*TWB(I) +1.4E-3* (TWB (I) **2.0) 

BW (I) =-2 . 3E-4+6 . 75  E-6*TWB (I) -2 .5E-8* (TWB (I) **2.) 

GO  TO  22. 

5 AW(I) =-176 .E6  +2.9E6  * TWB(I) 

PW(I) =-l .0235+0 . 01918*TWB (I) 

PWS (I) =-l . 0235+0 . 01913*TWC 
FKW(I) =0 .305+6 .E-4 *TWB (I) 

BW (I ) =-0 .9E-4+0 . 3E-5*TWB (I) 

PEW (I) =11. 9 7-0. 976 *TWB (I) 

GO  TO  22 


6 


AW(I)=-212.E6  + 3.3  E6  * TWB(I) 


PW(I) =-l .5603  +0 . 025 1*TWB ( I ) 
■PWS  (I)  =-l  .5608+0 . 0251*  TWC 
FKV? ( I ) = 0.314+  5 .E-4  *TWB(I) 
BW ( I ) = 0 . 2E-5  *TWB ( I ) 

PRW (I) =10 .62-0 . 061  * TWB(I) 
GO  TO  22 

PW(I)  =-2.306  8+- . 0 3256  *TWB (I ) 
PWS (I)=-2 .3068+0 .03256*TWC 
FKW(I) =0 . 321+4 . 3E-4  *TWB (I) 

PRW (I) =7 .52-0 . 03*TWB (I) 

AW  (I)  =-292  .E6  + 4 . 1E6*TWB  (I) 

BW (I) =-l .E-4  + 3.E-6  *TWB (I) 

GO  TO  22 

PW(I)  =-3 . 3188+0 .0+176 *TWB (I ) 
PWS (I)=-3. 3188+0 . 04176*TWC 
FKW(I) =0 .3221+4 .2E-4*TWB (I) 
PRW (I) =8 .62-0 . 04*TWB ( I) 

AW (I) =-518 .E6+6 . IE 6 * TWB (I) 

BW ( I ) =-0 . 45E-4+2 ,5E-6*TWB (I) 
GO  TO  22 

PW(I)  =-4 .6688+0 .05  30 1*TWB ( I ) 
PWS (I) =-4 .6688+0 .053301*TWC 
FKW(I) =0 . 325  7+3 . 9E-4  *TWB ( I ) 


PRW(I) =8 .62-0 . 0 4 *TWB ( I ) 

AW (I) =-620 .E6  + 7 .E6  *TWB (I) 

BW (I) =9 .E-6  + 2.2E-6  *TWB (I) 

CONTINUE 
DO  23  J=l,  5 

Z ( J) = ( (FKA/FKW ( J) * ( (AA*PRA) / (AW ( J) *PRW ( J) ) ) **0 . 25 ) **0 . 8 
W(J)=Z  (J)+l. 

T (J)  =TW(J)  -TA 

DTW  ( J ) = ( Z (J)  *T(J)  )/W(J) 

DTA  ( J)  =DTW  ( J)  /Z  (J) 

DTW  ( J)  =ABS  (DTW  ( J)  ) 

V ( J)  = (32 . 2*BW(  J)  * (DTW  ( J)  /2.0)*4.0)**0.5 

LW( J)  = ( (C . 55  *FKW ( J) ) /2 . 0) * ( ( AW ( J ) *PRW (J) * 8 . 0 * DTW ( J ) ) ** 
0.25) 

UP  IS  NEGATIVE  AND  DOWN  IS  POSITIVE 
IF  (T  ( J)  .GE.  0.0)  GO  TO  23 
V(J)=-V(J) 

CONTINUE 

UT (1) =1 .0  /((  1 . 0/HW (1) ) + ( 1 . 0/HO) ) 

UT (2 ) =1 . / ( ( 1 . /HW ( 2 ) ) + (1 ./HO) +.052) 

UT ( 3 ) = 1.0/  ((1. 0/HW (3))  + ( 1 . 0/HO) +0.864) 

UT ( 4 ) =1 .0/ ( ( 1 . 0/HW ( 4 ) ) + ( 1 . 0/HO) +0.864) 

UT (5) =1 . 0/ ( ( 1 . 0/HW (5 ) ) + (1 . 0/HO)  + 0,978) 
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DO  24  K=1 , 5 

US (K) =1 . 0/ ( (1 . 0/HW (K) ) + (1 . O/HO) +0.0524) 

QRAD(K)=0 .0 
24  CONTINUE 

QIN ( 1 ) =Q  * AT  *APA 
QIN (2) =Q*AT*TSP* . 85 
QIN (3) =Q*AT*TSSP* . 85 
QIN ( 4 ) =Q  * AT  *TBP  * . 8 5 
QIN (5 ) = Q*AT*  TBBP* . 85 
DO  28  JJ=1 , 5 

QRAD (JJ) =AT* . 1714E-8E . 96* ( (TAVG (JJ)+459.7)**4.- 
1 (TA+459. 7) **4 . (* (24 .-HRS) 

TAVG  (JJ)=TA+(  (QIN  (JJ)  -QEV(JJ)  -QRAD(JJ)  )/(24.*  (AT 
1*  UT (JJ)+AS*US (JJ)  ) ) ) 

QCS (JJ)=US (JJ) *AS* (TAVG ( JJ) -TA) *HRS 
QCT  ( JJ)  -UT  ( jj)  *AT*  (TAVG  ( JJ)  -TA)  *HRS 
QCON (JJ) =QCT (JJ) +QCS (JJ) 

28  CONTINUE 

DO  27  11=1,  5 

C ( II) =ABS (TWA ( II ) -TAVG ( II ) ) 

IF  (C(l'I)  .LE.  0.1)  GO  TO  2 7 
TWB(II)  = (TW(II)  +TAVG  (II)  ) /2  . 

TWA(II)=TAVG(II) 
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GO  TO  11 
27  • ' CONTINUE 
G (I) =0 . 0 

WRITE  (6,111)  (QCT(KK)  ,QCS(KK)  ,QCON(KK)  , KK=1,5) 

WRITE  (6 ,104)  (QIN (L)  , QEV(L)  ,QRAD(L)  ,TAVG(L)  , L=l,5) 
WRITE  (6 ,106) T (1) , DTW (1) ,DTA(1) , DELT (1) ,V(1) 

WRITE  (6 ,108)  T (2)  ,DTW(2)  ,DTA(2)  , DE  LT  ( 2 ) ,V(2) 

WRITE  (6 ,108)  T (3)  ,DTW(3)  ,DTA(3)  , DELT  (3)  ,V(3) 

WRITE  (6 , 108 ) T (4 ) ,DTW(4)  ,DTA(4)  , DELT  ( 4 ) ,V(4) 

WRITE  (6 ,108) T (5)  ,DTW(5)  ,DTA(5)  , DE  LT ( 5 ) ,V(5) 

WRITE  (6 ,10  8)  HW  ( 1)  ,UT(1)  ,US(1)  ,G(1) 

WRITE  (6 ,107)HW(2)  ,UT  (2)  ,US  (2) 

WRITE  (6,107)  HW  (3)  ,UT(3)  , US  ( 3 ) 

WRITE  (6,107)  HW  ( 4 ) ,UT(4)  ,US(4) 

WRITE  (6 ,107)HW(5) ,UT(5) ,US (5) 

100  FORMAT  (6F10.5) 

102  FORMAT  (1X,5F20.5) 

103  FORMAT  (//4F20.5) 

104  FORMAT  (/////3E2C . 6 .F2 0 . 6 ) 

105  FORMAT  (6F15.5) 

106  FORMAT  (//3F20 . 5 , 2E20 . 6) 

107  FORMAT  (3F20.5) 

108  FORMAT  ( 3F2 0 . 5 , 2E 20 . 6 ) 
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111  FORMAT  (3E20.6) 
GO  TO  10 


END 


APPENDIX  D 


TIME  VARIATION  OF  AVAILABLE  SOLAR  ENERGY  (H  ) 

su 


WITH  TIME 


Date 

<J> 

Solar 

Time 

I 

dn 

I 

dn 

X 

0 .87* 

I 

dir 

I 

dif 

HSU  - (Iair 

+ W 

1/1 

78.5' 

8 

4 

202 

176 

34 

0 

2x34  = 

68 

68.5 

9 

3 

262 

228 

84 

16 

2x100  = 

200 

60.5 

10 

2 

287 

250 

123 

20 

2x143  = 

286 

55 

11 

1 

299 

261 

149 

21 

2x170  = 

340 

53 

12 

12 

304 

264 

159 

22 

1x181  = 

181 

1075 

2/1 

75 

8 

4 

230 

200 

52 

14 

2x66  = 132 

64.5 

9 

3 

273 

237 

104 

19 

2x123  = 246 

55.5 

10 

2 

29  7 

258 

146 

22 

2x168  =336 

49.5 

11 

1 

307 

267 

173 

23 

2x196  = 392 

47.5 

12 

12 

310 

270 

183 

24 

1x207  = 207 
1313 

3/1 

69 

57.5 
48 

40.5 
38 

8 

9 

10 
11 
12 

4 

3 

2 

1 

12 

252 

287 

304 

313 

316 

xO . 90  * 
226 
258 
274 
282 
284 

81 

138 

183 

214 

224 

18 

22 

24 

26 

26 

2x99 

2x160 

2x207 

2x240 

1x250 

= 

198 

320 

414 

480 

250 

1662 

4/1 

75 

7 

5 

205 

184 

48 

17 

2x65 

130 

62 

8 

4 

262 

236 

101 

23 

2x124 

= 

248 

49.5 

9 

3 

285 

256 

166 

26 

2x192 

= 

384 

38 

10 

2 

300 

270 

213 

28 

2x241 

= 

482 

29 

11 

1 

307 

276 

242 

29 

2x271 

=3 

542 

25.5 

12 

12 

312 

281 

254 

29 

1x283 

= 

283 

2069 

*The 

values  of 

0.87  and 

0.90  are 

clearness 

factors  as 

reported 

in  [37]. 
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Appendix  D continued 


Date 

<f> 

Solar 

Time 

I 

dn 

T 

dn 

X 

0.90* 

I 

dir 

I 

dif 

H = 
su  dir 

+ I . ) 
dir 

5/1 

70 

7 5 

210 

' 189 

65 

22 

2x87  = 

174 

57 

8 4 

252 

227 

124 

26 

2x150  = 

300 

44 

9 3 

275 

247 

178 

29 

2x207  = 

414 

31.5 

10  2 

287 

258 

220 

31 

2x251  = 

502 

20.5 

11  1 

295 

266 

249 

31 

2x280  = 

560 

15 

12  12 

298 

268 

259 

32 

1x291  = 

291 

2241 

6/1 

67 

54 

41 

28 

16 

8 

7 5 

8 4 

9 3 

10  2 
11  1 
12  12 

209 

244 

266 

277 

284 

288 

188 

220 

240 

249 

256 

259 

70 

130 

181 

220 

246 

256 

24 

28 

31 

32 

33 
33 

2x94  = 
2x158  = 
2x212  = 
2x252  = 
2x279  = 
1x289  = 

= 188 
= 316 
= 424 
= 504 
= 558 
= 289 
2279 

7/1 

66.5 

7 5 

205 

184 

73 

25 

2x148  = 

= 296 

53.5 

8 4 

240 

216 

128 

29 

2x157  = 

= 314 

40.5 

9 3 

260 

234 

178 

32 

2x210  = 

= 420 

27 

10  2 

273 

246 

219 

33 

2x252  = 

= 504 

15 

11  1 

278 

250 

242 

33 

2x275  = 

= 550 

7 

12  12 

282 

254 

252 

34 

1x286  = 

= 286 

2370 

8/1 

68.5 

7 

5 

197 

177 

65 

24 

2x89  =178 

55.5 

8 

4 

237 

213 

121 

28 

2x149  = 298 

42.5 

9 

3 

260 

234 

172 

31 

2x203  = 406 

30 

10 

2 

273 

246 

213 

32 

2x245  = 490 

IS 

11 

1 

278 

250 

238 

33 

2x271  = 542 

12 

12 

12 

282 

254 

248 

33 

1x281  = 281 
2195 

*The 

values  of 

0.87  and 

0.90 

are  clearness 

factors  as  reported 

in  [37] 
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Appendix  D continued 


Date 

<P 

Solar 

Time 

c 

TD 

M 

I 

dn 

I 

dir 

I 

dif 

H 

su 

<1  . 
dir 

X 

0.90* 

+ I . ) 
dif 

9/1 

73 

7 

5 

188 

169 

49 

21 

2x70 

= 140 

60 

8 

4 

240 

216 

108 

25 

2x133 

= 266 

47 

9 

3 

264 

238 

162 

28 

2x190 

= 380 

35 

10 

2 

277 

249 

204 

31 

2x235 

= 470 

25.5 

11 

1 

284 

256 

2 31 

32 

2x263 

= 526 

21 .5 

12 

12 

288 

259 

241 

32 

1x273 

= 273 

2055 

xO  .87 

10/1 

78.5 

7 

5 

173 

151 

30 

14 

2x44 

= 88 

66 

8 

4 

237 

206 

84 

20 

2x104 

= 208 

54 

9 

3 

267 

232 

136 

25 

2x161 

= 322 

44 

10 

2 

285 

248 

179 

28 

2x207 

= 414 

36 

11 

1 

294 

256 

207 

29 

2x236 

= 472 

29 

12 

12 

297 

258 

226 

29 

1x255 

= 255 

1759 

H/1 

72.5 
62 

52.5 

46.5 
44 

8 4 

9 3 

10  2 
11  1 
12  12 

226 

266 

286 

297 

300 

197 

232 

248 

258 

261 

59 

109 

151 

178 

216 

15 

20 

24 

26 

26 

2x74  = 148 

2x129  = 258 
2x175  = 350 
2x204  = 408 
1x242  = 242 
1406 

12/1 

77.5 

8 4 

207 

180 

39 

2x39  = 78 

67.5 

9 3 

263 

229 

88 

17 

2x105  = 210 

59 

10  2 

287 

250 

129 

21 

2x150  = 300 

53.5 

11  1 

297 

258 

154 

22 

2x176  =352 

51.5 

12  12 

302 

262 

16  3 

23 

1x186  = 186 

1126 

*The  values  of  0.87  and-  0.90  are  clearness  factors  as  reported 


in  [37] 


APPENDIX  E 


APPARENT  TRANSMISSIVITIES  OF  ABSORPTIVITY  OF  WATER, 
SINGLE  SEARS  PLASTIC,  AND  DOUBLE  SEARS  PLASTIC 


Date 

Solar 

Time 

H 

su 

P 

a 

a*H 

su 

a 

T 

sp 

T 

sp 
x H 

si 

T 

sp 

T 

ssp 

T 

ssp 
x H 

su 

T 

ssp 

1/1 

8 4 

68 

0.72 

49 

0.31 

21 

0.46 

31 

9 3 

200 

0.82 

164 

0.71 

142 

0.635 

127 

10  2 

286 

0.86 

246 

0.855 

0.80 

229 

0.767 

0.715 

204 

0 .695 

11  1 

340 

0 . 88 

300 

0.83 

282 

0 .735 

250 

12  12 

181 

0.88 

159 

0.835 

151 

0.742 

134 

1075 

918 

825 

746 

2/1 

8 4 

132 

0.77 

101 

0.535 

70 

0.52 

68 

9 3 

246 

0.83 

204 

0.76  5 

188 

0.68 

167 

10  2 

336 

0.88 

296 

0.865 

0.83 

279 

0.795 

0.74 

248 

0.71 

11  1 

392 

0.39 

349 

0.85 

333 

0.75 

294 

12  12 

207 

0.90 

136 

0.85 

176 

0.755 

156 

1313 

1136 

1046 

933 

3/1 

8 4 

198 

0.82 

162 

0 .705 

140 

0 .635 

126 

9 3 

320 

0.87 

273 

0.82 

262 

0.73 

234 

10  2 

414 

0.89 

368 

0.89 

0.85 

352 

0.829 

0.75 

311 

0.742 

11  1 

4 80 

0.92 

442 

0.855 

411 

0.765 

367 

12  12 

250 

0.92 

230 

0.86 

215 

0.77 

193 

1662 

1480 

1330 

1231 

VI 

7 5 

130 

0.77 

100 

0.535 

69 

0.52 

68 

8 4 

248 

0.85 

211 

0.79 

196 

0.705 

175 

9 3 

384 

0.89 

342 

0.85 

326 

0.75 

288 

10  2 

482 

C .92 

443 

0.902 

0.86 

414 

0.329 

0.77 

371 

0 .747 

11  1 

542 

0.93 

503 

0.86 

466 

0.775 

420 

12  12 

283 

0.9  35 

264 

0.865 

24  5 

0.78 

221 

2069 

1863 

1716 

1543 

5/1 

7 5 

174 

0.81 

141 

0.69 

120 

0.625 

109 

8 4 

300 

0.87 

261 

0.82 

246 

0.73 

219 

9 3 

414 

0.905 

375 

0.85 

352 

0.75 

315 

10  2 

502 

0.925 

464 

0.912 

0.86 

432 

0.842 

0.775 

389 

0.757 

11  1 

560 

0.94 

527 

0.865 

484 

0.78 

437 

12  12 

291 

0.94 

274 

0.865 

252 

0.78 

227 

2241 

2042 

1886 

6/1 

7 5 

188 

i 

0.83 

156 

0.735 

138 

0 .66 

124 

8 4 

316 

0.88 

278 

0.83 

262 

0.74 

234 

9 3 

424 

0.91 

386 

0.916 

0.855 

36  2 

0.847 

0.765 

324 

0.762 

10  2 

504 

0.93 

469 

0.86 

433 

0.775 

390 

11  1 

558 

0.94 

525 

0.865 

483 

0.78 

436 

12  12 

289 

0.94 

272 

0.865 

250 

0.78 

226 

2279 

2086 

1928 

1734 
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Date 

Solar 

Time 

H 

su 

a 

a-H 

su 

a 

T 

sp 

I 

sp 
x H 

su 

T 

'sp 

T 

ssp 

T 

ssp 
x H 

su 

Tssp 

7/1 

7 

5 

296 

0.83 

246 

0.745 

221 

0.67 

198 

8 

4 

314 

0.88 

276 

0.835 

262 

0.74 

232  ' 

9 

3 

420 

0.91 

382 

0.911 

0.855 

359 

0.843 

0.765 

321 

0.758 

10 

2 

504 

0.93 

469 

0.86 

433 

0.775 

391 

11 

1 

550 

0.94 

517 

0.865 

476 

0.78 

4 29 

12 

12 

• 286 

0.94 

269 

- 

0.865 

24  7 

0.78 

223 

2370 

2159 

1998 

1792 

8/1 

7 

5 

178 

0.81 

144 

0.71 

126 

0 .635 

113 

8 

4 

29  8 

0.37 

259 

0.85 

247 

0.735 

219 

9 

3 

406 

0.905 

36  7 

0.908 

0.855 

346 

0.844 

0.765 

310 

0.758 

10 

2 

490 

0.9  25 

453 

0.86 

422 

0.775 

380 

11 

1 

542 

0.94 

508 

0.865 

468 

0.78 

423 

12 

12 

281 

0.94 

264 

0.865 

243 

0.78" 

219 

2195 

1995 

1852 

1664 

VO 

\ 

M 

7 

5 

140 

0.775 

109 

0.62 

87 

0.57 

80 

8 

4 

266 

0.86 

228 

0.805 

214 

0.715 

190 

9 

3 

380 

0.90 

342 

0.90  2 

0.85 

323 

0.75 

285 

10 

2 

470 

0.92 

433 

0.86 

404 

0.835 

0.77 

36  2 

0.750 

11 

1 

526 

0.93 

488 

0.86 

453 

0.78 

411 

12 

12 

273 

0.94 

256 

0.865 

236 

0.78 

213 

2055 

1856 

1717 

1541 

10/1 

7 

5 

88 

0.72 

63 

8 

4 

208 

0.83 

173 

0.75 

156 

0.67 

139 

9 

3 

322 

0.37 

280 

0.83 

267 

0.74 

238 

10 

2 

414 

0.905 

375 

0.888 

0.85 

352 

0.837 

0.76 

315 

0.752 

11 

1 

472 

0.92 

434 

0.86 

406 

0.77 

364 

12 

12 

255 

C .925 

236 

0.86 

219 

0.775 

198 

1759 

1561 

1400 

1254 

11/1 

8 

4 

148 

0.77 

114 

0.63 

93 

0.58 

86 

9 

3 

258 

0.85 

219 

0.79 

204 

0.705 

182 

10 

2 

350 

0.88 

308 

0.873 

0.84 

294 

0.814 

0.745 

260 

0.726 

11 

1 

408 

0.90 

367 

0.85 

347 

0.755 

308 

12 

12 

242 

0.90 

218 

0.85 

206 

0.76 

184 

1406 

1226 

1144 

1020 

12/1 

8 

4 

78 

0.73 

57 

0.31 

24 

0.46 

36 

9 

3 

210 

0.82 

172 

0.73 

153 

0 .65 

136 

10 

2 

300 

0.865 

259 

0.855 

0.81 

24  3 

0.772 

0.72 

216 

0.70 

11 

1 

352 

0.88 

310 

0.83 

29  2 

0.74 

260 

12 

12 

186 

0.88 

164 

0.845 

157 

0.745 

139 

962 

86  9 

787 

APPENDIX  F 


APPARENT  TRANSMISSIVITIES  OF  BUBBLE  PLASTIC  COVER 


Date 

Solar 

Tine 

H 

su 

Tbp 

T.  -H 
bp  su 

Tbp 

Tbbp 

T,  , *H 
bbp  su 

Tbbp 

7/1 

7 5 

296 

0.655 

194 

0.548 

162 

8 4 

314 

0.728 

229 

0.618 

194 

9 3 

4 20 

0.768 

322 

0.762 

0.658 

276 

0 .656 

10  2 

504 

0.786 

396 

0.68 

343 

11  1 

550 

0.795 

437 

0.69 

380 

12  12 

286 

0.79  8 

'228 

0.692 

198 

2370 

1806 

1553 

8/1 

7 5 

178 

0.615 

109 

0.51 

91 

8 4 

298 

0.72 

214 

0.61 

182 

9 3 

406 

0.76  2 

309 

0.761 

0.653 

265 

0 .654 

10  2 

490 

0.783 

384 

0.676 

331 

11  1 

542 

0.793 

429 

0.687 

372 

12  12 

281 

0.797 

224 

0.69 

194 

2195 

1669 

14  35 

9/1 

7 5 

140 

0.47 

66 

0.45 

63 

8 4 

266 

0.695 

185 

0.587 

156 

9 3 

380 

0.75 

285 

‘ 0.746 

0.64 

243 

0.643 

10  2 

470 

0.778 

366 

0.67 

315 

11  1 

526 

0.79 

415 

0.68 

358 

• 12  12 

273 

0.791 

216 

0.685 

187 

2055 

15  33 

1322 

10/1 

7 5 

88 

0.2 

17 

0.2 

17 

8 4 

208 

0.645 

134 

0.538 

112 

9 3 

322 

0.727 

234 

0.722 

0.617 

199 

0.617 

10  2 

414 

0.76 

315 

0.65 

26  9 

11  1 

472 

0.775 

366 

0.667 

315 

12  12 

255 

0.786 

205 

0.677 

173 

1759 

1271 

1085 

11/1 

8 4 

148 

0.47 

70 

0.45 

67 

9 3 

258 

0.681 

176 

0.572 

147 

10  2 

350 

0.732 

256 

0.706 

0.6  22 

218 

11  1 

408 

0.752 

30  7 

0.642 

26  2 

12  12 

242 

0.76 

184 

0.65 

157 

1406 

993 

851 

12/1 

8 4 

78 

0.2 

15 

0.2 

15 

9 3 

210 

0.63 

132 

0.525 

110 

10  2 

300 

0.70 

210 

0.667 

0.592 

177 

0.566 

11  1 

352 

0.728 

256 

0.618 

218 

12  12 

186 

0.736 

137 

0.626 

116 

1126 

750 

636 

216 
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Date 

Solar 

Time 

H 

su 

Tbp 

T.  #H 
bp  su 

Tbp 

Tbbp 

T.  , -H 
bbp  su 

Tbbp 

1/1 

8 4 

68 

0.2 

13.6 

0.2 

13.6 

9 3 

200 

0.615 

123 

0.51 

10  2 

10  2 

286 

0.695 

199 

0.663 

0.585 

167 

0.561 

11  1 

340 

0.722 

245 

0.613 

208 

12  12 

181 

0.732 

132 

0.622 

112 

1075 

713 

603 

2/1 

8 4 

132 

0.5 

66 

0.4 

53 

9 3 

246 

0.66 

16  2 

0.551 

135 

10  2 

336 

0.72 

242 

0.696 

0.61 

205 

0.588 

11  1 

392 

0.742 

290 

0.632 

248 

12  12 

207 

0.75 

155 

0.64 

132 

1313 

915 

773 

3/1 

8 

9 

10 
11 
12 

4 

3 

2 

1 

12 

198 
320 
414 
4 80 
250 
1662 

0.61 

0.71 

0.748 

0.768 

0.773 

121 

227 

310 

369 

193 

1220 

0.734 

0.508 

0.60 

0.638 

0.658 

0.663 

101 

192 

264 

316 

166 

1039 

0.624 

4/1 

7 

5 

130' 

0.5 

65 

0.4 

52 

8 

4 

248 

0.681 

169 

0.572 

142 

9 

3 

384 

0.742 

285 

0.632 

242 

10 

2 

482 

0.773 

372 

0.744 

0.653 

319 

0.636 

11 

1 

542 

0.786 

426 

0.677 

36  7 

12 

12 

283 

0.79 

224 

0.68 

192 

2069 

1541 

1314 

5/1 

7 

8 

9 

10 
11 
12 

5 

4 

3 

2 

1 

12 

174 

300 

414 

502 

560 

291 

2241 

0.595 

0.71 

0.76 

0.782 

0.792 

0.795 

103 

213 

314 

392 

443 

231 

1696 

0.757 

0.50 

0.60 

0.65 

0.675 

0.685 

0.69 

87 

180 

259 

339 

384 

201 

1460 

0.651 

6/1 

7 

5 

188 

0.63 

118 

0.53 

100 

8 

4 

316 

0.727 

230 

0.617 

195 

9 

3 

424 

0.766 

325 

0.761 

9.657 

279 

0.658 

10 

2 

504 

0.787 

397 

0.678 

342 

11 

1 

558 

0.795 

444 

0.69 

385 

12 

12 

289 

0.79  8 

231 

0.692 

200 

2279 


1745 


1501 


APPENDIX  G 


A SAMPLE  OF  MEASURED  WATER  TEMPERATURES 
FOR  UNCOVERED  POOL; 

DATA  WAS  TAKEN  ON  12/1/72 


8 Equal  Spaces  at  18  In.  Each 
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■n 

<o 

w 


5 

a 

o 

w 


Measured  Water  Temperatures  (F) 
Time:  0600  Hours 


Equal  Spaces  at  18  In.  Each 
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in 


Measured  Water  Temperatures  (F) 
Time:  1200  Hours 


APPENDIX  H 


SOIL  TEMPERATURES  AND  AMBIENT  AIR 
VS.  TIME,  REFERENCE  [36] 


Surface 
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(£)  san^raaduiai 
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